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The effects of random geometric imperfections on the transmission of the 
TEoi wave in circular waveguide are studied; the necessary theory of guides 
with known arbitrary imperfections is first developed. The TE i transmis- 
sion statistics are determined in terms of the statistics of the various types 
of geometric imperfections. Both discrete mode converters — i.e., localized 
imperfections such as tilts, offsets, or diameter changes at joints between 
pipes that are perfect right-circular cylinders — and continuous geometric 
imperfections — such as straight ness deviation, diameter variation, ellip- 
ticily, etc., that vary smoothly with distance along the guide — are con- 
sidered. The average, variance, power spectrum, and probability distribu- 
tion of the TEoi loss-frequency curve are discussed. 

Continuous straightness deviation (of the individual pipes of the guide) 
appears to be the most serious tolerance in present copper waveguide, and 
a significant factor in helix guide as well. The power spectrum of the straight- 
ness deviation is all-important in determining the TE m loss due to mode 
conversion. Fourier components of straightness deviation having wave- 
lengths between roughly 1.4 and 4,4 feet arc the significant ones for the 
present 2-inch I.D. guide operated in a frequency band from 35 to 90 kmc. 
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I. INTRODUCTION 

Long distance waveguide transmission via the TE i mode in circular 
waveguide is an attractive goal because the theoretical TE i heat loss 
decreases monotonically as the operating frequency increases. As is 
well known, operating frequencies far above the TEoi cutoff frequency 
are required to realize sufficiently low heat loss and delay distortion, 
so that the guide must operate far into the multimode region. Thus, 
considering a typical case, a 2-inch I.D. perfect copper circular guide 
operating at 55 kmc will have a theoretical TEoi heat loss of 1.54 db/ 
mile; but this guide will propagate 223 additional modes, which we 
call spurious modes, at this frequency. 

The TEoi transmission loss will approach the theoretical TE i heat 
loss in a copper waveguide only if the waveguide is a geometrically 
perfect right-circular cylinder over its entire length. Any departure 
from this ideal geometry will couple the TE i mode to some of the 
spurious modes 2 to 7 the net effect of this coupling will be to increase 
the TEoi transmission loss above the theoretical heat loss, and to cause 
the TEoi transmission loss to vary with frequency. ' 

Two types of geometric imperfections are of interest: 

(a) Intentional deformations introduced in the guide for various 
reasons, e.g., to go around corners, 3,6 ,7 ' 9 to taper from a small guide to 
a larger one, 10 etc. Mode conversion effects control the design of such 
devices, but we will not discuss them further. 

(b) Random geometric imperfections arising during the manufac- 
turing or the laying of the guide; these imperfections will increase the 
TEoi loss and cause it to vary randomly with frequency. The study of 
such effects is the purpose of the present paper. 

The transmission characteristics of multimode waveguide with such 
random imperfections may be improved by adding heat loss to the 
spurious modes while keeping the TE i heat loss low, i.e., close to its 
value for ideal copper guide. Examples are helix waveguide 11 and copper 
waveguide with a thin lossy dielectric lining. ' This additional spuri- 
ous mode loss will reduce the TEoi loss fluctuations with frequency, 
but will not reduce the average TEoi loss. (In contrast, for large in- 
tentional bends, it is desirable to alter the phase constant of one par- 
ticular copper guide mode without increasing the heat loss to either 
TEoi or any of the spurious modes. 3,14 ) 
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The present paper is concerned primarily with determining the 
statistics of the TEoi transmission in terms of the statistics of the vari- 
ous geometric imperfections. The results of this analysis will indicate 
the required tolerances on the various types of geometric imperfections 
in different types of guide as a function of the allowable transmission 
degradation. Finally, the computed transmission statistics will be use- 
ful in determining the over-all degradation in various possible communi- 
cation systems using imperfect waveguide as a transmission medium. 

II. THEORY OF GUIDES WITH KNOWN IMPERFECTIONS 

In this section we summarize the theory of circular waveguide with 
known geometric imperfections. These results yield, at least in princi- 
ple, the TEoi transfer function for a circular waveguide with an arbi- 
trary, known departure from perfect geometry. If the various geometric 
imperfections are assumed to be random processes, then the TEoi trans- 
fer function will also be a random process. In Sections III and IV these 
results for known imperfections are used to determine the TEoi trans- 
mission statistics in terms of the statistics of the various geometric 
imperfections. 

We first require a solution to Maxwell's equations, in terms of the 
normal modes of the guide in question, for boundary conditions given 
by different types of geometric imperfections in various types of guide. 
Several people have studied these problems over the past ten years; 
we will give below a brief description of some of this work. 

Transmission of TEoi through bends was first studied by M. Jou- 
guet, 6 and by S. O. Rice in unpublished work. S. E. Miller made use 
of these results to devise several methods for transmitting TE i around 
intentional bends. 

S. P. Morgan first computed via perturbation theory the first-order 
spurious modes scattered from a unit incident TEoi wave by small, 
abrupt tilts, offsets, and diameter changes in ideal lossless metallic 
guide. 15 These results were derived independently and published by S. 
Iiguchi. 16 Such discontinuities will often be called discrete mode con- 
verters, because the guide possesses perfect cylindrical geometry except 
at isolated, discrete points along its axis. 

Next, Morgan determined the first-order spurious modes scattered 
by an arbitrary small continuous deformation of the surface of an 
ideal lossless metallic guide, again via perturbation theory. These 
results may be used to determine the first-order spurious modes excited 
by small continuous diameter variation, straightness deviation, ellip- 
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ticity, and higher-order cross-sectional deformations of metallic guide. 
In addition, Morgan used these results to evaluate the expected value 
of the additional TE i loss due to mode conversion caused by random 
distortion of the guide for a rather special mathematical model of the 
guide distortion. 

Finally, Morgan applied the generalized telegraphist's equations of 
Schelkunoff 4 to the problem of lossless metallic waveguide with an 
arbitrary curvature of its axis; 2 if desired, the dielectric constant of the 
material filling the guide may be an arbitrary function of position. 
By this analysis, Maxwell's equations are reduced to an infinite set of 
coupled differential equations, the coupled line equations, 17 in which the 
dependent variables are the complex mode amplitudes of the normal 
modes of the unperturbed metallic guide. In principle, the coupled line 
equations provide an exact description of a lossless metallic waveguide 
with an arbitrary straightness deviation, and are not subject to the 
severe restrictions of the perturbation theory which was described in 
the preceding paragraph. Thus, if this infinite set of differential equations 
could be solved for an arbitrary straightness deviation, we would have 
an exact solution for Maxwell's equations for the particular deformed 
guide. As will appear below, solutions to these equations in the general 
case are not available, and useful results are obtained only by applying 
perturbation theory of one form or another to the coupled line equa- 
tions; 18 however, the equations themselves are an exact description of 
the field problem. 

Since the loss in real copper guide is low, we expect that the above 
coupling coefficients, which strictly speaking apply only to lossless 
metallic guide, will provide a good approximation for copper guide, and 
that the coupled line equations for lossless metallic guide need be modi- 
fied only by changing the propagation constants for the various modes 
to take account of the small losses actually present in copper guide. 

Equivalent results have been derived by H. G. Unger for various 
geometric imperfections in helix waveguide via the generalized tele- 
graphist's equations; 4 these analyses have been carried out both in 
terms of metallic guide modes and helix guide modes. 11 Unger has stud- 
ied straightness deviations 3 and cross-sectional deformations 10 ' 20 in 
helix, as well as winding imperfections 21 in helix, and has given propaga- 
tion constants and coupling coefficients for the various cases. 22 He has 
similarly studied continuous diameter variations (tapers), 10 which have 
identical behavior in both helix and metallic guide. 

Thus, the study of random geometric imperfections in copper or helix 
waveguide has been reduced to the study of solutions to an infinite set 
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of differential equations with random coupling coefficients. However, 
this latter problem is a formidable one for which there is no really 
satisfactory treatment except in rather special cases. In the present 
paper, we use perturbation theory to approximate the solution to the 
coupled line equations; however, there is so far no rigorous justification 
for this approach. The convergence of this approximate solution is 
discussed elsewhere in this issue for the idealized case in which there 
are only two modes (rather than an infinite number of modes). 18 Even 
in this simple case, we do not know how good an approximation the 
perturbation solution provides. 

It is obvious that the various results of S. P. Morgan for metallic 
guide must be related to each other, even though they may appear 
somewhat dissimilar. In this section, after first developing the necessary 
theory for a long guide with many discrete mode converters, we show 
how the coupled line equations and Morgan's results for discrete mode 
converters may be derived from each other, and how perturbation 
theory derived from either the discrete case or from the coupled line 
equations yields equivalent results to perturbation theory applied 
directly to the field equations. 5 This discussion is intended to provide 
a better physical understanding of the coupled line equations themselves, 
as well as of the approximate solution that we use, than would be ob- 
tained by merely a formal treatment based entirely on the coupled 
line equations. We will often simplify the problem by including only a 
few of the spurious modes (sometimes only a single spurious mode), 
in addition to the TE i signal mode. While this procedure is useful in 
studying some of the basic problems, it of course does not provide a 
rigorous treatment for the real problem, which involves an infinite 
number of modes. However, experimental results for copper guides 
show that often only one or two spurious modes are present with sig- 
nificant magnitude, 8 and thus provide some additional justification 
for the study of the idealized problem. 

2.1 Scattering Matrices of Discrete Mode Converters 

While S. P. Morgan has performed a field analysis for discrete mode 
converters in lossless metallic guide, many of the general properties of 
the scattering matrices for discrete mode converters may be derived 
from conservation of energy, reciprocity, and the symmetry properties 
of the different mode converters. We give such a discussion in the present 
section, making use of Morgan's results 15,16 where necessary. Most of 
the discussion will be confined to cylindrical guides of infinite conduc- 
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Fig. 1 — Generalized mode converter. 

tivity. Wo choose the particular case of the discrete tilt to illustrate the 
general approach; briefer discussions are given for offsets and diameter 
changes. 



2.1.1 General Properties of Scattering Matrices 

Consider the general mode converter illustrated in Fig. 1. On the 
left of the obstacle we have the modes 0l,\l,2l, • • ■, and on the right 
the corresponding modes 0r,lr,2r,3r, •••, where the letters I and r 
stand for "left" and "right" respectively. The subscript zero will de- 
note the TEoi mode and the other subscripts will denote the spurious 
modes. This convention will be used throughout this paper. We assume 
that there are the same number of propagating modes on each side of 
the obstacle, thus ruling out cases in which one of the spurious modes 
is below cutoff on one side of the obstacle, above cutoff on the other 
side. Cutoff modes are neglected throughout this analysis. 

Denoting the normalized complex amplitudes of the modes incident 
on the obstacle from the right and left as I nr and /„/ respectively, and 
of those leaving the obstacle as R nT and R,a respectively, the matrix 
equations 23 of the obstacle may be written: 
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Morgan's results show that except very close to the cutoff of a spuri- 
ous mode, the power scattered from TEoi to the forward modes greatly 
exceeds the power scattered to the backward modes for "small" mode 
converters. Consequently, in the following treatment we shall neglect 
all reflected waves. Using this assumption, and the fact that S must he 
symmetric ( reciprocity) , 
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where S denotes the transpose of the submatrix S. Thus, 

Rl = SI r , 
Rr = Sh . 

If the obstacle is assumed lossless 

SS* = [1], 



(3) 

(4a) 
(4b) 

(5) 



where the * denotes the complex conjugate and [1] denotes the unit 
matrix. From (2) we thus have 



SS* = [1] or equivalent^ SS* = [1]. 



(6) 



2.1.2 Scattering Matrix for a Tilt 

Consider the tilt of Pig. 2(a). Neglecting reflected waves, this mode 
converter may be characterized by the matrices of (3) or (4). Since the 
tilt is symmetric about the plane A -A' of Fig. 2(a), its input and output 
terminals may be interchanged without altering its behavior. From this 
fact and (4), we have 

S = S; (7) 

the matrix »S' is symmetric for a tilt. If we consider a tilt in perfectly 
conducting guide, so that energy is conserved, we have from ((>) and (7) 



SS* = S*S = [1]. 



(8) 
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d »o 

(b) 

Fig. 2 — Waveguide with tilt. 

Next, consider a cascade connection of the tilt of Fig. 2(a) and an 
identical tilt rotated 180°, as shown in Fig. 2(b); we may consider that 
the first tilt has an angle +«, the second an angle —a. Let the scatter- 
ing matrix for the first tilt be S, for the second (rotated) tilt be S'. Then 
as the distance d between the two tilts approaches zero, the over-all 
scattering matrix of the two tilts becomes S'S (neglecting reflected 
waves). But it is clear from physical considerations that this cascade 
connection of two opposite tilts of equal magnitude must be equivalent 
to a straight piece of guide (of zero length); thus we must have 

S'S = SS' = [!]. (9) 

From (8) and (9) we have in the lossless case 

S* = S'. (10) 

By utilizing the symmetry properties of the various modes involved, 
further restrictions on the elements of »S' are easily found. For the sake 
of dcfiniteness, consider the case where the only modes considered are 
the signal mode, TE i , and the first -order forward spurious modes scat- 
tered by a discrete tilt from TE W , i.e., one polarization of the forward 
TMh and TK|„, modes, denoted by TM» + and TEi m + , as shown by 
Morgan. 1 "' 16 It is obvious by symmetry that only the (linear) polariza- 
tion having an asymmetric transverse field distribution with respect 
to the plane of the tilt will be excited by an incident TEw wave." While 
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the backward TEi„, modes, denoted by TEi m - , are excited to first order 
by a tilt, their magnitudes are much smaller than the magnitudes of 
the corresponding forward spurious modes, as stated earlier, and conse- 
quently we neglect the backward modes for the present. We recall the 
convention of (1), i.e. the top elements in the R and / column vectors 
of (4), having the subscript 0, denote the TEoi signal mode; the other 
elements denote the various spurious modes. Let the elements of the S 
and S' matrices be Sij and s,/ respectively. Thus, for example, if a unit 
TEoi wave is incident on the tilt of Fig. 2(a) from the left, a TE i wave 
of (complex) amplitude %> will emerge on the right. We now observe 
that rotation of the tilt by 180° leaves the TEoi mode unaffected, but 
reverses the sign of the field components of all of the spurious modes, 
since their field components vary as cos <p or sin <p. Consequently, the 
matrix components s,/ for the rotated tilt arc related to the matrix 
components s.-y as follows : 



i = 0, j ?± 0. 
**'' i 5*0, j = 0. 



sn' = 



+«.-/; 



i = o, j = o. 
i ?* o, j * o. 



(ii) 



From (10) we have 



««* = •*/■ (12) 

Equations (11) and (12) thus yield 



i = 0, j 9* 0. 
pure imaginary; { _^ Q . _ Q 



Sij = 

pure real; 



i = 0, j = 0. 
i * 0, j 9* 0. 



(13] 



The coupling coefficients between TE i and the spurious modes are pure 
imaginary, while all other matrix components are pure real. 

Summarizing the above results, the scattering matrix S of a discrete 
tilt in lossless metallic guide must satisfy the following relations, if we 
include only TE i and the propagating first-order forward spurious 
modes : 

S = S, s^ = sji. (14) 

SS* = [1]. (15) 
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i = 0, j ^ 0. 
pure imaginary; ' 

(16) 

. i = 0, j = 0. 

pure real; . ^ ^ Q 

As an example, let us determine the form of S for the case where only 
a single spurious mode is considered, in addition to the TEoi signal mode. 
(For a guide with a large intentional bend, the most significant spurious 
mode might be one polarization of forward TMu . For a guide with a 
small random straightness deviation confined to a single plane, the most 
significant spurious mode might be one polarization of forward TE12 , as 
discussed in Section 2.3 below.) Then (14) to ( 16) yield for the S-matrix 



S = 



(17) 



vi - c- jc 

L jc vi - c 2 J 

For small tilts in lossless metallic guide, Morgan 15,16 has given the 
coupling coefficients s , [or jc of (17)] in terms of the tilt angle a of 
Fig. 2(a) to first order in a, as follows; 

So; -jC t (j) + -a+ ••• , JVO, (18) 

where the first j on the right-hand side of (18) represents y/ — 1, the 
subscripts j' indicate the spurious mode. C, 0) + is a constant depending 
on the (forward) spurious mode; formulas and numerical values at a 
frequency of 55 kmc in 2-inch diameter guide for the Ct 's are given in 
Appendix A. In addition, the coupling coefficients C t ~ to the correspond- 
ing backward spurious modes are also given in Appendix A; as indicated 
above, these arc much smaller than the forward mode coupling coeffi- 
cients. 

Consider a tilt of angle ai followed by a tilt of angle «•» ; as the distance 
between the two tilts approaches zero, it is obvious that the structure 
approaches simply a single tilt of angle «i + a 2 • If S(a) is the matrix 
for a tilt of angle a, then the ^-matrix must satisfy the further require- 
ment that 

S(ai)-S(a->) = S( ai + a 2 ), (19) 

for every value of «i and as , where we again neglect reflected modes. 
Consider again the idealized case where only a single spurious mode is 
allowed, in addition to the TE i signal mode. We show in Section 2.3.3 
that the "exact" matrix for a large tilt for this idealized two-mode case 
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is given by 



8(a) = 



cos C t a j sin C t a 
j sin C t a cos C t a J 



(20) 



where C t is Morgan's coupling coefficient for the spurious mode in ques- 
tion. The matrix of (20) can readily be seen to satisfy the consistency 
condition of (19), and to approach the results of (17) and (18) for 
C t a « 1. However, the reader should be warned again that (20) will 
not be valid for large tilts in the physical case, because such large tilts 
will excite many spurious modes with significant magnitude, and hence 
can not be described in terms of only two modes. 

The above results for tilts in metallic waveguide include only one 
polarization of each of the forward spurious modes, i.e., TEi» and 
TMn + . We wish to extend these results to include both polarizations of 
each of these spurious modes so that we will be able to treat a long line 
containing arbitrary tilts with arbitrary angular orientation, i.e., not 
confined to a single plane. 

We first write the general results of (4b) as follows, dropping as un- 
necessary the subscripts / and r since we will always assume that all 
modes travel in the forward direction, from left to right. 



R] = IS]- II 



R] = 



BJ - 



pol 




~/ol 


R x ] 


, n = 


/J 


_AJ_ 




J V L 



Ru 
R* x 



R«] = 



Riu 
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, ZJ- 




, h] = 


hy 



(21) 



(22a) 



(22b) 



In these and following matrix relations, we adopt the convenient nota- 
tion that column vectors (n X 1 matrices) are denoted by the symbol ]. 
row vectors (1 X n matrices) by symbol ,_,, and square matrices by 
the symbol [ ] where it is not obvious from the context that something 
else is intended (e.g., the column vectors on the right-hand sides of (22a) 
and (22b)). The top elements in (22a), R u and /,, , represent the trans- 
mitted and incident TE„i wave. R x ], /,] and /?„], I„] are column vectors 
whose elements represent the two orthogonal (linear) polarizations of 
each of the spurious modes. 
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Fig. 3 — Geometry of a tilt in waveguide. 

Consider the tilt of Fig. 3, in which the two polarizations of each of the 
TEi„, + and TMn + spurious modes are defined along the x and y axes. 
To use the previous treatment including only a single polarization of each 
spurious mode, we must set 6 = in Fig. 3, so that the plane of the tilt 
lies in one of the planes defining the spurious mode polarization. Then, 
using the notation of (22), we may write the results for the ^-polarization 
as follows: 

Ro -%0 

_R,}_ ~ JC] 



r [D) 






(23) 



C] = 



JCj = Soy 



(24) 



The components of the column vector of (24) are given to first order by 
(18). In (23), we have partitioned the »S'-matrix to conform to the par- 
titioning of the R and / column vectors. The restrictions of (14) to 
(16) become respectively, in terms of the submatrices of (23): 

[D] = [D]. (25) 

+ !>/= 1, 



Sou 



:%>- [D])-c] = o], 



(26) 



C1-.0.+ [Df = [1]. 
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Soo 



C], [D] pure real. 



(27) 



In (26), [1] denotes the unit matrix. 

Now by symmetry an incident TE i wave in this particular tilt (6 = 0) 
will not excite the ^-polarization of any of the spurious modes as stated 
earlier, and by reciprocity the ^-polarization of any of the spurious 
modes will not excite TEoi . Further, the x- and ?/-polarizations of all of 
the spurious modes are uncoupled from each other by symmetry. There- 
fore, we may expand the matrix relation of (23) to include the ?/-polariza- 
tion of the spurious modes as follows : 



RA 



$00 ! i£ ! 0, 






jC) ! [D] [ [0] 



(28) 



where (18) and (24) to (27) still apply. 

Now let us rotate the axes by an angle — 6 with respect to the plane 
of the tilt, as shown in Fig. 3, and write the field in the guide in terms 
of modes referred to these new axes. The geometry of the rotated tilt 
may be specified precisely in the following way. Imagine that before 
the guide is tilted, lines <p = constant are drawn on the surface of the 
guide parallel to the guide axis. Then the tilt, of orientation 0, is con- 
structed at a specified point on the axis by tilting the guide in the plane 
defined by the axis and the <p = 6 line, by an angle a. If 



p = distance from the (tilted) guide axis in a plane perpen- 
dicular to the guide axis, 
s = distance measured along the (tilted) guide axis, 



(29) 



the three coordinates p, <p, s constitute "bent cylindrical coordinates," 
as used by Morgan for continuous bends confined to a single plane. 
In subsequent analysis for a guide with many tilts of arbitrary orienta- 
tion, we will adopt the convention that a > while is unrestricted; 
in contrast, for a guide with tilts confined to a single plane, we will set 
= and allow a to be unrestricted. 

Now let Ro] and l ] denote the fields with respect to the old axes and 
Re] and Ie] the fields with respect to the rotated axes. First, we note 
that TEoi is the same in both sets of coordinates. Next, we note that 
the field components of all of the spurious modes under consideration, 
i.e., TEi m + and TM U + , vary as cos <p or sin <p. If we call ipo the "old" and 
<pe the "rotated" angular coordinates, the transformation is simply 

W = <po + e. (30) 
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For each spurious mode this yields in terms of wave amplitudes 



R, 



lA.J 



cos sin 
—sin cos 



'R 



*t 



J \_Ryg J 



(31) 



with a similar result for the J's. Note that (31) applies to a single spuri- 
ous mode only (two polarizations); the index denoting the particular 
mode has been omitted for convenience. Thus, the R's (or 7's) in (31) 
are single (complex) numbers, and not column vectors. We may write 
the corresponding general rotation matrix including TEoi and all spurious 
modes in the following convenient form: 

Ro] = [M]-R e ], Jo] = \M}-Iel (32a) 



i - i ! i 

[M] = 0] | cos0-[l] ! sin0-[l] 

_0] ; -sin0-[l]|cos0-[l] J 
M = M~\ 



(32b) 



(32c) 



R and J in (32a) are as given in (22); the [l]'s in (32b) represent the 
unit matrix. We note from (32c) that ilJ is an orthogonal matrix. 
Now rewriting (28) as 



Ro] = [S ]-h], 
we substitute the relations of (32a) into (33) to obtain 
[M]-R e ] = [SoHM]-Ig] 

Ro] = liUrMSoHiVJ-Je] 
= [M]-[S ]-[M]-Ie}. 
Therefore, 



(33) 



(34) 





Re] = [S e ]-I 8 ], 




(35a) 


[So] = [M]-[S ]-[M]. 




(35b) 


Substituting [So] from (28) and [M] from (32b), we have 






[&] = 








Boo j cos • C j sin • C 




(36) 




j cos 6-C] ! cos 2 0-[D] + sin 2 8-[l] ! sin cos 0([D] - 


[1]) 




_ j sin 0-C] ! sin cos 6([D] - [1]) sin 2 0-[D] + cos 2 


Mil _ 
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as the scattering matrix for a tilt in lossless metallic guide, having an 
orientation of 6 radians with respect to the axes defining the polariza- 
tion of the spurious modes. C] remains as defined in (24), with com- 
ponents given to first order by (18), and the restrictions of (25) to (27) 
still apply. 

The above transformation may readily be seen to yield the correct 
results in a few simple cases: 

1. e = t. [S T ] is readily seen to be identical to the matrix [S'] de- 
scribed in connection with (11). 

2. 6 = 7r/2. Here we see that R x ] e=0 = R y \ =Tli , Ry\ B=0 = — #x] s=ir/2 , 
and similarly for the /'s, as is obvious from geometric considerations. 

The coupling coefficients between TE i and the x- and ^-polarizations 
of the/ 1 * spurious mode are from (36), (24) and (18) 

soj.x = jCt U) + -(a cos B), (37a) 

so,* = jC lU) + -(asm d). (37b) 

These general results simplify so that they may be simply expressed in 
terms of fixed x, y, z coordinates in the following special case : 

1 . The angular deviation of the guide axis from the z-axis is small. 

2. The unit vector perpendicular to the guide axis and lying in the 
plane defined by the guide axis and the <p = line is almost parallel to 
the .r-axis. These conditions insure that the <p = line drawn on the 
guide will remain almost parallel to the x-z plane. Subject to these 
conditions, we have approximately 

a cos 6 « a x , (38a) 

a sin 6 W a y , (38b) 

where a x and a v are the angles made by the projections of the guide 
axes, adjacent to the tilt, in the x-z and y-z planes respectively. 
These results are readily verified for the case shown in Fig. 3. They are 
derived in Appendix B. Substitution of (38) into (37) [and thence into 
the matrix of (36)] will greatly simplify certain later calculations. 

2.1.3 Scattering Matrix for an Offset and a Diameter Change 

The scattering matrix for an offset may be found by similar methods 
as used in Section 2.1.2 above for a tilt. Consider the offset of Fig. 4(a), 
with scattering matrix S; as above we assume the offset is small, so that 
backward modes may be neglected and only the propagating first-order 
forward spurious modes (i.e., TEi m + ) 15,16 need be considered. In Fig. 
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JL-. 



(a) 



r d *i 



(b) 

Fig. 4 — Waveguide with offset. 



4(b), this offset is cascaded with an identical offset, which may be 
derived from the original offset in two ways: 

1 . Reverse input and output terminals. 

2. Rotate the first, offset 180°. 

From 1 , the matrix for the second offset is S, the transpose of S. From 
2, the matrix for the second offset is »S", as defined in (11) above, by an 
argument similar to that used for the tilt. From these facts, as the dis- 
tance d separating the offsets approaches zero, 



SS = SS = [1], 
S'S = SS' = []]. 



(39a) 
(39b) 



Assuming lossless metallic guide, from ((>) 

SS* = SS* = [1]. (40) 

From (39) and (40) we have for an offset in lossless metallic guide: 

to -m. (4i) 



s,-j = pure real = 



i = 0, j * 0. 
+Sjl ' i * 0, j * 0. 



(42) 
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For a single spurious mode the matrix S for an offset becomes 



s = 



— c- 
■c 



Vl - c 2 _ 



(43) 



corresponding to (17) for a tilt. For small offsets in lossless metallic 
guide, Morgan 15,16 has given the coupling coefficients s y [or c of (43)] 
in terms of the offset b of Fig. 4(a) to first order in b, as follows: 



Soj — Cou) '" ~r 



7*0. 



(44) 



Cou) + is a constant depending on the spurious mode; formulas and nu- 
merical values for the C ±'s at a frequency of 55 kmc in 2-inch diameter 
guide are given in Appendix A. 

Analogous results to those of (19) to (38) for tilts are readily found for 
offsets, but will not be discussed in detail here. In particular, the scatter- 
ing matrix for an offset with an arbitrary angular orientation is found in 
the same way as given in (2 1 ) to ( 37 ) , making use of the rotation operator 
of (32). The geometry for offsets of arbitrary orientation is much simpler 
than for tilts. As before, imagine that lines <p = constant are drawn on 
the surface of the initially perfect guide, parallel to the guide axis. Then 
the offset of Fig. 5 is constructed at a specified point on the guide axis 
by translating the guide a distance b in the <p = plane without rotating 
the guide, so that corresponding <p- lines on the two sides of the offset 
are separated by a distance b. For a guide with many offsets of arbitrary 
orientation we take b > with unrestricted, while for a guide with 
offsets confined to a single plane we set = and allow 6 to be unre- 
stricted, as in the case of tilts. Then using the notation of Section 2.1.2, 
the scattering matrix S 9 for an offset of orientation 6 is given as follows: 



Se = 



Sao 


cos d-C^ 


sin 6-C^ 


-COS0-C] 


Qoa 2 d-[D] + sin 2 0-[l] 


sin 6 cos d([D] - [1]) 


-sin0-C] 


sin0cos0([.D] - [1]) 


sin 2 0-[D] + cos 2 0-[l] 



[45) 



where 



C] = 



Cj = Sqj 



(46) 
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bj, = b cos e 
by = b sin e 



Fig. 5 — Geometry of an offset in waveguide. 

The components of C] are given to first order by (44). Equation (45) 
(for offsets) has been written in such a way that the restrictions of 
(41) and (42) become identical to those of (25) to (27) (for tilts). 

The coupling coefficients between TE i and the x- and //-polarizations 
of the j th spurious mode are from (44) to (46) 



S 0j.y — C (j) • 



(b cos 0), 
(b sin 0). 



(47a) 
(47b) 



From Fig. 5 it is readily seen that b x = b cos and b u = b sin are ex- 
actly the x- and //-components of the offset, in analogy to the approxi- 
mate results of (38) for tilts. 

A similar treatment may be applied to a discrete diameter change. 
Here the spurious modes are the higher order TEo m ; 15 again for small 
offsets only the forward modes need be considered. This case differs in 
one fundamental respect from that of the offset. For an offset the TE i 
signal mode and the TEi„, spurious modes have a different angular de- 
pendence; for a diameter change the TE 01 signal mode and the higher 
order TE 0m spurious modes have the same (i.e., no) angular dependence. 
Thus, for a diameter change the signal and all spurious modes are coupled 
to each other to first order. 

First by conservation of energy (true in helix as well as copper guide, 
because the TEo„, modes are the same in both) (6) yields 



ss* = SS* = [1]. 



(48) 
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Next, consider the cascade connection of two identical diameter changes 
connected back-to-back. As in the case of an offset, this yields 

SS = SS = [1]. (49) 

From (48) and (49) we find that 

S = S* = pure real. ■ (50) 

However, there is no operation corresponding to the 180° rotation, used 
for tilts and offsets, because the signal and spurious modes have the 
same symmetry in the present case. While we have obtained no further 
information from general considerations than contained in (49) and 
(50), for the case of a single spurious mode the scattering matrix for a 
diameter change is identical to that of (43) for an offset. 

For small diameter changes, Morgan 10,16 gives the coupling coefficients 
Saj in terms of the change in radius Ar as follows: 

so; = C dlj) + -Ar+ ••• , 3 > 0. (51) 

da)* is given in Appendix A. 

There is of course only a single polarization of each spurious mode in 
this case, and consequently, there is no analysis in the present case cor- 
responding to those for tilts and offsets with two polarizations. 

2.1.4 Discrete Mode Converters in Helix Guide 

While the above results for diameter changes apply equally well to 
both helix and copper guide, those for tilts and offsets apply to only 
copper guide (strictly speaking, ideal lossless metallic guide). This is so 
because (6) no longer holds true in helix; energy is not conserved in 
helix, and in addition the various normal modes of helix are not even 
orthogonal with respect to power. 22 The coupling coefficients for discrete 
mode converters in helix have been obtained by linger. 18 However, one 
useful result is readily obtained from general considerations without 
performing a detailed field analysis. 

We show that at a discrete tilt or offset, the TE i transfer coefficient 
soo (the upper left-hand element in the scattering matrices for a tilt or 
offset given in Sections 2.1.2 and 2.1.3 above) is identical in both copper 
and helix guides, if we neglect backward modes and include all forward 
modes. This fact will permit the average TE i loss in helix guide with 
discrete random tilts and offsets to be readily calculated in terms of 
similar results for copper guide. 

The above statement is proven as follows. First consider an offset 
or tilt as a joint between a helix and a copper guide, with a unit TE„i 
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wave incident from the helix. If backward waves are neglected, the for- 
ward waves in the copper guide will be identical to those in a copper- 
copper joint because the incident TE i is the same in helix as it would 
be in copper. Therefore, the TE i transfer coefficient s o is the same for 
this helix-copper joint as for a copper-copper joint. Now reverse input 
and output terminals, so that the TE M is incident from the copper guide, 
and TEui and spurious modes travel away from the joint in helix. The 
spurious modes are now quite different than before, since they must be 
normal modes of helix guide; however, by reciprocity s 00 must remain 
the same. Finally, we may replace the copper guide containing the inci- 
dent TEoi with helix without further altering the fields in any way. 
Thus if backward modes are neglected, s q is identical in helix and copper 
guides with identical tilts or offsets. 

This conclusion has been verified experimentally by the authors in 
2-inch diameter helix guide at a frequency of 55 kmc. It has also been 
verified by a field analysis by H. G. Unger. 20 

2.2 The Discrete Case — Single. Spurious Mode 

We next apply the results of Section 2.1 to the study of long guides 
with many discrete mode converters separated by guide sections that 
are ideal, i.e., geometrically perfect right-circular cylinders. We restrict 
our attention to the case of a single spurious mode, in addition to the 
TEui signal mode. (If the spurious mode is polarized, such as TEi m + , 
we consider only one of its linear polarizations.) 

The over-all transmission matrix for such a guide with N discrete 
mode converters consists of a product of 2JV matrices, one matrix for 
each mode converter (as given in the preceding section) and one diagonal 
matrix for each section of ideal guide. These matrix results are then used 
to derive a perturbation theory, valid when the mode converters are 
sufficiently small, that greatly simplifies further calculations. 

2.2.1 Matrix Analysis 

Consider the guide illustrated schematically in Fig. 6. This guide con- 
sists of N discrete mode converters separated by N sections of ideal 
copper guide, of length h- . Since we consider only a single spurious mode, 
if the discrete mode converters are tilts or offsets they must lie in a single 
plane (taken to be the <p = plane in the notation of Section 2.1). We 
seek the response of such a guide to a unit input TE 0] wave. 

The scattering matrix Sk for the A" 1 mode converter is given in Section 
2.1 for the case of a single spurious mode, i.e., (17) and (18) for tilts, 
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Io(o)- 
Ii(o)- 



r — L2 ~"i 

|*"Li->| 



s, 



3 



s k 



Sk+i Sn 



Io(N) 
I, (N) 



M, 



M, 



M k 



M„ 



m=i 

Sk = SCATTERING MATRIX FOR k TH MODE CONVERTER 

Mk = SCATTERING MATRIX FOR k TH MODE CONVERTER 
PLUS FOLLOWING SECTION OF GUIDE 

I = SIGNAL MODE (TE 01 ) NORMALIZED AMPLITUDE 
I, = SPURIOUS MODE NORMALIZED AMPLITUDE 

Fig. 6 — Waveguide line with discrete mode converters separated by ideal 
guide sections. 

(43) and (44) for offsets, (43) and (51) for diameter changes. In the 
ideal guide sections connecting the discrete mode converters, the signal 
and spurious modes propagate independently with their respective 
propagation constants; the scattering matrix W k corresponding to the 
k th section of ideal guide, of length h as shown in Fig. 6, is given by 

where r is the propagation constant for the TEoi signal mode, I\ the 
propagation constant for the spurious mode. Thus, the scattering matrix 
M k for the k th mode converter plus the (following) fc th section of guide, 
as shown in Fig. 6, is 

M k = W k -S k . 



W k = 



e 



For tilts M k is given as 



r e -*v* 



M k = 



Vl - ci? 



je- r ° lk -c k 



,-iVfc 



e ' w -\/l - c fc 2 _ 



with c k given by (18). For offsets and diameter changes 

c k - 



M k = 



ryj'.fc.y'fTr 



e ° -c k 



-e- r ^-c k 



™- r i'* 



Vl - c fc 2 _ 



(53) 



(54) 



(55) 
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with c k given by (44) and (51 ) respectively. 

The over-all transmission matrix T for the entire guide, N sections 
long, will he 



T = \l m l n ] = II ,1/, = M.-M^ ■ ■ ■ M,-M x , (5 



(•») 



and the output mode amplitudes will be given in terms of the input mode 
amplitudes by 



UN) 
Ii(N) 



= T 



7 (0)" 

/i(0) 



(57) 



We assume the guide is excited by a pure TEoi wave of unit amplitude 
and zero phase; the initial conditions on (57) become 



/o(0) = 1, 7,(0) = 0. 

Then, 

h(N) = r 00 /o(0) = Ton, TEoigain. 

Ii(N) = 7'i(,/u(0) = 7\o , TEiu-spnrious mode transfer 

coefficient. 



(58) 

(59) 
(60) 



It will subsequently be convenient to normalize these quantities as 
follows: 



_-r u /-» 



T m = o ° 4ff -G ; Go = TEw normalized gain. (61) 

Tiu = e~ ,r ' N -Gi ; G x = TE i-spurious mode normalized transfer 



coefficient. 



(62) 



In each case, the propagation factor of the corresponding mode has been 
removed; as shown in Fig. 6, L N is the total length of the A^ sections of 
guide being considered. Since TE i has a lower heat loss (i.e., a = Re r ) 
than any other mode, in a physical guide Go ^ 1. 

The exact solution above is of limited value both because of its com- 
plexity and also because, as discussed in Section 2.1, the available ex- 
pressions for the coupling coefficients of discrete mode converters are 
valid only to first order. Consequently, we seek approximate expressions 
for Too and Tm , or equivalently G and Gj , valid when the coupling 
coefficients are sufficiently small that the guide departs only slightly 
from ideal. Under these conditions, it will be shown that G\ is of first 
order and G„ departs from unity only to second order. 
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It is convenient to write 

G = 1 - P , (03) 

where p, the (complex) departure of G from unity, will be of second 
order. We further define the TE i normalized magnitude g, complex 
loss A, loss (in nepers) .4, and phase as follows: 



G = ge je = c" A = c-V e = 1 - p; 
A = A - jQ. 



(64) 



Then if the coupling coefficients are sufficiently small so that | p | « 1, 
we have to second order: 

gttl - Re p. (05a) 

A w p. (65b) 

A « Re p. (65c) 

« - Im p. (65d) 

2.2.2 Perturbation Theory 

Consider the transmission matrix T of (56), with the M k given by 
(54) or (55). Let us expand the square root in the diagonal elements of 
(54) or (55) in a power series as follows; 

vr^" 2 = i - W - W + • • • . (66) 

It is apparent from the rules of matrix multiplication that the com- 
ponents of T may be expressed as power series in the c fc 's. Since each c k 
may be expressed as a power series in the appropriate geometric parame- 
ter, with the first term given by (18), (44) or (51), we can thus obtain 
expressions for Gi and G'o as power series in the geometric parameters. 
The first-order results of (18), (44) and (51) are sufficient to give #i 
to first order, G to second order; if the mode converters are sufficiently 
small we may hope that these results will give a valid approximation 
for the TEoi gain and the TE„i-spurious mode transfer coefficient. 

We first determine the first few terms of expansions for T u , and Too . 
For convenience, we write for Mk in (56) 

~A k e," 

>. B*. 

where A k , t k , h and B k are determined by comparison with (54) or 



M k = 



(67) 
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(")")). For small coupling, A k and B k will have magnitudes a little less 
than 1, while e k and 8 k will have magnitudes much smaller than 1. In 
writing out 7% and Ton , we group terms according to the number of 
"small" quantities (i.e., t's and 6's) they contain. Thus, we have: 

.V 

r M = Ea, ••• .i,-,s,b, +1 ••• Bs 

i=l 

+ E E E -1. ■••-l I -i5 1 B 1+1 ..- (68) 

i'=l ; = i+i *=/+l 

■ • • Bj-iejAj+i • • • A k -i8kBk+i ■ ■ • B N + ■ ■ ■ . 

N—l N 

7'oo = Ai • ■ • A N + E E Ai • • • -1 1 _ 1 5,S 1+1 • • • B>_ie^ >+1 • • • .4* 



s— :i JV — 2 .V— 1 tf 

+ E E E E iii • • ■ -4/- 1 S,£, + i • • • Bj^ejAj^ • ■ • (69) 
,-=1 ,=,+1 k=i+l l=k+l 

■ • • Ah-ifaBk+i ■ • • Bt-ieiAi+i ■ • • A N + ■ • • • 

Consider first (68) for T w . Referring to (67), (66), (54) and (55), 
the successive terms of (68) (of which we have written down the first 
two) contain components of the following orders in the c k : 
term order of components 

1 C, C , C , ■ • • 

,, 3 5 7 

1 c ,c ,c , • • • 



Equation (18), (44) or (51 ) gives the c* to first order in terms of the 
tilts a k , offsets b k , or diameter changes Ar* . Therefore, we may obtain 
7'ki to first order in the geometric parameters by retaining only the com- 
ponents of order c, which occur in only the first term of (68) (the single 
summation), and using Morgan's first -order coupling coefficients. In- 
deed, since we do not have the coupling coefficients to more than first 
order, it would be totally unjustified to retain any additional components 
in (68 ), and in particular, to seek an exact result for the matrix multipli- 
cation of (56) via numerical techniques or otherwise, without first ob- 
taining the coupling coefficients c k to higher order. 

Similarly, the successive terms of (69) for T u0 contain components of 
the following orders in the c k : 

term order of components 

- 2 4 

1 c , c , c , • • • 

,-. 2 4 6 

2 c ,c ,c , ■•■ 



1056 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1962 

We obtain Too to second order in the geometric parameters by retaining 
only the components of order c and c in the first term and the com- 
ponents of order c" in the second term, and using Morgan's first -order 
coupling coefficients. As above, any additional terms would be unjusti- 
fied. We note that there are no first -order terms present in T m . 

We thus obtain the following approximate results for G and G\ , the 
normalized transmission parameters defined in (61) to (64), from (66) 
to (69), (54), (55) and (18), (44) and (51). In these results we use the 
differential propagation constant AI\ defined as 

Ar = Aa + j'AjS = r - r t . (70a; 

Aa = a - ai , A/3 = 0o - ft . (70b) 

The real and imaginary parts of Ar, Aa and A/3 respectively, are called 
the differential attenuation constant and the differential phase constant 
respectively. Since the TEoi signal mode has lower heat loss than any 
other mode, Aa < throughout the present paper. The geometry of the 
guide is shown in Fig. 6; L k -i is the length of guide up to the k th mode 
converter, l k the distance between the k th and the (k + l) th mode con- 
verters. Xi denotes Morgan's first-order approximation for the coupling 
coefficient of the i th mode converter. 

For tilts and offsets the diameter of the different guide sections is of 
course identical, so that AT is strictly constant. This is not true for di- 
ameter changes, so that strictly speaking we should include in the analy- 
sis the fact that AT changes from section to section. However, we assume 
that the total range of the guide diameter is very small, centered about 
its nominal value, and neglect the small changes in Ar in all following 
analysis, both for the discrete case and for the continuous case. This ap- 
proximation is not necessary; it would be possible to include the varia- 
tion of AT in the present analysis without great difficulty. However, we 
choose to ignore this effect without careful study in the interests of 
simplicity. 

i. G — Tilts, Offsets and Diameter Changes 
Go = 1 - p. (71a) 

N 
P W 2 



X> 2 +£ £ w^'-*-*-' 1 



1 ' (71b) 

N N 






._,. Afl Ly_ ,-/,,_, I 



= *E*/+ E^ rt E^- 



KNOWN IMPERFECTIONS 1057 

Special case — equally spaced mode converters 
h = k , L, = kk . (71c) 

,«iE..v+ E E .f,.r/ r '°"-" 

1=1 i=l /=.■+! 

= ieew^^*' 1 <7id) 

.=1 ;=1 

N ff—1 

E ** + E 

i=l A=l 

C-fli, tilts. 
Xi = C„-bi, offsets. (71e) 

Cd-Ar,-, diameter changes. 
m. Gi — rtto 

(?i«iE^- 4ri '-'. (72a) 

i=i 

Special case — equally spaced mode converters 

h = 7 , L k = kl ■ (72b) 

N 

-Ar/n(t-l) 



ft«i2>* • (72c) 

i=i 

Xi = Crot { . (72d) 

m. G t — Offsets and Diameter Changes 

Cn^-Z x*-" 1 "- 1 . (73a) 



Special case — equally spaced mode converters 

/* = /„, L k = kU. (73b) 

d«- i>, c - Ar, ° ( *- i) . (73c) 

i=i 

C -6,-, offsets, 
.r,- = (73d) 

Cd • Ar,- , diameter changes. 
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2.2.3 Discussion 

The above results for tilts and offsets apply strictly to only lossless 
metallic guide, for which A« = 0. However, we expect on intuitive 
grounds that they provide a satisfactory approximation for real copper 
guide. The results for diameter changes apply equally well to both cop- 
per and helix; for tilts and offsets in helix we must use the helix guide 
coupling coefficients, but otherwise the analysis is the same. 

The results of (71 ) to (73) yield the TE i signal mode to second order 
(there are no first-order terms), the spurious mode to first order, in the 
appropriate geometric parameters. In each case we have determined the 
first correction term to the solution for a geometrically perfect guide 
with all Ck = 0, i.e., G (z) = 1, Gi(z) = where z is distance measured 
along the guide axis. 

A rough physical interpretation may be given for the approximate 
results of (71) to (73). The spurious mode at any point in the guide is 
regarded as a sum of waves arising at each mode converter. Each of 
these waves is computed by assuming an incident TEoi wave at each 
mode converter identical to the TE i wave that would be present in a 
perfect guide, with all Ck = 0, and further assuming that the converted 
spurious mode wave is unaffected by subsequent mode converters, i.e., 
propagates as it would in perfect guide. The TE i signal mode is re- 
garded as the sum of three components: 

1. The TEoi wave that would be present in a perfect guide. 

2. The signal lost from TE i at each mode converter, assuming an in- 
cident TEoi wave identical to that in perfect guide and no incident 
spurious mode wave. 

3. A sum of waves reconverted to TEoi from the spurious mode at 
each mode converter. Each of these waves is computed by taking the 
approximate spurious mode as computed above, and assuming that the 
reconverted TE i wave is unaffected by subsequent mode converters. 
It is easy to see that this component may be expressed as a sum over 
all pairs of mode converters. 

The approximate solutions of (71 ) to (73 ) may be regarded as the initial 
terms of power series expansions in the Xi's (or equivalently the a,'s, 
bi% or Art's). It is reasonable to assume that if the xfa are sufficiently 
small these power series will converge sufficiently rapidly so that their 
first terms will provide valid approximations; under these conditions we 
would expect that | p | « 1 , | Gi \ « 1 . However, to use these results with 
confidence we must have bounds on the errors introduced by these ap- 
proximations that will give quantitative information on the waj r in 
which the solutions of (71) to (73) approximate the true solution. We 
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postpone consideration of these questions to Section 2.3 below, where 
we study the continuous case, since we will show there that the discrete 
and continuous cases are closely related to each other. It turns out that 
although we can give bounds on the higher terms of the series expansions 
that are in a sense the best possible, we still lack sufficiently precise in- 
formation to make any useful statement on the way in which the pertur- 
bation solutions approximate the true solutions. Consequently, confi- 
dence in the accuracy of practical calculations using the results of (71 ) to 
(73) can at present be justified only on intuitive grounds. 

In the region where we hope the perturbation solution is valid, 
\p\ « 1. Under these conditions, the TE i normalized magnitude g, 
loss A, and phase 8 are given approximately by (05), where p is given 
by (71). If | p | is not much smaller than 1, it would seem that the ap- 
proximation of (65) must be invalid. However, it will be shown in Sec- 
tion 2.3 on the continuous case that (05b), (05c), and (05d) together 
with (71 ) remain plausible for a wide class of interesting cases where 
the magnitude of the right-hand side of (71b) or (71d) becomes much 
greater than 1 and (05a) fails, although here again the justification is 
no more rigorous than that for (71 ) to (73) when | p | « 1, | G\ | « 1. 
This extension is important, for otherwise we should be limited to con- 
sidering only cases where the total loss | A |, and hence, A and | 8 |, 
are small. 

2.3 The Continuous Case — Single Spurious Mode 

In this section, we study mode conversion caused by distributed geo- 
metric imperfections, such as continuous curvature of the guide axis or 
continuous variation of the guide diameter. We discuss briefly the gen- 
eral telegraphist's equations or coupled line equations for the general 
case, but again restrict the detailed treatment to the case of a single 
spurious mode. (If the spurious mode is not TE „, , we consider only one 
of its linear polarizations.) 

It is obvious that a close relationship must exist between corresponding 
discrete and continuous cases. A guide with an arbitrary curved axis 
may be regarded as the limit of a guide with many small discrete tilts; 
similarly, a guide with a continuous diameter variation may be regarded 
as the limit of a guide with many small discrete diameter changes. We 
will show how the matrix equations for the discrete case and the differ- 
ential equations (coupled line equations) for the continuous case may 
be obtained from each other by suitable limiting processes. 

We next discuss the perturbation theory for the coupled line equations, 
and show its relationship to the discrete perturbation theory of Section 
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2.2. An improved perturbation theory for the TEoi complex loss 
A = A — j& [see (04)] is given that permits the treatment of many 
eases of practical interest. 

For moderate lengths of copper guide the differential loss is small com- 
pared to 1, | Aa | L « 1, where L is the guide length. We obtain an ap- 
proximate treatment for this case by setting Aa = 0, as in lossless me- 
tallic guide. The TEoi loss A for this case may be expressed in terms of 
the Fourier coefficients of the coupling coefficient c(z); this treatment 
makes evident the relationship between the power spectrum of the 
coupling coefficient and the corresponding transmission statistics of the 

TEoi loss. 

Finally, we discuss Morgan's coupling coefficients for general continu- 
ous cross-sectional deformations of copper guide, and the relationship 
between the various copper guide coupling coefficients. 

2.3.1 Generalized Telegraphist's Equations 

By means of the generalized telegraphist's equations, Maxwell's 
equations for the fields in a deformed guide may be expressed in terms 
of the normal modes of the undeformed guide; 

^= Ek /„. (74) 

dZ n 

I m represents the normalized complex amplitude of the m th mode and 
the summation is extended over all modes. K mm represents the propaga- 
tion constant of the m th mode; k„„, , with m ^ n, represents the coupling 
coefficient between the n th and m th modes. The n mn are functions of the 
geometry of the guide imperfection. Of course, (74) is in rather general 
symbolic form; in particular problems a more specific notation is often 
used to denote the various modes. For example, each subscript in (74) 
often becomes a double subscript, to conform to the usual ways of in- 
dexing waveguide modes. Square brackets enclosing a subscript pair 
are often used to indicate a TE mode, round brackets a TM mode; some 
sort of notation, such as ± superscripts or different symbols, is used to 
differentiate between forward and backward modes. 2 

The K mn have been evaluated for a variety of different types of imper- 
fections in a variety of guides, as discussed in the introduction to Sec- 
tion II. Morgan has given the K m „ for curved metallic waveguide filled 
with an inhomogeneous dielectric. 2 For a homogeneous lossless dielectric, 
the propagation constants K mm are simply equal to the corresponding 
propagation constants in undeformed metallic guide; 

Kmm = ±i/3„. (75) 
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Morgan finds that curvature couples the TE i + to only TM U + and 
TEinA (This should not be taken to mean that no other spurious modes 
will be excited in a curved guide with a homogeneous dielectric; such 
other modes will indeed arise, e.g., because of the coupling between TEi„, 
and TE 2m . However, in good guides these other modes will normally be 
negligible compared to the first -order spurious modes.) Using his results 
the k's of (74) involving the TE i + mode may be expressed in terms of 
the CVs of (18) and Appendix A as follows, for a curved circular guide 
with a homogeneous dielectric: 

iC ± 

K[01] + llm]-t = K[lm]±|01]+ = ' "' , (76a) 

iC + 
K[oi]+<in+ = K(u)+[oi]+ = ' D '/ (11 \ • (76b) 

R{z) 

R(z) in (76a) and (76b) is the radius of curvature of the axis of the 
curved guide, with z measured along the curved axis; the curvature of 
the axis is confined to a single plane. Similarly, linger 10 has evaluated 
the K m „ for a tapered metallic waveguide as follows :f 

1 fcoJta /3„„ ± fa*. dr(z) 

rh ,r ~ from 2 V#>,„/3o„ dz (77) 

m ^ n for the upper signs. 

In (77), r(z) is the radius of the guide, the fc's are the Bcssel roots of 
the modes in question (see Appendix A). When the incident mode is the 
TE ()1 , (77) becomes 

K[oi]+[omi± = Cdim* • -^n~> m ^ !> (78) 

ill terms of d of (51 ) and Appendix A. Finally, we note from (77) that 
if both modes are forward modes, 

K[0ml + [0; 1 ]+ = — K[0n] + [0m]+ , (79) 

which, as we shall subsequently see, must be true by (43). 

Various geometric imperfections in helix guide have been similarly 
treated by Unger. 3192021 

The generalized telegraphist's equations of (74) contain an infinite 
number of modes. To simplify the problem, we wish to approximate 
the true situation with a finite number of modes; indeed, in much of the 
present study we include only two modes, the TEoi signal mode and the 

fThe sign of the second term inside the summation in (40) of Reference 10 
appears to be in error. 
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most important spurious mode. It is obvious that this approximation 
must distort the results in some respects; however, we hope that if the 
significant spurious modes are included, any errors will be small. Thus, 
in (74) we will set all /'s equal to zero except those corresponding to the 
TEoi + signal mode and one or a few of the forward spurious modes. 
Other ways of approximating the true situation with a finite number of 
modes might of course be possible, but we choose this one as the simplest. 
In (70) and (78) above we observe that the coupling coefficients for 
continuous bends and diameter changes may be expressed in terms of 
the coupling coefficients for discrete tilts and diameter changes. This, 
plus the obvious fact that a continuous deformation may be regarded 
as the limit of many small, closely spaced discrete imperfections, sug- 
gests a closer study of the relation between the discrete and continuous 
cases. In Sections 2.3.2 and 2.3.3 below we study the correspondence be- 
tween the coupled line equations for the continuous case and the matrix 
equations of Section 2.2.1 above for the discrete case. 

2.3.2 Transformation from the Discrete to the Continuous Case 

Consider a lossless metallic guide containing many small discrete mode 
converters spaced an equal distance Az apart. Fig. 7 illustrates three 
cases of interest — tilts, offsets, and diameter changes — with of course 
a greatly exaggerated vertical scale. c k is in each case the conversion 
coefficient in the corresponding matrix, (17) or (43) for tilts or offsets 
and diameter changes respectively; c k is given by (18), (44), and (51) 
for tilts, offsets, and diameter changes respectively. Then it is clear that 
by the proper limiting process, in which Az — > and the discontinuities 
become smaller and more closely spaced, we may approach the continu- 
ous deformations illustrated in Fig. 7 by the dotted lines. 

The limiting continuous deformation of Fig. 7(a) corresponds to a 
guide with a continuously varying curvature of its axis. The continuous 
deformation of Fig. 7(b) corresponds to a waveguide made of very thin 
circular punchings which may slide with respect to each other like a 
stack of cards. Finally, the continuous deformation of Fig. 7(c) corre- 
sponds to a continuous taper. The case of Fig. 7(a) will be important in 
determining the effect of random straight ness deviation on the TEoi 
transmission. For the continuous bend and offset of Figs. 7(a) and 7(b) 
the guide axis must lie in a single plane if the analysis is to be restricted 
to a single polarization of the most important spurious mode (TE, 2 ). 

For purposes of illustration we consider the transition to the continu- 
ous bend, illustrated in Fig. 7(a). A larger-scale drawing for this case 
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Fig. 7 — Waveguides with tilts, offsets, and diameter changes. 

is shown in Fig. 8; only the center line of the continuously bent guide 
and of* the tilt approximation to it are shown, for the two adjacent sec- 
tions lying on either side of the /."' tilt. 

Referring to Fig. 6, we set the distance between discrete tilts l k equal 
to Az. Let h(z) and I x (z) denote the TE (n and spurious mode complex 
amplitudes just to the left of the k th tilt; then I {) (z + Az) and h(z + Az) 
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Fig. 8 — Fig. 7(a) in larger scale. 

denote the corresponding complex mode amplitudes just to the left of 
the (k + l) th tilt. Then taking (17) for the scattering matrix of the tilt, 
and setting l k = Az in (52), (53), and (54), we have 



7o(z + Az)' 
h(z + A*) 



r -r n Az /: 5 • -r„Az 

e ° VI - c fc 2 je c k 



-r, Az 



'IM' 



r VlA ° Vi - c k \ JM_ 



je c k e 

where c k is the conversion coefficient for the k th tilt. From ( 18' 

Ck = CfCik + • • • , 



(80) 



(81 



where C t is Morgan's coupling coefficient for tilts for the spurious mode 
in question, given in Appendix A. From Fig. 8 



Thus, 



AZ 4- 
ti k 



Cfc = r— AZ + 



(82) 



(83) 
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where the dots in (82) and (83 ) represent terms of higher than first order 
in Az. 

Expanding the exponentials and the square roots in the matrix of (80) 
in power series [see (66)], making use of (83), and writing the right- 
hand sides of the resulting equations as power series in Az, we have 

/„(* + Az) = h(z) + |-r„/ (z) +j£/i(«)JAb+ ••• 

(84) 
Zi(* + A») = /,(«) + -T l I l (g)+j^lM A2+ .... 



We transfer the terms I a (z) and h(z) to the left-hand sides of these 
equations and divide by Az. Then as Az — ► 0, R k —* R(z), where R(z) 
is the radius of curvature of the guide axis, and (84) becomes 

I '(z) = -YMz) +jc{z)h(z), 

(85) 
/,'(*) = jc(z)h(z) - Whiz), 

q 

c(z) = j-j— coupling coefficient , (86) 

H\z) 

where the primes denote differentiation with respect to z, and R(z) is 
the radius of curvature of the guide axis. 

We now compare (85) and (86) with (74) to (76); in the case of lossless 
metallic guide, where T = j(3 , T t = jfa , the two sets of equations arc 
identical if we retain only two modes in the results of (74) to (76). Thus 
by taking the proper limiting form of Morgan's results for discrete tilts 
we arrive at Morgan's results for continuous bends obtained via the 
generalized telegraphist's equations, in the two-mode case. It is easy 
to see that this method extends readily to additional first-order spurious 
modes. It should in principle be possible to include as many modes as 
desired in this type of discussion, but such calculations have not been 
actually carried out. We will be content in the present paper to take 
the two-mode model as suggestive of these more general results. 

As in the discrete case, we assume the guide is excited by a pure TE i 
wave of unit amplitude. Thus, the initial conditions on the differential 
equations (85) become 

/o(0) = 1, A(0) = 0, (87) 

corresponding to (58) for the discrete case. Then h(z) will be the TEoi 
gain, h(z) the TE w -spurious mode transfer coefficient, corresponding to 
(59) and (60). 
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It proves convenient to follow the normalization used in the discrete 
case [see (61) and (62)], removing the propagation factor of each mode 
in ideal guide; 

/,(*) = e- r ° ! -G (z) (88) 

h{z) = e- v ^-(h{z). (89) 

#0(2) and G\(z) are again the TE01 normalized gain and the TE or 
spurious mode normalized transfer coefficient. Substituting the trans- 
formation of (88) and (89) into (85), we have: 

Go'(z) = jc(z)e AT 'G 1 (z) 

C ri '(z) = jc(z)e- Arz G l) (z) (90) 

AT = r„ — r x , differential propagation constant. 

C 

c(z) = jjt-t, coupling coefficient. (91) 

l\\z) 

The initial conditions of (87) become 

Go(O) = 1, (?i(0) = 0. (92) 

A similar treatment may be given for the limiting case of many dis- 
crete offsets or diameter changes, illustrated in Figs. 7(b) and 7(c) re- 
spectively. For continuous offsets we find : 

/„'(*) = -r /o(2) +c(*)/iOb), 

h'(z) = -c(z)h(z) - tjm. 

c(z) = Co-x'(z), coupling coefficient. (94) 

C is Morgan's offset coupling coefficient for the spurious mode in ques- 
tion, given in Appendix A. x(z) is the transverse displacement of the 
guide axis, as illustrated in Fig. 7(b); primes of course denote differen- 
tiation with respect to z. Equations (93) and (94) must agree with the 
generalized telegraphist's equations for a guide with a continuous offset, 
when restricted to two modes. Using the transformation of (88) and 
(89), (93) and (94) become: 

Go'(z) = c(z)e AV V l (z), 

(95) 

<V(«) = -c(z)c-* Tz G (z), 

AT = To — Ti ; differential propagation constant. 

(96) 
c(z) = Ca-x'(z), coupling coefficient. 
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The continuous diameter change may be similarly treated by replacing 
x'(z) by r'(z), C by C d in (93) to (96), where r(z) is the guide radius. 

2.3.3 Transformation from the Continuous to the Discrete Case 

We next determine the scattering matrix for a discrete mode converter 
by considering the limiting form of the generalized telegraphist's equa- 
tions as the continuous coupling coefficient approaches a 5-f unction. For 
purposes of illustration we consider the case of a discrete tilt, regarded 
as the limiting form of a continuous bend. We again restrict our treat- 
ment to only two modes, the TE i signal mode and a single spurious 
mode, for simplicity. 

Fig. 9 shows a tilt of angle a obtained as the limit of a continuous bend. 
From (8(5) or (91) the continuous coupling coefficient is given by 



c(z) = — , 



, 2 < ; 

< z < A; 



0, A < z. 
lim c(z) = C t a-8(z 



(97a) 



(97b) 



Thus the coupling coefficient is zero for z < and z > A, constant in the 
region < z < A, z being measured along the guide axis as usual. Now 
exact solutions for (85) or (90) are known for only a few special c(z); 
indeed this is why we are forced to use approximate solutions of various 
type's in the present paper. One important case in which exact solutions 
are available is the present case e(z) = constant; 1 this is obvious from 
the fact that the equations (85) then become simply two simultaneous 
differential equations with constant coefficients, whose solutions are well 
known. 

The solution for (85) with e(z) = c is given as follows: 



h(z] 



= T{z) 



U0) 

A(0) 



T(z) = 



/ r„ + r, \ 

expl - — ^ — z ) 



K+ - K 
-K-e« r,2)z '+K 



(-ir/2)-v 



B UTm ^-B 



(-Ar/2).v 



ur/2)r,/ 



(AI72W 



(-sr/2) Zy / 



(98a) 



(98b) 



- K-e 



(-Ar/s)*v 
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LIMITING DISCRETE TILT 

Fig. 9 — Tilt of angle a as limit of continuous bend. 



K ± = -j 



A±V 



K+K- = - 1 



Co 

Ar 



•■+ 



K + - K- 

V = A 1 



. co/Ar 



m- 



:98c) 



(98d) 



(98e: 



These results are applied to the situation of Fig. 9 and (97) by setting 

z-» A, 



(99) 



Co 



Thus 



r/o(An 

U(A)J 



= T(A) 



7o(0)" 
/i(0). 



(100) 



Now in (97) and Fig. 9 let us keep a fixed and let A — > 0, R — > 0; then 
the curved guide of Fig. 9 will approach a discrete tilt of angle a, and 
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the matrix T should approach the scattering matrix of this discrete tilt, 
From (98) and (99) we have: 



limKfc = ±1 

A-0 

lim F W = 9 

A-,0 A + — /V_ Z 

,. / r + r, A 

hm exp I — — A 1 = 1 

a-o \ 2 / 

Then from (98) to (101) we have 

|~ COS Ci a j sin C t a 
lim T(A) = T(0) = 
a-o \_j sin Cta cos C t a _ 



:ioi) 



102a) 



for 



\\mc(z) = Cta-S(z) (102b) 



A-0 



as the scattering matrix of a discrete tilt of angle a. This is the result 
given in (20) as the "exact" two-mode scattering matrix for a finite tilt. 
We have thus found the scattering matrix for a two-mode, finite tilt 
starting from the scattering matrix for an infinitesimal tilt by first pass- 
ing to the limit of continuous mode conversion and then transforming 
back to the discrete case by allowing the continuous coupling coefficient 
to approach a 5-function. Alternately, we have found the scattering 
matrix that satisfies the following requirements (for a single forward 
spurious mode ) : 

1 . Conservation of energy and symmetry. 

2. Agrees with Morgan's small tilt theory [see (17) and (18)]. 

3. Satisfies an additional requirement for finite tilts, so that the 
matrix for the sum of two tilts equals the product of the matrices of the 
individual tilts [see (19)]. In practice the above matrix for a finite tilt 
will of course not provide an exact description for large tilts because 
of the presence of additional spurious modes. 

An offset and a diameter change may be similarly treated. For a finite 
offset of magnitude b we have corresponding to (102) as the "exact" 
two-mode scattering matrix 
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cos C b sin C b 



lim T(A) = T(0) = 



.i-0 



•sin o 6 cos C„b 



(103) 



Diameter changes may be similarly treated if they are small enough so 
that the variation of AT may be neglected. 

The present analysis readily extends to include additional first-order 
spurious modes. It should in principle be possible to generalize it further 
to include as many modes as desired, but such calculations have not ac- 
tually been carried out. 

The astute reader may have noted a potential difficulty with the 
analysis above. This potential paradox may be stated as follows. The 
above discussion shows that the continuous bend and the discrete tilt 
of Fig. 9 are approximately equivalent to each other so long as the fol- 
lowing conditions are satisfied: 

2Cta »1. (104a) 



Ar l -a 

Ar 



A«l. (104b) 

z 

But if (104a) is not satisfied it is no longer obvious that the discrete 
tilt and continuous bend remain approximately equivalent. Stated phys- 
ically, consider a short continuous bend, with length A fixed and small 

enough so that - — — - • A « 1. Then we should expect on physical grounds 

that this continuous bend is approximately equivalent to its correspond- 
ing discrete tilt for all angles a. However the above analysis seems to 
guarantee this equivalence only for large a, and not in an obvious way 
for small a. 

While it is not obvious, it is true that this equivalence, remains valid 
for small a, or more precisely when (104a) is violated but (104b) re- 
mains true. To study this matter let us consider the case of a guide with 
zero differential attenuation constant, Aa = 0, so that 

Ar = jA@. (105) 

For brevity we restrict our attention to the TE i-spurious mode con- 
version coefficient, which for convenience we call TV (e.g., see the label- 
ing of matrix components in (56)). For the discrete tilt 

T 10 discrete = j sill C t Ct ■ eX P (- ^ ^ ^ A (106) 
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For the continuous bend the results of (1)8) may bo written [using ( 105)] 
as follows: 



■* 10 continuous — J 

J 



.4'feh-] 



( 2C\a 
\A0-A 



exp^-^-^A); (107) 



ff x) = vi + Z 



2C,a 
Xow when .c » 1, /(.r) % 1; thus when , , » 1 [see (104a)] 

then from (IOC) and (107) we have Tu, discrete ~ 7'io continuous , as stated 

above. 

2C ' a , , 

Xext consider the ease where i A/ o i A ' s not largo compared to 1. 

From (10C) and (107) we see that T 10 discrete ~ Ti continuous if 

'(f'l)- f><<1 - <108) 

If , A fl t . g A « 1, then noting that f(x) «- for .r « 1, (108) becomes 

L^1a«1, (109) 

2C,a 

which is simply the condition of (104b). If , , = 1, noting that 

/(l) = V2, (108) becomes 

V'2C,a«l, (110) 

i.e., the conversion coefficient must be moderately small. However a 
more detailed calculation for this case shows that 

0.952 io discrete = J 10 continuous ^ ■* 10 discrete (Ilia) 

for 

C t a ^ 0.55 radian, sin C,a ^ 0.523. (111b) 

The condition of (111b) is far less restrictive than that of (110); (111) 

2C,a 

states that for . "t; — ; = 1, the continuous bend and its corresponding 

I A/3 | • A 

discrete tilt will have conversion coefficients that differ by less than 5 
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per cent provided that the coupling coefficient to the spurious mode is 
less than 0.523, which is quite large. 

These considerations render it plausible that the corresponding dis- 
crete and continuous cases approximate each other well if (104b) is 
satisfied, for both large and small tilts. 

2.3.4 Perturbation Theory for the Coupled Line Equations 

Exact solutions to the coupled line equations, (85) and (86) or (90) 
and (91) [or (93) and (94), (95) and (96)], are known in only a very 
few special cases, i.e., for a few special c(z) . In subsequent work we will 
require a solution for an arbitrary c(z) so that we may eventually treat 
the case of random c(z) . Therefore it is necessary to consider approximate 
solutions to the coupled line equations, as we have done in the discrete 
case. 

These or similar equations arise in the design of a variety of devices 
in which the coupling between two modes is of interest, such as direc- 
tional couplers, tapers or impedance transformers, and bends. In domi- 
nant mode transmission systems the signal mode is the forward TEM 
wave, the spurious mode the reflected TEM wave. 24 ' 26,28,27 Directional 
couplers in a variety of waveguide systems have been similarly stud- 
ied. 17,28 Tapers and bends in a variety of waveguides have been studied; 
several recent papers deal with this type of problem when the principal 
wave is the TEoi mode in circular guide' ' as noted earlier. In much 
of this work an approximate solution is used which is similar or identical 
to the one employed here. This approximate solution has also been used 
to study randomly distributed nonuniformities in ordinary (dominant 

, N . • • ,• 29,30 

mode) transmission lines. 

The approximate solution that we use is most readily found via Pi- 
card's method of successive approximations. '' This is discussed in de- 
tail in a companion paper, 1 where it is shown that the solution to (90) and 
(91) may be written as follows: 

G (z) = £ g 0M (z), (112a) 

»=0 

<?i(«) = E g«nM. di2b) 

n=0 

<7<><n)(z) = j I c(s) e A1 Vi(»-i)(s) ds, n ^ 1. (113a) 

giwiz) = j f c(s) <T A, Wn(s) ds, n ^ 1. (113b) 
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(113c) 



Pino) (2) = 1, 9uoM = 0. 

It is readily seen that 

0o<„)(z) = 0, wodd (114a) 

gw,)(z) = 0, rc even. (114b) 

Bounds on the terms of the series solution of ( 1 12) are given as follows: 



f I c(«) I ds\ 



90.n)(z) 



= 



n\ 



n even. 



n odd. 



(115a) 



f/K«)U, 



= 



< - 



c(s) I ds 



n even. 



n odd. 



(115b) 



11: 



These bounds are in a sense the best possible, since cases are known 
where the equalities in (115a) and (115b) are satisfied. 

Now suppose that the series solutions for G (z) and Gi(z) given in 
(112) and (113) converge so rapidly that only the first nonzero terms 
that depend on the coupling coefficient c(z) need be retained. From (114) 
we see that the n = 1 term of (112a) for G (z) is identically zero; there- 
fore we retain in this equation only the n = 2 term in addition to the 
n = term, which is simply 1. The first nonzero term in (112b) for 
Gi(z) is the n = 1 term. We again use the notation of (63) and (04): 



Go = l — p = gc = e = e 
A = A - j®. 



-V»; 



116) 



Then we have: 



P (z) « f c(s) e M ' s ds [ c(t) e- Ml dt. (117) 

J • (i 

Gi(*)» j f c(s) e'^'ds. (118) 
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Since we have assumed that c(z) is sufficiently small so that \ p\ « 1, 
we have to second order from (65) and (117): 

gttl - Re// (119a) 

A « -In (?„ = // (119b) 

A « -In | Go | = Re// (119c) 

fl«ZG„= -Im // (119d) 
where // is shorthand for 

// = [ c(s) e sr " ds f c(t) e _ * n r// (120a) 

= / c M '"duf c(s)c(s + u)ds (120b) 

•Ml ''(I 

= l f f c(s)c(t) e mi "" ds d*. (120c) 

2 «'o •'o 

These results for the continuous bend are analogous to those of (63) to 
(65) and (71) and (72) for the case of discrete tilts. The continuous off- 
set and diameter change are of course readily handled in the same way, 
and the above relations hold with only minor modifications. In particu- 
lar, for the corresponding solutions to (95) and (96) for the continuous 
offset and diameter change thej in (1 13a) is dropped, thej in (113b) re- 
placed by —1, and the j in (118) replaced by —1 (compare (73) for the 
discrete case) ; the remainder of the equations in the present section arc 
unaltered. (For the continuous diameter change we have neglected the 
variation of AT.) 

These perturbation results for the continuous case may be regarded, 
as in the discrete case, as the first terms of power series expansions in the 
geometric parameter characterizing the deformed guide. Again, it is 
reasonable to assume that if the deformation is sufficiently small, these 
power series will converge so rapidly that their first terms will provide 
valid approximations. Equation (115) gives bounds on the magnitudes 

of the higher terms; for guides short enough so that / | c(s) | ds « .1, 

these bounds converge rapidly. But the fact that the bounds converge 
rapidly does not guarantee that the terms themselves converge rapidly. 
For example, it might be possible for the second term of (112a) to be 
much smaller than its bound, while at the same time the fourth term 
was close to its bound, so that the fourth term would be comparable to 
the second term, even though the bound on the fourth term was much 



KNOWN IMPERFECTION'S 



1075 



Bmaller than the bound on the second term. (On the other hand this 
might not be possible; we just do not know.) Consequently we have no 
precise information on the way in which these perturbation solutions 
approximate the true solutions; the accuracy of these approximate solu- 
tions, which are used throughout the following statistical analysis of 
random guides, must be accepted largely on faith at present. 

In Section 2.3.5 we show how the discrete and continuous perturbation 
theories may be obtained as limiting cases of each other. Section 2.3.6 
discusses a modified perturbation theory, which makes it appear plausible 
that the results of ( 119b), (119c), and (119d) for the TE 01 complex loss 
A, loss .4 , and phase (-) hold true in a wide range of important cases where 
| J'J | » 1, so that the other approximations of (117), (118), and (119a) 
fail. This extension is important in permitting the analysis of long guides. 

2.3.5 Transformation between Discrete and Continuous Perturbation 
Theory 

In this section we consider the relationship between the discrete and 
continuous perturbation theories; we select the case of discrete tilts 
and continuous bends for purposes of illustration. 

Let us first consider the transformation from the discrete to the con- 
tinuous case, as illustrated in Fig. 7(a); but instead of considering a 
single differential section of line with a single discrete mode converter, 
us in Section 2.3.2, we consider the entire line. Let the tilts be equally 
spaced, with spacing Az, and further let the position of the fc th tilt be 
Zk ; thus 

z k = AAz. (121) 

Lei the angle of the k th tilt be a l: . From Fig. 8 we have, as in (82), 

*«£+•••, (122) 

where the dots indicate terms of higher order in Az. R k is defined in Fig. 
8; in the limit as Az —* 0, Ri. approaches the radius of curvature of the 
guide axis at z k . 

Consider first the spurious mode normalized transfer coefficient, Gi . 
From ( 72 ) we have, setting h = Az and neglecting terms of higher order 
than Az, and substituting (122) for a k , 

Cn^jt^Az-e-^-^. (123) 

Setting 
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Zi = iAz = St, Az = As, Re = R(Si), (124) 

(12.3) becomes 

G^jt^-e-^'e^'As. (125) 

,=i lC(Si) 

We now let As — > and N — > °° in such a way that iVAs remains con- 
stant, 

NAs = L, (126) 

where L is the length of line under consideration. Then R(si) — > /r"(s), 
the radius of curvature of the guide axis, e T ' ' — > 1, and the summation 
of (125) becomes the. following integral. 

**"£&!'-"*■ (,27) 

Finally, noting (86) or (91), (127) becomes 

Gtttj f c(s) e-^'ds, (128) 

■'0 

which is identical to the Picard approximate solution for the coupled 
ine equations for the continuous case, given in (118). 

Similarly the approximate solution for p given in (71) for the discrete 
case can be shown to approach the Picard approximation for the con- 
tinuous case, given in (117). Setting k = Az in (71) and substituting 
(122), we have on neglecting higher-order terms in Az 

?~ltfe A*Y + £ E ^£§^ « ir ° W>i! - (129) 

Setting 

Zi = iAz = Si , Azi = As, Ri = R(Si), 

(130) 
Z j = jAz = lj, Azj=At, Rj = R(tj), 

(129) becomes 

We now let. As — > and At — > in such a way that NAs and NAt remain 
constant, 

NAs = L, NAt = L, (132) 

where L is again the total length of line. Then R(Si) — > R(s), R(U) — * 
R(t), where R(s) and R(t) are the radius of curvature of the guide axis. 
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The first term of (131) (the single summation) approaches zero; the 
second, however, approaches a double integral to yield 

Using (86) or (91), and further interchanging the order of integration 
and the integration variables, 

p « [ L c(s) e AVt ds f c(t) e _Ari dt, (134) 

''U ■'0 

which is identical to the Picard approximation for the solution to the 
coupled line equations, given in (117). From (120) we see that the two 
other forms of (7 Id) also have their equivalents in the continuous case. 
Next let us reverse the above process, and go from the continuous 
approximate solutions back to the discrete ones. Consider a line of 
length Nlo with N equally spaced tilts a distance /„ apart. The k th tilt 
is located at z k , where 

z k = (k - 1)1,, (135) 

and has an angle a k . (The form of (135) was chosen to be consistent 
with the notation of Fig. 6 and Section 2.2.) From (97b) or (102b) we 
may write the continuous coupling coefficient for this case as 

c(z) = C t jt, aic -8[z- (k - 1)M- (136) 

We now substitute (136) into the approximate solution for the con- 
tinuous case, and derive the approximate solution for the discrete case. 
Substituting (136) into (118), we find 

GitJW^ECW-^ 1 '-"' , (137) 

i=i 

in agreement, with (72c) and (72d). Substituting (136) into (117), we 
have 



p(Nk) 






I 



C.Eo/'*- (J - \)U\e- Mt dt 
o ;=i 



= [ °C t Z«.-5[.s - (» - l)l ]e* r 'ds 



(138) 



N 



■2 C t aj-ll{s - (j - l)/o] e 

;=i 



-M'(j-\)i 
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where it(.r) is the unit step function, 

1, x > 
u(x) = (139) 

0, x < 0. 

Continuing, 

p(OTo) « EC 1 a/ r,H, "f I C j a/ r,H "« + if; (Ca,)'. (140) 

»-2 1-1 i=l 

The second term of (140) (the single summation) was obtained via the 
relation 



/ 



6 

f(x)u(x)8(x) dx = |/(0), a < < b. (141) 



Alternately, the result of (140) may be obtained by regarding the 
5-funetions as the limit of some continuous functions [e.g., (07a)], 
taking the limit after the integrations have been performed. Finally, 
interchanging the order of summation and the summation indices, (140) 
becomes 

p(Nh) = \ E (C t ai) 2 + £ E (CtaiKCta,)*™**-*, (142) 

i-1 i=l j=i + l 

which agrees with (7 Id) for the discrete case. We see that the single 
summation, which ''disappeared" in going from the discrete to the 
continuous case, has satisfactorily "reappeared." The alternate forms 
for p given in (120) may similarly be transformed to their discrete 
equivalents in (71d) via the substitution of (13(5) and appropriate 
manipulation of 8-functions. 

The corresponding analysis for offsets and diameter changes is readily 
performed. 

One consequence of the results of the present section is that bounds 
derived for the approximate solution in the continuous case may be 
directly applied to the approximate solutions given in Section II for 
the discrete case. 

2.3.6 Logarithmic Form of the Coupled Line Equations, and Improved 
A pproximate Solution 

The perturbation results of Section 2.3.4, given in ( 1 16) to (120), were 
expected on intuitive grounds to be valid for short lines, whose additional 
loss due to mode conversion is small. In particular, the results of (117) 
and (119) depend on having | p | « 1, | // | « 1. These relations may 
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again be regarded as the first terms in power series expansions for the 
various quantities. 

As the length of guide increases we expect the mode conversion loss 
to increase; for example in Sections III and IV below it is shown that 
for all random guides that we deal with, < Re J J > is proportional to the 
total length of the guide. It is clear that as the length of guide increases, 
eventually | // | » 1. Under these conditions we no longer obtain a 
valid approximation from the first two terms of (112a); many terms 
become significant, and (117) and (118) are no longer valid. It would 
seem that the results of (119) are also invalid for long guides. However 
while (119a) is certainly invalid, wc will see that it is plausible that 
(119b), (119c), and (119d) will remain good approximations for long 
guides with large mode conversion loss if the differential attenuation 
constant is large enough in a certain sense compared to the coupling 
coefficient. The detailed mathematical analysis for this problem is given 
in a companion paper; 18 however, it is not difficult to see on physical 
grounds that something of this sort is to be expected. 

Consider a long guide of length L with a large enough differential 
loss so that the total differential loss | Aa \ z in a short section of length 
z is large; 

|Aa|z»l, z«L. (143) 

Now let this guide be divided up into M equal sections of length z by 
ideal mode filters, so that 

Mz = L. (144) 

An ideal mode filter by definition transmits TEoi without loss or phase 
shift, and suppresses all other modes completely. (Practical mode filters 
may consist of a section of helix guide that has low loss to TE W , high 
loss to all other spurious modes other than the higher order TE 0m .) In 
addition to the requirement of (143), further assume that each section 
z is short enough so that the perturbation results of (116) to (120) do 
apply to the individual sections; for example, we might require from 
(115) that 

f" \c(s)\ds«l. (145) 

It is more or less obvious on physical grounds that under these condi- 
tions the insertion of ideal mode filters will not alter the over-all TEoi 
loss significantly, because the spurious mode level is not likely to build 
up to a significant magnitude. 
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Denote quantities applying to the k tb section by the symbol k in the 
following equation; quantities applying to the over-all guide will be 
written without subscripts. Then the over-all TE i transmission param- 
eters for the guide with ideal mode filters are given in terms of the TE i 
parameters for individual sections as follows : 



A 







M 

= n^ , 

fc=i 


M 

= Ukg, 


1/ 

= I>, 

fc-=i 


= E^, 

k=l 


M 



(146) 



Then, for example, A may be written from (119), (120) and (146) 
as follows : 

M rkz pa 

A « £ / da dt-c(s)c(t) e m "- l) . (147) 

fc=l J(k-l)z J (k-l)z 

In (147) the double integral has been taken from (120a); the other two 
versions of this integral could of course be used equally well. Now it is 
more or less obvious by inspection of the integrand and the limits and 
from physical considerations that (147) is approximately equivalent to 

A « ( ' \ls f ' dt-c(s)c(l) e AV(s -'\ (148) 

In other words, the results of (119b), (119c), and (119d) should remain 
valid for the whole line, as stated above. Since the requirement in the 
above crude argument is that the mode conversion loss must be small in 
a section for which the total differential loss is high, we would intuitively 
expect that for (119b), (119c), and (119d) to hold for a long line, the 
ratio of the coupling coefficient \ c(z) \ to the differential attenuation 
constant | Aa | must be small in some sense which wo have not yet 
attempted to define. 

While the above argument may be physical!}' appealing, it is ob- 
viously desirable to put these statements on firmer ground; this is done 
elsewhere in this issue. 1 We summarize here briefly some of the results 
of this investigation. 
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Under certain conditions we may write A (2) as a series expansion; 
A(z) = f c(s) e Ar ° ds f c(t) e _Ari dt + £ \ n (z). (149) 

•'O •'0 n=2 

In (149) we have written out only the first term explicitly; we see that 
this is identical to the previous approximation of (119b) and (120a). 
Bounds on the higher terms are given in detail in Ret*. 18; we give here 
some slightly poorer but simpler bounds. 

I Xi(a) I = f c(s) e Ars ds f c(t) e~ ATt dt g K [ | c(s) | ds; 

I -X.(f) I ^ Zv 3 (2.22;>/i 2 ) n - 2 f | c(s) \ ds; n ^ 2, (150) 

•'o 

^ K ^ 0.3; 

where 7v is defined by 

( I c(s) I e Aa{z ' s) ds ^ K for every 2^0. (151) 

Finally, convergence is guaranteed only for K < 0.455; a case is known 
where the series of (149) diverges for K > 0.5. Equation (151) requires 
that I c(z) I be uniformly small in a certain sense, with respect to | Aa |. 

Once again we take as an approximate solution the first term of (149). 
The bounds in Ref. 18 are again almost the best possible in the same 
sense as in the case of Section 2.3.4; i.e., cases are known where the 
higher-order terms are almost as large as their bounds. In cases of inter- 
est K « 1, and the bounds converge rapidly. However, as in the case 
of Section 2.3.4, this is not sufficient to guarantee that the terms them- 
selves converge rapidly; and here again we lack precise information on 
the way in which the first term (perturbation solution) approximates 
the true solution. 

Because of the relationship between discrete and continuous cases, 
similar statements can be made regarding the results of (65) and (71) 
for the discrete case. 

2.3.7 TEai Loss in Terms of Fourier Coefficients of c(z) when Aa = 

If the differential attenuation constant is equal to zero, 

Aa = 0, (152) 

the above perturbation results for the loss A [see (119c) and (120)] 
may be further simplified. This case is of interest as an approximation 
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to the .situation present in real copper guides. In practical copper guides 
mode niters, consisting of a length of helix guide, must be inserted at 
regular intervals, perhaps a few hundred feet apart, for reasons that 
will be further discussed in Sections III and IV below. If the mode 
filters may be assumed ideal, we consider each section separately and 
simply add their individual TE i loss and phase, i.e., A, A, and 0, as in 
(140). If for each copper guide section of length L the total differential 
loss | Aa | L is small, 



Aa | L « 1 , 



(153) 



on intuitive grounds we approximate the true solution by setting Aa = 
in the various approximate solutions of (116) to (120). 
Thus, let us set Aa = in (119c) and (120c) to obtain 

A « iRe f f c(s)c(t) e^ 1 " 1 ds dt (154) 

2 Jo J o 

for a guide of length L. Recalling that c(s) and c(t) are real in ideal 
metallic guide, (154) yields 

A w I Re [ L f L c(s)c(t) e jmt - 8) ds dt 
2 ^o ""o 

= - [ L r c(s)c(t) e im, - 8) ds dt 
2 •'o *'o 

= \ I" c(s)e- JA **ds f c(t) 
2 •'o J o 

1 I C L 
= l\ c(s) 



(155) 



„-?A8« 



.+iAfl< 



dt 



ds 



Summarizing, 



A 



1 I f L 
III c(s) 



-iAS* 



ds 



Aa = 0. 



We note further from (119a) and (116) that 

L 



| Go | = g » 1 - 1 \[ c(s) 

From (157) we have to second order 



„-jA0« 



ds 



<?« 



- / c(s) 



n -i*P* 



ds 



(156) 



(157) 



(158) 
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But from (118), for Aa = 



Gi 



[ c(s) e-^'ds 



Therefore from (158) and (159) we have to second order 
I G I 2 + | G, f « 1 , Aa = 0. 



(159) 



100) 



In Appendix C it is shown that (160) must hold exactly for Aa = 
[see (C-13)]. This is something like conservation of energy; in fact if 
a = ai = 0, (C-13) is precisely conservation of energy. 

We note in passing that a similar result to that of (156) is readily 
found for the discrete case of Section II. Proceeding in an analogous 
way from (65c) and the second form of (71d), we find for the discrete 
case 



.1 



£ « 



-}&Plai 



Aa = 0. 



101 



Similarly the result of (160) is readily seen to hold true to second order, 
and by Appendix C must also be true exactly. 

Equation (156) states that the loss in nepers A is one half the square 
of the magnitude of the Fourier transform of the coupling coefficient, 
with transform variable A/3, the differential propagation constant. Since 
we deal with the case Aa = 0, the logarithmic perturbation theory of 
Section 2.3.6 does not indicate that (156) is valid for long lines; conse- 
quently, this approximate result will remain valid only for short line 

sections, perhaps subject to a condition such as / | c(s) | ds <C 1. 

It will prove convenient to rewrite (156) in a slightly modified form. 
Define 



/ = 



r(A0L/ 



■J ii 



c(s) e 



-j&Pe 



ds. 



We have 



n = 



c(s) e~^ s ds 



(162) 



(163) 



/ is closely related to the approximate expression for the spurious mode 
transfer coefficient, given in (118). From (162) and (163), (156) may 
be written as 



Atthll 



(164) 
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It turns out to be useful to expand the coupling coefficient c{z) in a 
Fourier series as follows :f 

c(z) = £ c„ e j2 " i:!L ; c_„ = c*. (165) 

Jt=— 00 

Substituting (165) into (162), we have after some algebraic simplifica- 
tion : 

M/3L \ 

sin ir ( — n 1 

/ = l £ c„(-ir — At [ — r^- ( 166 : 

n=-=o /A/3L 

Now, (166) has a rather striking form. Assume for a moment that 
the only variation of / (and hence A) with the frequency/ of the radia- 
tion in the guide occurs through A/3. If we take the independent variable 
in (166) to be proportional to A/3, then (166) is simply the sampling 
theorem representation of a complex band-limited function, 33,3 i.e., 
a complex function whose real and imaginary parts are each band- 
limited. Taking the dimensionless quantity — — as the independent 

variable, I[-^— ) will contain no frequencies | v | greater than \.% By 

(164) the loss A is proportional to the square of the magnitude of /, 
or alternately to the sum of the squares of the magnitudes of the real 
and imaginary parts of /; therefore the TE i loss A, regarded as a func- 
tion of the normalized independent variable -=— , will contain no fre- 
quencies greater than 1. If for the time being we neglect any variation 
of the c„'s with A/3 in (166), then il-^—J is determined by its values 
at the sample points. At the n th sample point 



t Here c„ is the n th complex Fourier coefficient of the continuous coupling 
coefficient c(z). In contrast, in Section 2.2 above we have used c* to represent the 
coupling coefficient of the k tb discrete mode converter. In the following work the 
meaning will always be clear from the context. 

t Here, and often in the sequel, we use the word frequency to denote the inde- 
pendent variable v of the Fourier transform of some quantity of interest. In the 

present case we consider the Fourier transform oil f -5 — 1 with respect to the 

independent variable -Jj—; since -i— is dimensionless, the corresponding Fourier 

Zir Air 

transform variable v is also dimensionless. Thus, if £T(*') is the Fourier transform 

of / (^), *M = / Jft) «-«"««■ 
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A^ = *.?*„£. (167) 

2tt ' A/3 n 

The n th sample point occurs at the free-space wavelength or frequency 
(of (he radiation in the guide) where the beat wavelength B is equal to 

the mechanical wavelength - associated with the n th Fourier coefficient 

n 

of the coupling coefficient. The value of / at the n* b sample point is 

I{n) = LcX-l)". (168) 

By (1G2) and (118) the spurious mode transfer coefficient at the ?i th 
sample point for the continuous bend becomes 

Gt(n) ttjLc n , (169) 

with an analogous result for the continuous offset and diameter change. 
I(n) and G\(n) are determined only by the w th Fourier coefficient. At 

intermediate values of -— -, I [ — — ) is determined by interpolating 
2r \ It / 

Sill 7" 

between the sample values with ' — — functions, as shown in (1(56). 

Now the object of the present calculations is to determine the loss A 
as a function of the frequency /, and later to determine the statistics of 
the loss-frequency curves for guides with random imperfections. How- 
ever, A/3 is approximately proportional to the free-space wavelength X 
of the radiation in the guide, and thus inversely proportional to the 
frequency/, if the two modes involved (the TE i signal mode and the 
spurious mode) are both far from cutoff. Therefore for analytical pur- 
poses it is more convenient to choose an independent variable propor- 
tional to the free-space wavelength X rather than to the frequency /, 
and this is what is done in the present paper. If a single spurious mode 
is under consideration, it is most convenient to choose the dimensionless 

parameter -^— , which is approximately proportional to X, as the inde- 
2t 

pendent variable. If more than one spurious mode is being considered, 
we see in Section 2.4 below that the loss A is given by a sum of terms 
of the form given in (156), with of course different c's and A/3's for each 
spurious mode. While the different A/3's are all approximately propor- 
tional to X, they have different constants of proportionality, and hence 
it is perhaps most convenient to take X itself as the independent variable. 
In practice we will consider only fairly narrow percentage bandwidths 
(although the absolute bandwidths will be enormous compared with 
conventional communication channels) at any one time. Therefore the 
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fact that we take X as the independent variable instead of / should 
cause no serious inconvenience. 

The sampling theorem interpretation of (16G) requires that the c„ , 
defined by (165), be independent of X (or frequency) and hence A/3; 
this may or may not be true. Since c(z) is given by the product of one of 
Morgan's coupling coefficients and a geometric parameter [see (86) or 
(94)], the X-dependence of the c„'s is identical to the X-dependence of 
C t , Co , or Cd • From the equations for these three coefficients given in 
Appendix A we see that far from cutoff, C and Cd are approximately 
independent of X, as desired, but that C t varies approximately inversely 
with X. In Section 2.3.8 below we shall introduce additional coupling 
coefficients of Morgan, S [n „,] , which permit a similar treatment for 
general continuous cross-sectional deformations of copper guide; the 
S[„m] vary approximately directly with X. The geometric parameters 
associated with C a and C d are the derivatives with respect to distance z 
of offset x and radius a respectively; hence these geometric parameters 
are dimensionless. The geometric parameter associated with Ct is curva- 
ture, which has the dimension of length -1 . The geometric parameters 
associated with the E [nw] all have the dimensions of length. As a con- 
venient device for recalling these facts, the exponent of X in stating the 
X-dependence of the coupling coefficient is the same as the exponent of 
length in stating the dimensions of its associated geometric parameter. 

We wish to retain the sampling theorem interpretation of (166) even 
in those cases where the c„'s are not independent of X and hence of A/3. 
Now in those cases where the. coupling coefficients are not approximately 
constant (i.e., C t , S[ nm j) the variation with X is quite slow. Since as 
mentioned above we need consider only narrow percentage bandwidths, 
the principal variation with X in (166) occurs through A/3, and not 
through the coupling coefficient. Consequently, over the moderate band- 
widths of interest we may evaluate Morgan's coefficients at the middle 
of the narrow band under consideration. 

There are several more elegant methods for deriving this approxima- 
tion that state in effect that c n is to be evaluated at the wavelength X„ 
(or frequency) corresponding to the n th sample point, given in (167), 
rather than at the operating wavelength, as implied in (166). This is 
certainly true at the sample points, by (168), and appears plausible in 
general. In Section 2.3.9 below we show, for example, that under certain 
reasonable conditions a guide with a given straightness deviation may 
be described equally well as either a continuous bend, for which the c„'s 
vary approximately inversely with X, or as a continuous offset, for which 
the c a 's arc approximately independent of X. By a suitable transforma- 
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tion it is always possible to change the geometric parameter to dimen- 
sionless form (in general, to the derivative with respect to z of a param- 
eter with the dimensions of length, such as radius, offset, ellipticity, etc.) 
so that the corresponding coupling coefficient is approximately inde- 
pendent of X (and hence A/3) far from cutoff. 

Alternately, this result may be derived by fairly simple manipulation 
of the sampling theorem. 

The final result, however, is that over the fairly narrow percentage 
bands of interest we are justified in neglecting the X-(or frequency) 
variation of the coupling coefficient (i.e., C t or E[„ m ]). 

The TEui loss A in terms of the Fourier coefficients c„ is, from (1G4) 
and (166), 



j 2 oo oo 

A «-±- X) £ c m c,*{-\) 



(m-n) 



/aql \ . /a/?l 

sin it I — m I sunr I — — — n 

\2* / \2x /_ (170) 



/A/3L \ /A0L 



At the sample points of /, given in (167), the TE„i loss becomes simply 

A(n)tt^\c n \\ (171) 

As mentioned earlier the bandwidth of A is twice the bandwidth of /; 
thus the Fourier transform of A with respect to the independent variable 

— — will contain no frequencies greater than 1. 

Therefore while A is also a band-limited function, it has twice as 
many sample points as /. Consequently A is not completely determined 
by its values given in ( 171 ) at the sample points of /. Simple examples 
are readily found of two different guides in which the TEoi loss is the 
same at the sample points of /, given in (167), but differs greatly be- 
tween these sample points. Thus consider the following two cases, in 
which all but two adjacent Fourier components of the coupling coeffi- 
cient are identicallv zero: 



(a) c(z) = 2 cos 2irk- \ + 2 cos 27r(A- + 1) \ . 

Li Li 

(b) c{z) = 2 cos 2irfc \ - 2 cos 2ir(k + 1) % . 

Li Li 



(172) 
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2.4 



* A 



c(z) = 2COS277k-£- + 2COS277(k + l)|- 
C(z ) = ZCOS2TTK-- -2COS 277(K+l) -^- 



k-4 




k+5 



Fig. 10 — TEoi loss for coupling coefficient with two sinusoidal components. 

The TEoi loss A is shown for these two cases in Fig. 10. While the losses 
are the same at the sample points of /, they differ markedly in between. 

While A is not completely determined by its values at the sample 
points of /, it is clear from (170) that the principal contribution to A in 
the region near the A- th sample point arises from values of m and n close 
to k. This relation between the TEoi loss and the coupling coefficient is 
of great importance, as it means that the additional TE M loss due to 
some particular spurious mode in a given frequency range is caused 
primarily by the components of the coupling coefficient which lie. in a 
corresponding (mechanical) spectral region. Thus, for example, loss due 
to TEoi — TE12 coupling in 2-inch guide between 50 and 60 kmc is 
caused primarily by straightness deviations with periods between 2 and 
2.4 feet, the range of beat wavelength B corresponding to this frequency 
band. 

Equation (170) is also useful in considering situations where only 
isolated Fourier coefficients of the coupling coefficient are nonzero; the 
TEoi loss A then becomes a series of isolated narrow peaks of shape 



sin .r 
x 



. There are at least two such cases of practical interest : 



i. A periodically supported copper guide that deforms elastically under 
its own weight (the "serpentine bend"), 3b with support period short com- 
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pared to the total length of guide between mode filters. The Fourier com- 
ponents occur only at the fundamental and harmonics of the support 
period; the amplitude of the Fourier coefficients decreases rapidly as 
the order of the harmonic increases. Thus the TEoi loss vs wavelength 
curve will have a series of equally spaced narrow loss peaks of rapidly 
decreasing magnitude. 

ii. Shuttle pulse measurements in a copper waveguide without mode 
Jitters. Here, because of the absence of mode filters, the pulse traverses 
an iterated structure. Since the Fourier series expansion of the coupling 
coefficient for any number of round trips must be identical to the ex- 
pansion for a single round I rip, for X round (rips only every X th Fourier 
coefficient will be nonzero. Again the TE i transmission will consist of 
equally spaced narrow loss peaks separated by wide pass bands of low 
loss. 

These two cases will not be discussed further in the present paper. 

The spacing between sample points of I is an important parameter in 
the analysis and in the interpretation of experimental transmission data. 
Since A/3 is proportional to the free-space wavelength X, 

M3 = DX, (173) 

where the constant D depends on the spurious mode. From (167) the 
wavelength corresponding to the //"' sample point, X„ , is given by 

D\ n L . 2im ,_.. 

— — = n, \„ = _ . (1/4) 

sir UL 

Therefore the sample point spacing in wavelength is 

AX„ = X„ + , - \„ = ^. (175) 

Thus from (173) and (175), noting the definition of beat wavelength B 
in (167), 

f-i- (176) 

where we have dropped the subscript n since the result of (175) is in- 
dependent of n. Since 



T~7' (177) 



we have finally 
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*/«£-? (178) 

for the sample point spacing, in either frequency or wavelength. We see 
from (175) that the sample points are equally spaced in wavelength, and 
consequently unequally spaced in frequency. However over a sufficiently 
narrow band the variation in sample point spacing in frequency will be 
small. The beat wavelength B between TEoi and all other TE modes is 
tabulated in Appendix D for the frequencies 50, 55, and 00 kme in 
2-inch diameter guide. 

2.3.8 Morgan's Coupling Coefficients for Small Cross- Sectional Deforma- 
tions in Lossless Metallic Guide 

Morgan has determined to first order the spurious modes excited by an 
incident TE i wave at a small arbitrary deformation of the cross section 
of a cylindrical metallic guide. This analysis must include the continu- 
ous offset and diameter change studied above when those deformations 
are small; it includes in addition higher-order deformations of the cross 
section, such as ellipticity, tri-foil, etc. 

Let the surface of the slightly deformed guide be given in cylindrical 
co-ordinates by 

r = a + p(<p,z), (179) 

CO 

p(<P,z) = £ [«»(*) eosnv? + b„(z) sin tup], b»(z) = 0, (180) 

z ^ 
pfot) = 0, a n {z) = 0, b„(z) = 0; (181) 

z ^ L 

where a is the radius of the perfect guide. p(<p,z) must be suitably small 
(we have omitted Morgan's small dimensionless parameter e). Then 
the n = term corresponds to a continuous diameter variation, the 
n = 1 term to a continuous offset, the n = 2 term to ellipticity, the 
n = 3 term to what has been called tri-foil, etc. By (181 ) it is assumed 
that the guide is distorted only in the interval < z < L. 

Morgan has shown that to first order a TE„i wave incident on the 
deformation of (179) to (181) excites the forward and backward TE, im 
modes. We denote the various modes as follows: 



7[om] + — forward TE„„, 

Z[o m ]~ — backward TE ( j,„ 
7[nm]" + , f[„m\ ±+ — two polarizations of forward TE nm 
7 ( „ m] ll_ , I[ nm ] X ~ — two polarizations of backward TE„, 



(182; 
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We denote the two orthogonal polarizations of each mode (except the 
TEom) by the symbols || and _i_, rather than the subscripts x and y used 
in Section 2.1 above. The geometric imperfections of Section 2.1 — offsets 
and tilts — could naturally be resolved into components along the x- 
and //-axes; with the more general deformations of ( 179) to (181 ) this is of 
course not. possible, and the |l and _|_ notation appears more natural. 
In the special case n = 1, using the geometry of Pig. 5, we may if we 
wish identify a x (z) with .r(z), b x (z) with y(z), and the || and j_ polari- 
zations of the resulting TE lm spurious modes with the subscripts x and 
II respectively of Section 2.1. 

We normalize the complex mode amplitudes in the usual way [see 
(88) and (89)]: 

/[0„] ± (Z) = C Tr """ F (7 (0mI ± (2) 

W*(«) = « :Fri - ,- Ci.., ,± (*) (183) 

/, ,**(«) - e Tr *-"'-(7 [ „ ml i± ( 2 ). 

In lossless metallic guide the propagation constants r [nm ] are of course 
pure imaginary, 

r,,,,,,, =^ [nm] . (184) 

Assuming a unit incident, TEoi wave, 

G m + (0) = 1, (185) 

the (7[„,„| become the normalized spurious mode complex transfer co- 
efficients. They are found to first order from Morgan's analysis, and 
given in a slightly modified notation as follows: 

— J3[„ OT ] 






'(L). 



a n (z) 



-Mov&nm)* 



where the coupling coefficient E[„,„] is given by 

1 k 0l k nm a 1 



— [mill 



e„ « y/k Hn ? - li 1 s/pmpin 



dz (186) 



(187) 



n = 

n > 1 



/.'(ft*) = 0. 



(188) 



It is worth noting that the angular index n of the spurious mode is 
identical to the subscript of the pertinent component of mechanical de- 
formation, o„ or b„ ; therefore, to first order the modes excited by TE ( ,i 
at a given deformation have the same angular symmetry as the defor- 
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mation. E[„„,] is the same for forward and backward modes, and for 
both polarizations; in contrast C ,C t , and C d of Appendix A are different 
for forward and backward modes. This contrast in behavior will be dis- 
cussed in Section 2.3.9. The S [nm] are given in Appendix D for the fre- 
quencies 50, 55, and 60 kmc, in 2-inch diameter guide. We note again 
that these results hold true only when there is no distortion of the guide 
at the ends, z = and z = L, by (181). 

Special comment is necessary for the G [0m j terms in (186), and in 
particular for the G m + term. The n = terms correspond to the TE 0m 
modes, which have only a single polarization. Since b Q (z) = in (181) , 
the || terms of (186) are the significant ones; and since we have only one 
polarization, the symbol || may be conveniently dropped. The G m] term, 
corresponding to the forward TEoi , requires special interpretation. It is 
to be considered a first-order correction to be added to the unperturbed 
solution, i.e., G m + = Litis readily seen that G m + in (186) represents 
a phase shift of the unperturbed TE i wave, caused by a change in the 
average diameter of the guide. 

By analogy with the previous analysis, we may use the results of 
(186) to (188) in any of the results of Section 2.3 above for the con- 
tinuous case— e.g., (116) to (120)— by substituting -S^a^z) or 
— Zi„m]b n (z) for the continuous coupling coefficient c(z). While these 
coupling coefficients were derived by Morgan for lossless metallic guide, 
we expect as before that these results provide a satisfactory approxima- 
tion for real copper guide by modifying the various propagation constants 
to take account of the small losses present in copper guide. The relation- 
ship of these results to previous ones for the continuous case is discussed 
in Section 2.3.9. 

2.3.9 Relationships between Various Metallic Guide Coupling Coefficients 

A sufficiently small straightness deviation of the guide axis may ob- 
viously be described equally well as a continuous bend, a continuous 
offset, or simply as a continuous displacement. For these three cases we 
have given above three different coupling coefficients — C t [m) , C [ m ] , 
and H[i m) — that yield, among other things, the first-order spurious 
mode transfer coefficient of a deformed guide. While the results of Sec- 
tion 2.3.4 — in particular (118) for continuous bends and its analog for 
continuous offsets and diameter changes — were derived for only a 
single forward spurious mode, the result for the first-order spurious 
mode transfer coefficient holds true in general, with of course the proper 
coupling coefficient for each of the spurious modes, as discussed in Sec- 
tion 2.4 below. Since a given deformation that is suitably small may be 
described in these three different ways, and since we must get the same 
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answer for the spurious mode output in each calculation, the three cou- 
pling coefficients must be related to each other. Similar statements ob- 
viously apply to the two coupling coefficients Cd [ml and E[om] for diameter 
variations. 

For a guide with a continuous bend, confined to the x-z plane, the 
TEi,,,* transfer coefficient is, from (118) and (91), to first order 

<W «i f £rC *-""*»» dz, (189) 

where H x (z) is the radius of curvature of the guide axis, and the super- 
script x indicates that the .r-polarization is under consideration. If the 
slope of the guide axis is small compared to unity, the distance along 
the perturbed guide axis is approximately equal to the distance along 
the unperturbed guide axis, so that we may regard z in (189) as meas- 
ured along a fixed rectangular co-ordinate axis. Further we may approxi- 
mate the curvature as 

=\- K *"(*). (190) 

Then (189) may be written 

GW - " ~ JC<l»>* r **(«)<T J(fc,:ff, - ta dz. (191) 

Similarly for a guide with a continuous offset in the x-z plane the spuri- 
ous mode output is, from (96) and (118) with the,/ replaced by — 1, 

AW* ~ -aw* [ L x'(z)e- i( ^ 0im)! dz. (192) 

•'o 

If the ends of the guide are parallel, so that 

x'(0) = x'(L) = 0, (193) 

:iii(l if the perturbation is small, the spurious mode output must be 
identical in the two cases and the right-hand sides of (191) and (192) 
must be equal. Integrating ( 191 ) by parts, we find subject to (193) 

<W* % -(ft* =F ftjCiw* f x'(z)e- i(Po ^ i "' )z dz. (194) 
For ( 192) and ( 194) to be identical, we must have 

= /3oi=F/3i».. (195) 






Using (A- 1 ) and (A-2) of Appendix A, (19.5) is easily shown to be an 
identity. 
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Similarly a line with a continuous offset, but with no net offset be- 
tween its two ends, is described in terms of displacement of the guide 
axis by the results of (186). Thus setting 

ai (z) = .viz), (196) 

we have for the spurious mode output 

<W* « -JKd»i f x(z)e- m ^^' dz, ( 197) 

where we have replaced the superscript || in (186) by x. This result is 
valid only when the distortion vanishes at the ends of the guide; in this 
case this means that there must be no offset between the two ends of the 
guide, 

;t(0) = x(L) = 0. (198) 

Subject to (198), (197) must agree with the previous two results. Inte- 
grating (192) by parts, and using the condition of (198), 

(W* ~ -i(A» =F ft-)?*.* f x{z)e- m ™^' dz. (199) 

■mi 

For (197) and (199) to be identical, we must have 



Hi'"] 



'c- (JmTA.. (200) 



By (A-l) of Appendix A and (187), (200) is easily seen to be an iden- 
tity. 

An analogous study can be made for a continuous diameter change. 
From Sections 2.3.2 and 2.3.4 

(W = -<W f r'(z)e- i( ^°>» ): dz. (201) 

•Ml 

From (186) we have 

GWj* = - jK M f a,(z)e- j ^ ^ dz, (202) 

•'o 

valid when 

ao(0) = a (L) = 0. (203) 

Subject to (203), (201) and (202) must be identical. From (179) we 
have r'(z) = a '(z); substituting this into (201), integrating by parts, 
and using (203), 
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Ota* = -j(fa =F jSb»)CW f ao(z)e- mi:ffitm) ' dz. (204) 



For (202) and (204) to he identical, we must have 

S[0m 



C d , m ,± 



= /3..i =F/3« W . (205) 



Equation (A-4) of Appendix A and (187) show that (20")) is an identity. 

It is interesting that coupling coefficients found in such diverse ways 
are so simply related, and pleasing that these different calculations are 
all consistent with each other. 

A simple physical interpretation may be given for the fact that the 
mode conversion caused by a given type of deformation of the guide 
may be described by different coupling coefficients (with of course 
corresponding different geometric parameters). These different descrip- 
tions correspond to using different sets of normal modes to describe the 
fields within the deformed guide. 

For purposes of illustration we take the case of deviation from straight- 
ness of the guide axis. The coupling coefficient H[i„,] and (197) corre- 
spond to the normal modes of the original undeformcd guide. The cou- 
pling coefficient f 0[ „ (] and (192) correspond to the normal modes of the 
deformed guide with cross section perpendicular to the axis of the un- 
deformed guide. Finally, the coupling coefficient C t [ m ) and (191) corre- 
spond to the normal modes of the deformed guide with cross section 
perpendicular to the axis of the deformed guide. This process may be 
carried one step further, using the normal modes of a curved circular 
guide.'"' 

A similar discussion can of course be given for the two representations 
of mode conversion caused by changes in diameter [C,i\ m ) in (201) and 
£,„„,, in (202)1. 

2.4 Extension to Many Spurious Modes and Two Polarizations 

Most of the above results are readily extended to include additional 
spurious modes, coupled to the TE i signal mode to first order, and two 
polarizations of all such spurious modes (except TE 0m ). First, we may 
proceed via the scattering matrices for discrete mode converters includ- 
ing many spurious modes, and thence to the continuous case via limit- 
ing processes, as above. Alternately, we may take the generalized teleg- 
raphist's equations, including all modes, and find perturbation solutions 
as above, and thence treat the discrete case by limiting processes. The 
conclusions of such a study may be summarized as follows: 
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1. The first-order spurious mode transfer coefficients for each of the 

spurious modes are given as above [e.g., (72), (73), (118) or (186)], 
with of course the appropriate coupling coefficient and differential propa- 
gation constant for each of the modes. 

2. For the second-order TE i transmission parameters (Go , g , A , A , 
0) we must replace each term involving the coupling coefficient by a 
summation of similar terms, each corresponding to a particular spurious 
mode, with of course the appropriate coupling coefficient and differential 
propagation constant for each mode. For this purpose we regard the 
two orthogonal polarizations of a spurious mode as two distinct (but 
degenerate) spurious modes. Thus for example we write 

a«Sa M) 4~Z^, e«Ee«. (206) 

The quantities for the m th mode are given for example by (G5) and (71 ) 
or (161) for the discrete case, and (116) and (117), (119) and (120), or 
(156), (170) for the continuous case, with in each of these equations 
the subscript m attached to the coupling coefficients and differential 
propagation constants. 

Much of this analysis is a straightforward extension of the two-mode 
analysis above, and will not be discussed here. The study of convergence 
of the approximate solutions, analogous to that discussed in Sections 
2.3.4 and 2.3.6 and treated in Ref. 18, appears more difficult in the 
general multimode case, and little work has been done. Of course even 
in the simpler two-mode case we lack precise information on the validity 
of the approximate solutions, as discussed above; the same is certainly 
true here. 

From the results of Section 2.3.7 and (206), the contribution to the 
TEoi loss A of each of the spurious modes arises from a different portion 
of the (mechanical) spectrum of the geometric imperfection. For ex- 
ample, consider a guide whose only imperfection is straightness devia- 
tion; the most important spurious modes arc the forward TEi,„ . Con- 
sider the frequency band from 50-60 kmc in 2-inch I.I), waveguide. The 
beat wavelength range, which is equal to the range of mechanical wave- 
lengths of straightness deviation that contribute significantly to the 
TEoi loss component, is shown in Table I for each of these modes (see 
Appendix D). 

The coupling coefficients for each polarization of each spurious mode 
for discrete tilts and offsets are obtained from (37), (38) and (47) 
respectively; for discrete diameter changes the coupling coefficients are 
obtained from (51). Thus in (71e) each of the C"s should have a sub- 
script denoting the spurious mode, and a,- and 6, should have subscripts 
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Table I 



III 


Mechanical Wavelength Range for 5(1-60 kmc Band in Feet 


1 

2 
3 

4 
5 

6 
7 
8 


2.45264-2.04898 
1.99720-2.40812 
0.46389-0.56232 
0.21503-0.26289 
0.12293-0.15228 
0.07762-0.09815 
0.05134-0.06723 
0.03377-0.04763 



denoting the x or y components of tilt or offset, depending on the polari- 
zation of the spurious mode. 

For the cross-sectional deformations of Section 2.3.8, the continuous 
coupling coefficients for the two polarizations of the TE„„, mode are 

c { ]'*(*) = -Bb.io.Ca), (207a) 

ci i ±:fc (2) = -E[„ M ,&„(2), (207b) 

in terms of the geometric parameters of (180). The corresponding differ- 
ential propagation constants are for both cases 



ATi,,,,,]" — Tim] ^F T\„, 



;208) 



For the coefficients C o[m] for continuous offsets, we find from (47), by 
taking the limiting form of the discrete case, that the coupling coefficients 
for the two polarizations of the TEi n mode arc 



c llmi v± (z) = cwVGO, 



(209a) 
(209b) 



where as usual x(z) and y(z) denote the two rectangular components of 
the displacement of the guide axis from the 2-axis (Fig. 7b) and the 
prime denotes differentiation with respect to z. The differential propaga- 
tion constants are again given by (208). 

Similarly for a continuous diameter variation 



£(•■]*(*) = Cdim^-r'iz), 



(210) 



where r'(z) is the derivative of the guide radius, and we have again 
(208) for the differential propagation constants. 

The geometry for general continuous bends (not confined to a single 
plane) is considerably more complex, as might be anticipated from the 
latter part of Section 2.1.2, and requires special discussion. In general 
the guide axis may be a quite arbitrary 3-dimensional space curve; e.g.. 
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it may be pretzel-shaped. In such extreme cases we require the notation 
of differential geometry to describe the guide axis. We will not attempt 
to discuss the subject of differential geometry itself, 36 but merely use 
some of its simpler results. 

Let us treat the general continuous case as the limiting form of the 
discrete case, as in Section 2.3.2 where the curved axis of the guide was 
confined to a single plane. In the present more general case we use the 
bent cylindrical co-ordinates p,<p,s described in conjunction with (29). 
Note that distance along the (bent) guide axis will be called s in the 
present section (it has previously been called z), so that z can refer to 
distance along a fixed rectangular co-ordinate axis. The other co-ordi- 
nates p and <p arc as defined in Section 2.1.2. We particularly require 
the definition of the <p co-ordinate. Briefly, lines <p = constant are drawn 
on the surface of the guide before it is bent, parallel to the (straight) 
axis of the guide. After bending (perhaps into a pretzel shape) these 
now deformed y>-lincs furnish the ^-co-ordinates. 

Now consider briefly the differential geometry of twisted space 
curves. 36 Let r be the vector from a fixed origin to the curve. Three unit 
vectors characterize the geometric properties of a general space curve: 

t — tangent vector. 
p — principal normal vector. 
b — binormal vector. 
Then the following relations hold true: 

dT (211a) 



(211b) 
(211c) 
(211d) 

(211e) 



ds 


t. 




dt _ 
ds 


1 
P- 

P 




b = 


t X p. 




db _ 
ds 


1 




dv = 


U- 


u 



ds 



In (211) p is the radius of curvature,- the curvature; similarly r is 

P 

called the radius of torsion, - the torsion. The principal normal p lies 
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in the plane of the circle that best approximates the twisted curve in 
the neighborhood of a given point; the radius of this circle is of course 
the radius of curvature p. 

Now let us consider the limiting process in which a guide with a con- 
tinuous three-dimensional bend is considered as the limit of a guide with 
closely spaced discrete tilts. In the following let 0(s) represent the orien- 
tation of the elementary discrete tilt located at distance s along the 
(bent) guide axis, measured in the bent cylindrical co-ordinate system 
described above in the present section and in Section 2.1.2. Then in the 
limit as the discrete tilts become a continuous three-dimensional space 
curve, 0(s) is given by 

^f)--i. (212) 

as t 

From (37) the continuous coupling coefficients to the two polarizations 
of the TEi^ spurious mode are: 

. r±, v f , ± cos 0(s) , , 

C[i m] (s) = C t[m] — — - (213a) 

p(s) 

. u ±,' ) _ r ± sin 6(s) roiqvo 

C[lm (S) — L t [m] ■ 7— r — . (Z\.6D) 

P(«) 

Similarly for the TMu spurious mode: 

,. h-/_n _ ( , + . cos 0(a) , . 

p(s) 

C(,i) (*) = Cku) • — rr^- (214b) 

p(s) is the radius of curvature, determined from (2111)). 0(s) is found 
from (212) as the negative of the integral of the torsion - , determined 

T 

from (211d). 

The results of (213) and (214) are sufficiently general to include such 
things as pretzel-shaped waveguides and other unusual deformations. 
However since we seldom expect such things in practice, we seek to 
simplify these results. We are guided by the simplification in the discrete 
case that is discussed following (37). We assume that the angular devia- 
tion of the guide axis from the z-axis of an x, y, z rectangular co-ordinate 
system remains small ; 

| x'(z) I « 1, 

(215) 
//(*) « 1. 
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Further we assume that the tp = co-ordinate on the surface of the 
guide initially lies in the x-z plane, and that it remains close to this 
plane. Then from (38), by the usual transformation from the discrete, to 
the continuous case, we have 

s ph z, (216a) 

cos6(s) « ,"{z), (216b) 



sinfl(s) 



/(*), ( 21( * 



P(s) 

to be substituted into (213) and (214). The superscripts x and y of 
(213) and (214) in this special case correspond to linear polarizations 
defined with respect to the x and y axes of the fixed rectangular co- 
ordinate system; in general the superscripts of (213) and (214) corre- 
spond to the tp = and tp = x/2 planes in the bent cylindrical co-ordi- 
nate system. The approximations of (210) are found directly for the 
continuous case in Appendix E. 

III. THEORY OF GUIDES WITH RANDOM DISCRETE IMPERFECTIONS 

We now apply the results of Section II to the study of guides with 
random discrete mode converters. The following cases of practical 
interest will be discussed : 

*". Guides made of individual pipes that are perfect right circular cylinders 
of identical diameter an-d length 1 but that hare randomly imperfect joints, 
with either tilts or offsets. The first-order spurious modes in this case will 
be the TEi* family, with the most significant being the forward TE n 
and TE 12 , and for tilts the forward TM,, . Both polarizations of each 
spurious mode must of course be considered. 

ii. Guides made of individual pipes that are perfect right circular cylin- 
ders of identical length l with perfect joints (no tilts or offsets) but of 
slightly different random diameters. The spurious modes will be the TE ,„ 
family, with forward TE (I2 the most important. Each mode now has 
only a single polarization. 

The statistics of the TEoi loss will be determined in terms of the statis- 
tics of the guide imperfections. Only the case of individual pipes of 
identical length will be considered. 

The necessary results from Section II are summarized below in Sec- 
tion 3.1. Section 3.2 states the statistical models adopted for guides 
with random tilts, offsets, and diameter changes. In Section 3.3 the ex- 
pected value of the average TE i loss, and the power spectrum and total 
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power (expected value of the mean square value) of the TE i loss varia- 
tion about its average value, are determined in terms of the statistics of 
the random offsets, tilts, or diameter variations of the guide. In Section 
3.4 these results are extended to cover long guides with mode filters. 
Representative numerical examples are given in Section 3.5 for a 20- 
mile total guide length. Section 3.6 considers the application of certain 
of these results for tilts and offsets to helix guide. 

Most of the work of Section III will be confined to copper guide. As 
before, we assume that the coupling coefficients for ideal metallic guide 
provide a good approximation to those for copper guide, but modify 
the (pure imaginary) propagation constants for ideal metallic guide 
to account for the loss present in copper guide. Of course the analysis 
for diameter changes applies equally well to copper or helix guide. 

3.1 TEoi Loss — Summary of Previous Results 

We give below the TE i normalized loss (in nepers) A, written in a 
form suitable for the purposes of the present section. These relations are 
readily obtained from the results of Section II, and in particular from 
Section 2.4 and (37), (38), (47), (51), (01) to (65), (71d), (71e), 
(161), and (206). First, from (206) the total TE i loss A is given as 
the sum of terms due to the individual first-order spurious modes; 

A&ZA™. (217) 

m 

The term A m due to the ///"' spurious mode is given below in (218) and 

(219). In these and similar relations involving only a single spurious 

mode we omit the subscripts denoting the spurious mode on the coupling 

coefficients and differential propagation constants, in order to avoid 

unnecessarily cluttering up the equation. Therefore for each spurious 

mode 

. .v-i 

-1 « iUo + £ ,1, cos /,-A/J/o (218a) 

2 *=i 

where 

N-k 

A k = e^ al " Z (xhi+k* + ..:>..•,+/), ^ k ^ N - 1. (218b) 



Tilts 



X { * = Ci <Xi COS 0, « CtOi 1 

.c ( x = Ct a, sin 6i w Ci af. 



(218c) 
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Offsets: 

x? = C bi cos di = C b} 
Xi x = CJbi sin di = CJbi\ 



Diameter change : 



.r,-" = C d An 
*,* = o. 



(218d) 



(218e) 



The distinction between A k of (218a) and (218b), the k th Fourier co- 
efficient of the TEoi loss due to a single spurious mode, and A m of (217 ), 
denoting the TEoi loss due to the w th spurious mode, will always be 
clear from the context. In particular cases the subscripts indicating the 
spurious mode will always be enclosed in [ ] or ( ) to denote TE or 
TM spurious modes respectively. 

For the case of zero differential attenuation constant, an alternate 
form is sometimes useful. From (161) for each spurious mode we have 
(with mode subscripts again omitted) 



^ 



£ »' 



-i&filoi 



+5 



£a* e 



i -JW i 



Aa = 0. (219) 



The .r.'s are again given by (218c) to (218e). . 

In the above results \vc have chosen to group the two polarizations of 
each spurious mode (TEi„, and TM„ + ) for tilts and offsets in each term 
of (217). We use the symbols || and 1 to distinguish the two polariza- 
tions, rather than x and \j as in Section II. The length of the individual 
pipes is la , the total length of guide is L N = Nl ■ 

In these results the significant frequency dependence, at least over 
moderate bandwidths, occurs through the A/3, the differential phase con- 
stants between TEoi and the spurious modes. Far from cutoff (the case 
of greatest practical interest) the A/3 are approximately proportional to 
the free-space wavelength X. A great simplification in the subsequent 
analysis is obtained by neglecting the frequency dependence of Aa and 
the .r.'s, setting these quantities equal to their values at the middle of 
the band of interest. Then A m , the contribution to the TE i loss of con- 
version to the m th spurious mode, becomes a Fourier series, periodic in 
X, with random coefficients A k , as given in (218). A, the total TE i loss 
due to mode conversion, is given by a sum of periodic components of 
incommensurable periods, according to (217). 

Over moderate bandwidths, the effect of the frequency dependence of 
Aa is small far from cutoff. For offsets and diameter changes, the fre- 
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quency dependence of the coupling coefficients and hence the .r, is also 
small far from cutoff; however, for tilts the .r, vary approximately in- 
versely with X. Thus, the above approximation may appear questionable 
for tilts, but is used even here because of the resulting simplification in 
the analysis. [Compare the discussion in Section 2.3.7 following (169).] 
We thus regard the loss A as a function of the free-space wavelength 
X, and write 

A = A + 8A (220) 

where A is the average value over free-space wavelength X of the loss 
over some suitable band, and 5.4 is the deviation of the loss from its 
average value. The expected value of the average loss <A>, and the 
power spectrum of 8A and its total power or the expected value of its 
average mean square value <(8A)->, are easily found in terms of the 
statistics of the random Fourier coefficients Ak ; the bar again indicates 
an average over the free-space wavelength X. 

A more exact treatment of the loss statistics, one that avoids the 
above approximations and includes the frequency dependence of all 
quantities, is straightforward but lengthy. A brief discussion, given in 
Appendix F, of the statistics of each of the two terms of (218a) as 
functions of X, verifies the approximate analysis. 

3.2 Statistical ^Fo(l<^l of Guide 

3.2.1 Tilts and Offsets 

We assume that the parallel and perpendicular components of tilt or 
offset (a\ a x or b\ b ± ) are independent Gaussian random variables 
with zero mean and equal variance; tilts or offsets at different joints are 
assumed to be independent, and to have the same distribution. Then the 
magnitude of tilt or offset (a or b) will have a Rayleigh distribution and 
the orientation (8) will be uniformly distributed between and 2x, and 
these two quantities will also be independent. 

From (37), (38), (47), and (218c) and (218d) we have 

C t a l C t a x Cta 

•r" = , x* = • x = . (221) 

CJb 1 Cob* C b 

Thus identical statements to (hose above may be made about x\ x ± , x 
and 8. 

For convenience we state the following results in terms of .r's; any of 
these (Miuations obviously holds true if x is replaced by a or b throughout. 
The rms value of x is denoted by x; thus 
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Ct 2 <a> 
x s = <x*> = , ., , (222) 

C -<ff> 

where the symbols < > denote the expected value of a random variable. 
Then the various probability distributions at any joint may be written 
as follows: 

^^-S^ ^-^-(j)' (223) 
p(x\x x ) = p(x l )p(x x ) 



(224) 



pU)=|exp-(|)", s>0; p{B) = 1, 

< < 2tt 

p(^) = p(x)p(d). 

We have: 

< ;r »> = < ;l . A > = 0, <x'a;- L > = 0. 

<*.««>-* «a- «->-?. (225) 

2 4 

<,-> = V? jj, <.,'-'> = .f 2 , </> = 2x\ (226) 

Subscripts indicating the joint or mode converter have been omitted in 
the above relations. Since different joints are independent, for two differ- 
ent joints i and j we have 

<x.W> = <•'•■"•'•/> = <Xi X x i x > = 0, i * j. (227) 

Subscripts denoting the spurious mode have been omitted in the above 
equations. 

3.2.2 Diameter Changes 

The radius of the i ih pipe r t is given by 

r, = a + e, (228) 

where a is the average radius and «, a small random variation about the 
average. The e, are assumed to be independent Gaussian random vari- 
ables with zero mean and variance e". 

V (e;) = —^exp-^ (229) 
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<e,> = 0; <e, 2 > = t; <e,e,> = 0, i ^ j. (230) 

Then A/ 'i , the change in radius at the i th joint, is given by 

Ar { = a - e<_i, (231) 

and has a Gaussian distribution with zero mean and variance (Ar)" 

p(Ar,) = — L^ exp - i^ ; Ar = V5i (232) 

V27rAr 2(Ar) 

However adjacent joints are no longer independent; from (231) Ar, and 
Afj arc independent only if | j — i \ > 1. Therefore: 

<Ar,> = 0. 

<(Ar,) 2 > = (Ar) 2 = 2e 2 . 

, A A s (Ar) 2 , (233) 

<Ar,Ar, ± i> = — = -t . 

<Ar,A/v> = 0, 1./ - i | > 1. 

Equations (232) and (233) apply also to the corresponding r,, given 
by (218c) for a diameter change, for any spurious mode, by making the 
substitution Ar, — > x, and Ar — > x, where x is given by 

i = CjAr. (234) 

Subscripts denoting the spurious mode have again been omitted in the 
above equations. 

3.3 Statistics of the TE in Loss for a Single Section of Waveguide between 
Mode Filters 

3.3.1 Offsets 

For a copper waveguide section whose only defects are independent 
random offsets at the joints between pipes, the most important spurious 
modes will be the forward TE i2 and TEn , with the other forward TEi m 
modes contributing only a small amount to the total TE„i loss. Neglect- 
ing the frequency dependence of all quantities except the A/3, we have 
from (218a) and (220) for each spurious mode: 

A = A + 5.1, 

A = §.■!„, (235) 

5.-1 = £ A h cos kAplo , 
*:=i 
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where the random Fourier coefficients -4 A . are given in (218b), together 
with (218d). The subscript indicating the spurious mode has again been 

omitted. 

The following moments of At are easily determined from the results of 

Section 3.2.1: 

<Ao> - *' N > (236a) 

<yl„ 2 > = xN(N+ 1). 

<A k > = 0, <An = I (N - l:)c k2 * a '°; I* k g N - 1. (236b) 

<Ajbiii> = 0, ft * I. (236c) 

From (235) and (236a) the expected value of the average TE„i nor- 
malized loss due to each spurious mode is 

<I>=|V (237) 

From (236b) the (discrete) power spectrum of 5.4 is 

l\ = K4, 2 > - £ (N - ft) e mai0 , I ^ k£ N - 1. (238) 

From (236c) the different Fourier components of 5.4 are uncorrelated. 
The total power, or the expected value of the mean square value, of 5.1 
is found by summing over the discrete power spectrum P* . From (238) : 

<WT-> = E P* = f E (N - fc) e* 2Aa " 

i_ ,*"*"■ (239) 

AT - 



i 1 e 2Afl "° 



4 1 - 



■2Aah 



1 - c' V2A "'° 1 
1 - c ,Aa,i J ' 



Strictly speaking, the average over free-space wavelength X indicated by 
the bar in (237) and (239) must be taken either over a single funda- 
mental period of 8A (such that A/3/„ changes by 2tt) or over a band large 
compared to a single period. Let AX and A/ be the interval in wavelength 
and frequency, respectively, corresponding to a fundamental period. From 
(167) and (173) we find 

A/_AX = i* (240) 

/ A k 

where B is the beat wavelength, tabulated in Appendix D, and / is the 
distance between joints in the guide. (This result is similar to that of 
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(178) for the sample point spacing; in (178), however, the length L is 
the total length of guide. ) 

Equation (239) simplifies in two special cases of interest: 

1. Small differential loss over total length L N = Nl ■ 

<JSA) 2 > = % N(N - 1) ; -N2Aah « 1 . (241) 

8 

2. Large differential loss over total length L N , small differential loss 
over pipe length /„ . 

a* N -N2Aah » 1 

<(say> = *^--4t5 (242) 

8 -Aa/o -2Aa7 « 1 . 

We recall again that Aa is negative throughout the present treatment, 
in which the TE i signal mode has lower heat loss than any of the 

spurious modes. 

Referring to Section 3.2.1, <A> and V<(bA) 2 > are both proportional 
to the square of the rms offset at the joints between pipes. < A> is pr o- 
portional to L N , the length of the waveguide section. V<(5A) 2 > is 
initially proportional to L N , when L N is small enough so that the differ- 
ential loss may be neglected; for large L N it becomes proportional to 

\/T N . 

The power spectrum P k of (238) has its maximum value for k = 1, 
and decreases monotonically as k increases to N — 1. For small differ- 
ential loss P k is triangular; for large differential loss it is exponential. 
The "3-db bandwidth" of the power spectrum P k , the value of k for 
which P k is equal to half its maximum value, is related to the rate of 
variation of the TE ()1 loss component due to a particular spurious mode. 
We have: 

k 3db = ^; -N2Aak « 1 (243a) 

*** - -^T ! -^Aa/n » 1 . (243b) 

Making use of (167) and ( 173), we find 

Afadb _. AX; i( /6 B 



J A hSdbio 

for the interval in free-space wavelength or frequency (of the radiation 
in the guide) corresponding to the 3-db bandwidth of the TE i loss 
variation. 



1108 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1962 

The statistics of the total TE 01 loss due to mode conversion are given 
simply by summing over the contributions of the individual spurious 
modes. From (217) we have 

<A> = £ <A [m] >, (245) 

[m] 

<j8Ay-> = 2 <(5A [m] ) 2 >, t (246) 

where [m] indexes the TE ]ro + spurious modes. The individual terms in 
(245) and (246) for each spurious mode are given by (237) and (239) 
(or (241) and (242) in special cases), in which the subscript [m] has 
been omitted for convenience. The power spectra of the 5-4 |m] are given 
by (238). From (222) and (237), 

<A> = < -^-ZC oM \ (247) 

2, [m] 

A practical waveguide system will contain mode filters for the TEi,„ 
modes at a close enough spacing so that the differential losses for the 
important spurious modes are small in each section. For this special 
case we have from (241) 

<(5Z7 2 > = <b>N W ~ 1} £ C o[m] 4 ; -N2Aa [m] h « 1 . (248) 

» [m] 

In (247) and (248) <6 2 > is the mean square magnitude of offset, C [ m \ 
is Morgan's coupling coefficient between TE i and TE im for offsets, 
given in Appendix A, and Aa[,»] = a[ou — «[i.«] , the difference in at- 
tenuation constants of the TE ot and TE im + modes. 

Formulas and numerical values for the various coupling coefficients 
and beat wavelengths are given in Appendix A and D. For a frequency 
of 55 kmc and a 1-inch guide radius, (247) and (248) become, summing 
over the nine propagating TE Jm modes: 

<A> = ^^ (1.107 + 4.581 + 0.641 + 0.271 + • • •) 

= <b 2 > AT(3.519) ; <6 2 >ininch 2 . (249) 

<(5Ay 2> = < b2 >" N ^ ~ H (1.226 + 20.984 

o 

+ 0.411 +0.074 + •••) 



t No cross terms appear in the summation of (246) because, subject to the 
approximations of the present section that neglect the frequency variation of the 
coupling coefficients (C t , C„ , and Cd) and the Aa's, the Fourier components of the 
different &A\ m \ have incommensurable periods, and hence their total powers or 
mean square values may be simply added. The cross terms are treated exactly in' 
Appendix F. 
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= <bVN(N - 1) (2.842) 



V<(5/lp = <tf>VN{N - 1)(1.686); <6 2 > in inch 2 , 

( 2o0 ) 
-N2A<x lm] l « 1 . 

The most important terms are those due to TE U + and TEi 2 + . 

3.3.2 Tilts 

We next consider a copper guide whose only defects are random tilts 
at the joints between pipes. The spurious modes are the forward TE Jm 
with TEi2 + and TEn + the most important, as in the offset case above, 
and in addition the forward TM„. The effects of the TE lm + modes on 
the TEoi transmission are given by the results of Section 3.3.1 above 
[see (235) to (246)], using of course the appropriate coupling coefficients 
for tilts [see (218c)]- However, TM U + requires special consideration. 

Equations (235) and (236) apply to TM U + as well as to the TE, m 
modes. For TM n + , A/3 = Aa in copper guide. Thus the beat wavelength 
for TMn + is very long — 3195 feet in 2-inch I.D. pipe at 55 kmc — com- 
pared to the beat wavelengths of the TEi„, modes, and long even com- 
pared to the length of guide sections between mode filters in a practical 
waveguide system. t Thus, the bandwidths we will consider (e.g., 50-60 
kmc) are only a small portion of the fundamental period of A iU) + , the 
TMn + component of the TE 01 loss, as given in (240), and so the sum- 
mation of powers of Fourier components given in (239) is no longer 
valid in determining the mean square loss variation. In fact A ( i» + will 
be almost independent of frequency (except for the slow variation of 
coupling coefficient, which is inversely proportional to free-space wave- 
length X, neglected over moderate bandwidths in the present analysis). 

Thus, consider a section of copper guide short compared with the 
TM U + beat wavelength (3195 feet ). Then both Ac* and A/3 may be set 
equal to zero. From (219) the TM U + component of the TEoi loss is 

^ + -\[pf] + l[t^]. (25D 

This is simply equal to the TEoi loss due to TMu + at a single tilt equal 
to the net tilt between the input and output ends of the guide. This 
result is obvious from the fact that we have neglected both attenuation 
and phase shift in the relatively short sections of guide under considera- 
tion. The loss is therefore independent of the lengths of guide between 

t This is true only for unmodified copper guide. A thin dielectric coating will 
reduce the TEoi-TMu + heat wavelength to much smaller values; 121314 the present 
treatment of TMn + will obviously not apply in such a case. 
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discrete mode converters, and these lengths may thus be set equal to 
zero, yielding only a single tilt whose magnitude and orientation are 
equal to the net tilt between the ends of the guide. This result holds 
true for TM U + for any arbitrary continuous bend of the guide axis, 
which may be considered as the limit of a series of discrete tilts (Section 
2.3). 

To this approximation -4(in + is thus independent of frequency (ex- 
cept for the slow variation of coupling coef ficient with A). TM M + thus 
contributes to the average loss A but not to (8A) 2 . Thus (245) contains 
an extra term, while (246) remains unaltered. 

<A> = <I> TE + <A m) + > (252) 

<I> TE = 2 <^[«.]> ( 253 ) 

[m] 

<(5A)^> = Z <(5A [m] ) 2 > ( 254 ) 

[m] 

The index [m] again indexes the TEi m + modes. <i4 [m] > and <(SA [mi ) 2 > 
are again given by (237) and (239) (or (241) and (242) in special 
cases) and <A (U ) + > is also given by (237), with the appropriate coupling 
coefficients for tilts. Equations (238), (243), and (244) remain true 
for the TEi,„ + components of the TE 01 loss. The results analogous to 
(247) to (250) for tilts are summarized below. 

7 <OC~>N V~» it '!■ /<srr\ 

<A> TE = -—— Z. Ci W (-255) 

2 [m] 

<-l(..) + > = ~- C IW 2 (TM { „> + ) (256) 

For small differential loss (because of mode filters), 

<(8Ay> = <a ' yN( X ~ l) E C', [w| 4 ; -N2Aa lmi l {) « 1 . (257) 

<a> is the mean square magnitude of tilt. Substituting numerical values 
for the C t + 's from Appendix A, for a frequency of 55 kmc and a 1-inch 
guide radius, (255)-(257) become, summing over the nine propagating 
TEj m + modes: 

<I> TE = < ^L (29.465 + 81.085 -f 0.616 + 0.057 + • • •) 

= <a>N( 55.619); <a> in radian 2 . (258) 

Gl ai) + > = < a 2 >N( 14.598); <a> in radian 2 . (259) 
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<(jj)2> = < a>N( _ N - 1) ( 868 .i 8 + 0574.74 4- 0.38 + ■ ■ •) 

o 

- <a>*N(N - 1)(930.41) 

V<(SAY 2 > - <a2> VN(N - 1) (30.503); <a> in radian 2 , 

(200) 
-N2Aa [m] l « 1. 

Again the most important terms are those due to TE n + and TEi 2 + , and 
for the average loss, to TMu + in addition. 

For a given (copper) guide A (U ) + , the TM U + contribution to the 
(average) TE M loss is determined simply by the net tilt between the 
input and output of the guide, as discussed above. The present model, 
which assumes that the only imperfections are tilts at the joints between 
perfect pipes, seems grossly unrealistic as far as the effects of TMn are 
concerned for any practical guide, for at least two reasons. First, practical 
pipes will have long bows or gradual curvature of the guide axis; this 
factor will probably be much more important in determining the net 
tilt between the guide input and output than tilts of the very small 
angles of interest here. Second, practical guides may be subject to 
mechanical constraints of various types which will also introduce slow 
variations in curvature of the guide axis. 

In contrast, gradual curvature of the guide axis will have little effect 
on the TEim components of the TE,,i loss for reasons indicated in Section 
2.3.7 and to be discussed in detail in Section IV; this is so because only 
straight ness deviations whose wavelengths are approximately equal to 
the beat wavelengths of the important spurious modes contribute to 
the TEoi l (,ss hi copper guide. 

Consequently the effect of TMn* on the average loss has been stated 
separately for the particular model discussed here; it is given in (250) 
and (259) for whatever tutorial value it may have. As stated above, 
TMn + will have no significant effect on the variation of TEoi loss about 
its average value for the relatively short mode filter spaeings which 
must be used in a practical waveguide system. 

:<.:?.:* Diameter Changes 

Finally, consider a copper or helix guide whose only defects are 
random diameters of the individual pipes, which are perfect right circular 
cylinders and have no tilts or offsets at their joints. Again from (218a) 
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and (220) for each spurious (TE 0m ) mode we have, as in (235) for 
offsets or tilts: 



A = A + BA 

A = \Ao 

N-l 



(261) 



8 A = 2 -4/l-cos &A/3/ . 



*=i 



The random Fourier coefficients are again given by (218b), together 
with (218e). 

ff-k 

A k = e kAal " E.w +fc - (262) 

i=l 

The subscript indicating the spurious mode is again omitted. 

The moments of the A k must be slightly modified from those given 
in (236) for independent offsets and tilts, because of the correlation 
between adjacent diameter changes imposed by the present mathemati- 
cal model and because each spurious mode now has only one polarization 
rather than two. From (228) to (234): 



<ylo> = &N, <A„ 2 > = x(N 2 + 3N - 1). 



cl,>= -| (N - l)e* a \ 



<Ai> = I (N 2 + 5N - 8)e 2Aa '° 



<A k > = 0, <A, 2 > = x - [S(N - k) - l] e fc '-' Aa '° ; 

2 ^ k ^ N - i . 



(263a) 
(263b) 

(263c) 



<A*A» -*******'< 



-(N-l) 

+ UN - I) 





I - 


- k 


= 


ll 


I - 


- k 


= 


2 


I - 


- k 


> 


3_ 



,44. - 1 iSnl. 

I 



- h(N - l)(N + 4) 
+ *(iV - 2) 

n 



>: 



I = I 
I = 2 
I > 3 



, < /c < J, (263d) 



(263e) 
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From (236a) and (261) the expected value of the average loss due 
to each spurious mode is again 

Ol> = |\v, (264) 

as in (237) for offsets or tilts. However, the statistics of the Ak's, as 
given in (263), differ from those of (236) for tilts and offsets. Because 
<.!,> is no longer equal to zero, the expected value of the TE i loss will 
have a fundamental periodic component, in addition to a dc component. 
Consequently it is convenient to rewrite the first relation of (261) as 
follows: 

A = A + C-ii> cos A/3/o + 5.4'. (265) 

Thus, 

JV-l 

5/1' = (ill - C4i>) cos A/3/o + £ A k cos kA0l Q . (266) 

fc=2 

<Ai> is given in (263b). A is the average loss as before, <.4 X > cos A/3/ 
is a slowly varying sinusoidal component of loss whose period equals 
the fundamental period of the TE i loss (see (240)), and 5.4' includes 
the remaining random loss variations. The (discrete) power spectrum 
/Y of 8 A' is given by 

h(<A*> - Cl,> 2 ); k = 1. 
/V= 2 (267) 

i<A/> ; 2 ^k ^ N - 1. 

From (263d) the different Fourier components of 5/1' are now uncorre- 
lated only if their indices differ by three or more. This does not affect 
the calculation of the total power of 5.4', which remains simply the 



£-1 .4 (■ 2A„/ r , _ -V2Aa/ -| 



e 2*al 



I _ c 2Aa/ 



_ e (A'-l)2Aaf l 



i _ „y»— iizaaie I , K-N U e 2AaJ ( 



(268) 



Again for small or large differential loss, (268) simplifies: 
1 . Small differential loss over total length L s = Nl Q . 



..4 

x 



<(5.lT> = i- (3A r - 1 )(N - 1); -iV2Aa/„ « 1 . (269) 

o 
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2. Large differential loss over total length L N , small differential 
loss over pipe length / . 

4* W -N2Aal » 1 

<(*A') 2 > = | ^y; • (270) 

8, -Ao/o -2Aq/ «1 



<A> and <(8A') 2 > have the same general functional form as for tilts 
and offsets. The "3-db bandwidth" of P*' is approximately the same 
as given in (243)-(244) for tilts and offsets. 

The statistics of the total TEui loss due to mode conversion are now 
given by summing over the TE 0m + spurious modes. 

<A> = £ </[ w > (271) 

<{8Ay> = Z<(^W7 2 >. (272) 

The individual terms in (271) and (272) are given by (264) and (268) 
(or (269) and (270) in special cases). In addition, each spurious mode 
will contribute a single sinusoidal component to the TE i loss, of magni- 
tude given by the middle term of (265). From Section 3.2.2 and (264), 

<A> = ^yNj^Cunf (273) 

where <e 2 > is the mean square variation of pipe radius. 

In discussing the effects of the TEi„, + modes in the first two parts of 
this section, it was assumed that the line contained ideal mode filters 
at a close enough spacing so that the differential loss in each section 
could be neglected. However, a practical mode filter in 2-inch guide 
presently consists of a section of helix guide, which has a low loss for 
the TEii„, spurious modes (although it effectively suppresses all other 
spurious modes). In the present case it is therefore assumed that there 
will be no mode filters in the entire length of line between repeaters. 
Thus, the total differential loss will be large, and the approximate result 
of (270) yields 



<mt> _SjgNx^-, - N2 ^ h>>l - (274) 
2/ [m] -Aa [m] -2AaZ « 1- 

Substituting numerical values for the coupling coefficients and the 
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Aa's, for a frequency of 55 kmc and a 1-inch guide radius, (273) and 
(274) become: 

<I> = <e a > N (2.424 4- 0.771 + 0.394 4- 0.244 + • • •) 

= <e 2 > N(4.414) ; <e 2 > in inch 2 . (275) 

<(8A 7 y> = - <e J N X 10 4 (7.197 4- 0.274 4- 0.037 4- 0.008 + • • •) 

= < l2JL(u.2S x io 4 ) 

1(1 

, 'at <e"> in inch' 

V<(SA'T-> = <e' J > 4/ P- (335.9) ; 

V '° /» in feet . 

-N2Aa [m] lu » 1 , -2Aa [ml /„ « 1 . (276) 

In these results the summation has been extended over all of the propa- 
gating TEom 4 spurious modes, TE O i + -TE 09 + . The most important modes 
arc the first few TE ,„ + ; for <(5.1') 2 > only TE 02 + and TE I13 + are significant. 
Finally, the sinusoidal component of the TE 01 loss contributed by the 
TEo„, ' spurious mode is, from (265) and (263b), 

<.lii,„]>cos A/3(„,|/o = -<( 2 >(N — l)C rf[m ] 2 cos A0[ m] lo ; 

— A«[ „,]/() « 1. 

3.4 TEqi Loss Statistics of a Long Guide with Ideal Mode Filters 

Consider a long guide made up of M sections of imperfect guide of 
equal length and the same statistical parameters, separated by ideal 
mode filters. We must evaluate the over-all transmission statistics of 
the guide in terms of the transmission statistics of each section, given 
in Section 3.3. The transmission parameters of such a guide with ideal 
mode filters are given in (146). Since for the present we are concerned 
with only the over-all TE i loss A, we have 

A = Z k A. (278) 

k A is the total TE i loss, due to all spurious modes, of the k lh section of 
guide. From (220) we write 

k A = k A + 5,-1. (279) 

Further, for the loss of the entire guide (with mode filters) 
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A =1 + 8 A (280) 

where 

A = D*I, SA = EM- (281) 

We assume that each section of guide between mode filters is statisti- 
cally independent of all other sections. Then from (281) 

M 

<A> = £ < k A> = M< k A>, (282) 

M 

<T&Ay-> = Z <(maT 2 > = ^<(Ml 2 >, (283) 

where M is the number of sections of guide separated by ideal mode 
filters. 

Finally, from (278) and the central limit theorem, if M is large the 
over-all loss A , regarded as a function of free-space wavelength X, will 
be a Gaussian random process; this random process in general will not 
be stationary, although over the relatively narrow bands of interest it 
may often be assumed stationary. 

3.5 Numerical Examples 

In the present section several numerical examples arc presented to 
provide concrete illustrations of the above results. A 20-mile total guide 
length, made of 2-inch I.D. pipes 10 feet long, with equally spaced ideal 
mode filters, is considered in all cases; the operating frequency is taken 
to be 55 kmc. It is assumed that the mode filters have infinite loss for 
the TE,„, + and TM n + spurious modes, zero loss for the TEoi signal 
mode and the TE 0m + spurious modes. The results for tilts and offsets 
apply to copper guide; the results for diameter variations apply equally 
well to either copper or helix. 

3.5.1 Offsets 

Assume an rms offset such that the additional average loss due to 
mode conversion to the forward TE lm modes is 1 db/mile (compared 
to the theoretical TE i heat loss at 55 kmc of 1.54 db/mile). Two eases 
are considered (see Table II): (1) mode filter spacing such that the 
rms loss variation for the 20-mile line is 1 db, and (2) a 200-foot mode 
filter spacing. The formulas for zero differential loss are used for sim- 
plicity, since the differential loss for TE ]2 + , the most important spurious 
mode, will remain small. 
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Table II 



Average loss 

RMS total loss fluctuation for 20-mile line 

RMS offset 

Mode filter spacing 

A/ 3d b for TE.2+ 



Case 1 



1 db/mile 

1 db 

7.87 mils 

1160 feet 

209 mc 



Case 2 



1 db/mile 
0.407 db 
7.87 mils 
200 feet 
1211 mc 



3.5.2 Tilts 

Assume an rms tilt such that the additional average loss due to mode 
conversion to the forward TEi,„ modes is again 1 db/mile. The addi- 
tional average loss due to conversion to TM n + is stated separately 
because the present model is unrealistic as far as TM n + is concerned, 
for reasons stated in Section 3.3.2. Two cases are again considered (see 
Table III): (1) mode filter spacing such that the rms loss variation for 
the 20-mile line is 1 db, and (2) a 200-foot mode filter spacing. The 
formulas for zero differential loss are again used for simplicity. 



Table 


III 






Case 1 


Case 2 


Average loss TEi„, + modes 

Average loss TMn + 

RMS total loss fluctuation for 20-mile line 

RMS tilt 

Equivalent crack on one side of joint 

Mode filter spacing 

A/adbforTEi 2 + 


1 db/mile 
0.262 db/mile 
1 db 
0.114° 
3.96 mils 
890 feet 
272 mc 


1 db/mile 
0.262 db/mile 
0.465 db 
0.114° 
3.96 mils 
200 feet 
1211 mc 



The TEoi average loss due to TM n + conversion will depend only on 
the net angle between the input and output of a waveguide section 
between mode filters for a mode filter spacing short compared to 3195 
feet, the TMn + beat wavelength in copper guide. As discussed in Sec- 
tion 3.3.2, this angle (Table IV) will depend principally on long bows 
in the pipes and on the way in which the guide is laid. 





Table IV 








Mode filter spacing 


200 feet 


1000 feet 


Average loss, TM U + 


1 db/ 
mile 


0.1 db/ 
mile 


ldb/ 
mile 


0.1 db/ 
mile 


RMS net angle between input 


and out- 


1.00° 


0.316° 


2.23° 


0.706° 
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3.5.3 Diameter Changes 

For diameter variations there are no mode filters for the spurious 
TE 0m + modes. Therefore, the mechanical tolerance required to yield an 
rms loss fluctuation <(5A') 2 > of 1 db for the 20-mile line is determined 
(see Table V) ; the additional average loss will now be very small. 

Table V 



A verage loss 

<(8A') 2 > : RMS total loss fluctuation for 20-mile line 

RMS diameter variation 

A/ 3l ib for TE 2 + 



0.214 db/mile 

1 db 

6.50 mils 

11.3 mc 



Sinusoidal Components of TE ol Loss for 20-Mile Line 



Mode 


Peak-to-Peak Amplitude 


Period 


TE 2 + 
TE 03 + 
TE 4 + 
TE« B + 


4.690 db 
1.492 db 
0.762 db 
0.472 db 


4791 mc 

1832 mc 

974 mc 

600 mc 



3.6 Helix Guide 

While the above results for diameter changes apply to both copper 
and helix guide, the results for tilts and offsets apply only to copper 
guide. Equivalent results for helix would require the coupling coeffi- 
cients for the normal modes of the helix at tilts and offsets. However, a 
very simple argument shows that A, the average loss, will be identical 
in a helix and a copper guide which have identical tilts or offsets; A (n) , 
the TMn + component of the TE i loss, must now be included hi the 
copper pipe average loss in the case of tilts, as shown below. In addition, 
the spurious modes have such a high loss in helix that the TE i loss 
fluctuations will be very small. 

From (235), (218) and Section 2.2 we see that A is simply the sum 
of the TEoi signal losses at each individual discrete mode converter 
(tilt or offset), where by signal loss we mean —In s o , where s 00 is the 
TEoi transfer coefficient of the discrete tilt or offset. From Section 2.1.4 
Son is identical in copper and helix guide with equal tilts or offsets. 

Therefore the above results for the expected value of the average 
TEoi loss for tilts and offsets hold equally well for helix waveguide. The 
TEoi loss fluctuations in helix will be very small for these cases. 
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3.7 Conclusions 



Experimental copper waveguides have been built whose tolerances 
are far better than those of the numerical examples in Section 3.5. 
Since the average loss and the rms loss variation are proportional to the 
square of the rms tolerance, it is clear that tilts and offsets at joints 
and uniform diameter variations of the individual pipes will not con- 
tribute significantly to the observed TE ] loss in these waveguides. 
Consequently the additional TEoi loss observed in present experimental 
waveguides must be due principally to continuous mode conversion, 
and in particular to continuous random deviations from straightness 
of the individual copper pipes themselves. 8 The continuous case will be 
treated in Section IV. 

The added TEoi average loss due to TM U + conversion in copper wave- 
guide (unmodified by a dielectric lining or anything else) is a function 
only of the net angle between the input and output of each waveguide 
section between mode filters, for a reasonable mode filter spacing. The 
tolerance on this angle must be held to a few tenths of a degree to keep 
this loss component down to 0.1 db/mile. 

The present analysis has been restricted to equally spaced mode 
converters, i.e., individual pipes of the same length. If the pipe lengths 
are allowed to become random, instead of starting from (71d) we must 
start with (71a) and (71b). The TE„i loss due to a single spurious mode 
will still have a discrete power spectrum, but the discrete components 
will no longer be equally spaced, and consequently the TE i loss will 
no longer be periodic. The frequencies as well as the amplitudes of the 
discrete components must now be treated as random variables. Aside 
from these minor differences, the analysis should be similar and lead 
to similar results. 

We refer again to the treatment of Appendix F, where more exact 
expressions for the TEoi loss statistics in the discrete case are derived 
without neglecting the frequency (or X) dependence of the coupling 
coefficients and the differential attenuation constants, as in the above 
treatment. It is found that these approximations are valid for our 
present purposes. 

Finally we note that by means of the Kronecker product, it is possible 
to compute certain of the TEoi transmission statistics exactly — i.e., 
without using perturbation theory — for the case of statistically inde- 
pendent discrete mode converters. 8 This treatment requires that the 
individual conversion coefficients be known exactly; unfortunately only 
in the idealized two-mode case is the exact form of the coupling coeffi- 
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cients known. Such calculations may be used to check certain of the 
above approximate results in the two-mode case. 

IV. THEORY OF GUIDES WITH RANDOM CONTINUOUS IMPERFECTIONS 

This section applies the results of Section II to the study of multimode 
waveguides with random continuous mode conversion. Continuous 
mode conversion arises from gradual continuous changes in the geometric 
properties of the guide, such as curvature of the guide axis, variation 
of the guide diameter, or changes in the cross section of the guide such 
as ellipticity, etc., as opposed to the discrete case, studied in Section III. 
The statistics of the TEoi loss-frequency curve are determined in terms 
of the statistics of the different guide imperfections. In particular, the 
average TEoi loss and the rms value and the power spectrum of the 
TEoi loss variations are calculated. 

The most important practical application of these results to date has 
been to study the effects of random straightness deviations of the guide 
axis. Here we consider only small unintentional straightness deviations, 
either arising in the manufacturing process of the individual pipes 
themselves or resulting from the way in which the guide is laid. We 
exclude from consideration the case of large intentional bends (to go 
around corners), which couple TE 01 to the degenerate forward TMu 
mode. The spurious modes of interest here are thus the TE lm family, 
with the forward TE n and TE X2 the most important. The present analysis 
indicates that very small random straightness deviations in a certain 
spectral region (i.e., having mechanical wavelengths lying in a certain 
range), having an rms value of a fraction of a mil, are primarily responsi- 
ble for the observed departure of the TEoi transmission from its theoreti- 
cal value in present copper guide, 8,39 causing an increased average loss 
and random fluctuations about this average. In addition, the analysis 
indicates that random straightness deviations will be equally important 
in helix or dielectric coated waveguide in increasing the average TEoi 
loss; however, the high spurious mode loss in helix will effectively re- 
move the TEoi loss fluctuations. ' 

While random straightness deviations are the most important manu- 
facturing tolerance at present, the same methods are easily applied to 
study the effects of other tolerances of the guide. The present section 
will also consider random diameter changes, which produce the TE „, 
modes, random ellipticity, which gives rise to the TE 2m modes, and 
higher-order deformations of the cross section, which produce TE„ m 
modes of higher angular index. 

In order to specify the statistics of the guide, we assume that each 
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type of imperfection (e.g., deviation of the guide axis from straightness, 
diameter variation, ellipticity, etc.), regarded as a function of distance 
z along the guide axis, is a stationary Gaussian random process of known 
spectrum. The various continuous coupling coefficients are of course 
proportional to the geometric imperfections, and thus are also Gaussian 
random processes. 

The analysis of the continuous case is greatly simplified if the differ- 
ential loss between the TE nJ signal mode and the various spurious modes 
may be neglected over the lengths of interest. This approximation will 
be made throughout the present section. As discussed above, a practical 
system using copper guide will contain regularly spaced mode filters 
that have a high loss for all spurious modes except the TE 0m family. 
The mode filter spacing will be sufficiently small so that the differential 
losses may be neglected for all important spurious modes except the 
TE 0m family. 

For the TE 0m spurious modes the effective line length will be the total 
distance between repeaters; obviously the total differential loss is no 
longer negligible. However, the results for zero differential loss will be 
stated for this case to get at least a rough upper bound on the importance 
of diameter variations, for both copper and helix guide. 

4.1 TEm Loss — Summary of Previous Results 

We give in the present section the normalized TE i loss (in nepers) 
A , for the case in which the total differential loss (for each section of 
guide between mode filters) may be neglected, so that we may set Aa = 0. 
From (206) the total normalized TE 0] loss A is given by a sum of terms 
due to the individual first-order spurious modes; 

A &T,A m . (284) 

The .!„, are given by (285) to (288) below. These results are obtained 
from Section 2.3.7. We again omit subscripts denoting the spurious 
mode where no confusion will arise. 

A « \ | / | 2 , (285a) 

sin 7T I -3— — n 1 

The e n are the Fourier coefficients of the continuous coupling coefficient 
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c(z), defined by 



c(z) = E c n e i2 ' n " L ; c_ 7l = c n * . (286) 

n=— °° 



At the n th sample point, defined by 



^=7i or B = — = -, (28/) 

2tt A/? n 



we have 



/(n) = Lc„(-1)", (288a) 

A(n) tt\ \c n \ 2 . (288b) 

The coupling coefficient c(z) is given in terms of the various geometric 
parameters in .Section 2.3. 

As in the discrete case (Section 3.1), the principal frequency de- 
pendence in these results occurs through the A/3's, which far from cutoff 
are approximately proportional to the free-space wavelength X. Over 
the moderate fractional bandwidths of interest, any frequency depend- 
ence of the coupling coefficients will be slow and may be neglected. 
From the discussion of Section 2.3.7, any frequency dependence of the 
coupling coefficient may be taken into account in (285) by calculating 
c„ at the frequency corresponding to the n th sample point, rather than 
at the operating frequency. 

We regard the loss A as a function of the free-space wavelength X, 
and write A as follows; 

.4 = <A> + 8A. (280) 

We determine the expected value of the loss <A>, and the power spec- 
trum of 8A and its total power or mean square value <(8A) >, in terms 
of the power spectrum of the random coupling coefficient c(z) (and 
consequently of the random geometric imperfection of the guide). 



4.2 Statistics of Fourier Coefficients of c(z) 



37,40 



We assume that the geometric imperfection of the guide (e.g., devia- 
tion of the guide axis from straightness) is a stationary Gaussian random 
process with a known power spectrum. The continuous coupling coeffi- 
cient c(z) to the particular polarization of one of the spurious modes 
will be a similar random process, since the coupling coefficients are 
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simply proportional to the corresponding geometric imperfection. There- 
fore c(z) is a stationary Gaussian random process with a power spectrum 
<S(f). Thus, if R(t) is the covariance of c(z), 

R(r) = <c(z)c(z + r)>, (290) 

then 

S(f) = f «(t) e~ iUiT dr. (291) 

J— CO 

Consider the Fourier series expansion of c(z) over the interval 
< z < L, given in (165) and (280). 

00 

c(z) = £ e 1 « fl,,,,L i c_„ = c n *. 

„=-» (292) 

c„ = a„ 4- jo„ — | c„ I e 

The c„'s will be complex Gaussian random variables; i.e., a„ and b„ , the 
real and imaginary parts of c n , will be Gaussian random variables with 
zero mean. If L is sufficiently long, the a„'s and 6„'s become almost in- 
dependent, and hence the c„'s become almost independent complex 
Gaussian random variables. Thus, the magnitude and phase of each c„ 
are independent and have a Rayleigh and a uniform distribution re- 
spectively. The mean square value of the n Fourier coefficient <| c„ \ > 
is then simply related to the power spectrum £(f) of c(z). We have 
approximately for large L: 

<!*!•> -«.■-£■*©. (29.) 

<c n c„*> = 0, n ^ m. (294) 

The quantity c n denned in (293) is the rms magnitude of the n tb Fourier 
coefficient. The various probability distributions for the real and imagi- 
nary parts or for the magnitude and phase of the Fourier coefficients may 
be written approximately as follows: 



(295a) 



p(o„) = -j=rT- exp -(/) ; 

VttC,, \C„/ 

1 // \ 2 

p(b n ) = 77=7- exp -(7^) • 

y/TTCn \Cn/ 

p(a n ,b n ) = p(a n )p(b„). (295b) 
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/. in 2 | C n | (\ C„ 

p( | c„ | ) = ■ ' exp -1^- 

V Cn ' (296a) 

p(<Pn) = ^ , < *> < 27T. 

v(\c n \,<Pn) =p(\c n \)p(<p n ). (296b) 

Since the different Fourier coefficients are approximately independent, 
p(c n , c m ) = p(c n )p(c m ), | n | 5* | m |. (297) 

Finally, the first few moments of the | c„ | are of interest. 

<| c „ |> = y^L c„ ; <| c „ | 2 > = c„ 2 ; <|c„ | 4 > = 2c„ 4 . (298) 

The results of this section provide a good approximation for the prac- 
tical cases of interest in which L, the line length, is of the order of a few 
hundred feet and the power spectrum S(£) varies slowly in the range of 
interest, which includes mechanical wavelengths from a few feet to a 
few inches, depending on the spurious mode. These results become exact 
if the coupling coefficient c(z) has a white power spectrum (£(f) = 
constant, or equivalently R{t) « 5 (r), the unit impulse). 

4.3 TEoi Loss Statistics for a Single Section of Waveguide between Mode 
Filters 

4.3.1 Single Spurious Mode, Single Polarization 

From the relations of Section 4.2, H-^—j , given in (285b), will be a 

complex band-limited Gaussian random process; the real and imaginary 
parts of / are independent Gaussian random processes with flat power 
spectra over the range | v \ < i f33,34 Since by (285a) A is proportional 
to the square of the magnitude of / and is thus proportional to the sum 
of the squares of the real and imaginary parts of 7, the power spectrum 
of A may be determined from the well-known analysis for the response 
of a square law device to Gaussian noise. 37,40 The square of a Gaussian 
random process has, in addition to a dc component, a random component 
whose power spectrum is twice the convolution of the input power spec- 
trum with itself. Since the real and imaginary parts of I have flat 
band-limited power spectra over the range | v \ < i, the random com- 

t v again indicates the independent variable of the Fourier transform of /, or 
some other quantity of interest, with respect to the normalized independent vari- 
able -£- . See the footnote on p. 1084. 
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ponent of A ( -£— ) will therefore have a triangular power spectrum over 

the range | v \ < 1 . 

We first separate the TE i loss .4 as before [see (289)]. 

A = <A> + SA, (299) 

where as usual we omit subscripts denoting the spurious mode. Then 
from the results of Section 4.1 and 4.2, the expected value of the TE i 
loss is given by 

<A(K)> = f <| c n | 2 > = ~s(^j , (300) 

where from (174) 

K-%1 (301) 

is the free-space wavelength corresponding to the n th sample point and 
D is the constant relating the differential propagation constant A/3 to 
the free-space wavelength X [see (173)]. Substituting (287) and (301) 
into (300), and using the result to interpolate between the sample points, 

</1(x)> = < ^sg r x)^. S (L) > (302) 

relating the expected value of the TE i loss due to a single spurious mode 
(single polarization) to the power spectrum of the coupling coefficient 
between TE i and the spurious mode. B is the beat wavelength between 
TEoi and the spurious mode. 

The (continuous) power spectrum of 5.4 in the region close to X is 
given by 

«.>-£*£x)u- M > 

x (303) 

= fs-©(i-|H), I H < i- 

In deriving (303) we have tacitly assumed that S(£), the power spec- 
trum of the coupling coefficient c(z), varies only slowly in the region of 
interest, so that / and A are approximately stationary over moderate 
bandwidths; however, this is not a serious restriction. Equation (303) 
may be obtained either from the known results on the square of a 
Gaussian noise' 1 '' 10 or directly from .Sections 4.1 and 4.2. The total power, 
or the expected value of the mean square value, of 8A may be found by 
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integrating the power spectrum P(v). From (303) 
,- L 2 a2 /D A L 2 „ 2 /l 



<(8A)> = T^V = T s \i)> (304a) 

V<J8AW> = ^s(^ x) = %*(£) = <A>, (304b) 

for a single polarization of a single spurious mode. Alternately, the re- 
sults of (304) may be obtained directly (at the sample points) from 
(288b), (293) and (298); it is apparent that (302) and (304) hold for 
quite general power spectra *S(f). 

The power spectrum P{v) of 8 A is triangular; from (303) the 3-db 
bandwidth is 

vuh = i (305) 

The interval in free-space wavelength or frequency corresponding to 
the 3-db bandwidth of SA is thus 

A/3JI) ^ A\;i db _ 2B , ,-. _ _ v 

/ XL 

and is thus simply twice the sample point spacing (for /), given in 
(178). [Compare (306) with (244) and (243a).] 

Finally, we consider the probability distribution for A, considering 
for the present only a single polarization of a single spurious mode. As 
discussed at the beginning of the present section, A is the sum of the 
squares of two independent Gaussian random variables. Alternately, A 
may be regarded as the square of a Rayleigh-distributed random vari- 
able. Consequently, for a single polarization of a single spurious mode, 
A has an exponential probability density. 

pU ) = _Lexp-A, (307) 

where the average loss <A> is given in (302). We recall that this result 
(and all others of the present section) is based on the assumption of zero 
differential loss, Aa = 0. Equation (307) holds equally well for the cor- 
responding discrete case of Section III. 

4.3.2 Single Spurious Mode, Two Polarizations 

The above results are easily extended to two polarizations of the 
spurious mode. We assume that the two orthogonal components of the 
geometric imperfection giving rise to mode conversion are independent 
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Gaussian random processes with the same statistics. For example, con- 
sider random straightness deviations, which couple TE ( „ to the TE lm 
family but principally to the forward TE, 2 and TE„ . The position of 
the guide axis is specified by x(z) and y(z), the coordinates of the guide 
axis in the transverse plane as functions of distance along the axis z. 
We will assume that x(z) and y(z) are independent Gaussian random 
processes with identical power spectra. 

The coupling coefficients ^'(f) and c [m] x (z) between TE,„ and the 
two polarizations of the m th spurious mode will also be independent 
Gaussian random processes, since the coupling coefficients are propor- 
tional to the corresponding geometric imperfections. Thus, for straight- 
ness deviations we have from Section 2.4 for small deviations 

cw'OO = C, w •**(*), c w x (z) = (",[,„]• </"(*), (308) 

where the symbols || and _L distinguish the two polarizations of the 
spurious mode, rather than x and y as in Section II. Since x(z) and y(z) 
are independent Gaussian random processes, x"(z) and y"{z) and 
consequently c [m j 1 (2) and c [m] x {s) will also be independent Gaussian 
random processes. 

The TE„i loss .1|,„| for both polarizations of the m ih spurious mode is 
given by 

•4,,,,] - .I,,,,, 11 +-4 [ „ ll x , (309) 

where A [m] l and A[ m ] X are independent random processes with statistics 
given by the results of Section 4.3.1. Writing 

A [m] = <A lm] > + 8A [m] , (310) 

we have for the expected value of the TE01 loss 

<,!,„„> = <,!,„,/> + Oi w x >. (311) 

Since 8A [m] 11 and 5.4 [,„]*, the two ac loss components, are independent, 
their power spectra and total powers add. Denoting the (continuous) 
power spectra of 8A [m] , 6A lm] l and 8A imf by P lmJ (v), P [m] ] Hv), and 
P[m\ X {v) respectively, we have 

/%,](") = Pim] l (v) +P t *] X ("), (312) 

<(8A lmi )' 2 > = <(8A [mi ] f> + <(8A lm] x ) 2 >. (313) 

Since both polarizations are assumed to have identical statistics, we 
have from Section 4.3.1 : 

D A r .. ( 1 



<A M ( \)> =L'S M [^\) = L.fl w (_). (314) 



'"•■V^QJ-U-M), |»l<l 



(315) 



«*„)•> = g-S W *(^ x) = £-S w '( 1 ) (316a) 
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p M (*)=f-/w(^x)-(i-M) 

r; 2 

v^^> - £*.@? x) - &*-(£) - ^(«*) 

In these and all subsequent results £ [mJ (f) is the power spectrum of 
each of the orthogonal components of the coupling coefficient. 

The power spectrum P w (r) of course remains triangular, and the 
3-db bandwidth and the corresponding interval in free-space wave- 
length or frequency remains as given in (305) and (306) for a single 
polarization of the spurious mode. These latter quantities are the same 
as those for the corresponding discrete case, given in (243a) and (244). 

Finally, since A [m] i and A [m] x are independent random variables with 
the same probability distribution [see (307)], the probability distribu- 
tion for A [m] is simply the convolution of (307) with itself. 

/ 4 \ 4yl[,„] 2A[ m \ /017', 

pUw) = rr>s exp ~7aT7>- (3l7) 

This result holds true for discrete tilts and offsets for zero differential 
loss. 

4.3.3 Many Spurious Modes 

For many spurious modes, the total TE i loss A is given as a sum over 
loss components A lm] due to the different spurious modes. From (284) 

A=2> (m] . (318) 

[m] 

Writing the total TE 0l loss as before, 

.4 = <A> + 5.4. (319) 

Then we have 

<A> = Zol[„„>, (320a) 

[m] 

BA = £aA w . (320b) 

The average TE M loss is simply the sum of the contributions of each of 
the spurious modes. 
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From ( 320b ) we have for the variance of the total TEoi loss 

<( 5 ,1) 2 > = ££<*4 W M W >. (321) 

[m] [n] 

The terms <(8A lm] ) 2 > are given by (316a) or (304a). It would be most 
convenient if the different ac components 8A [m] were independent, so 
that the cross terms <5.4 [ „ j] 5.4 [ „]> could be neglected. 

For the spurious modes produced by a geometric imperfection of a 
given angular symmetry (e.g., the TEi„ , produced by straightness 
deviations) the different 8A [m] are not independent. The TE i loss com- 
ponent due to the w th spurious mode at one frequency is proportional to 
the TEoi loss component due to the n th spurious mode at a widely sep- 
arated frequency, as in the discrete case (Appendix F). Thus, knowledge 
of one of the 8 A [ml over a sufficiently wide frequency band is sufficient to 
determine all of the others. However, at the same frequency 5.4 [ml and 
8A [n] are almost uncorrelated, so that the cross terms in (321) may be 
neglected. 

8A lm] and 8A M are easily seen to be almost independent in a simple 
way. Consider a frequency which corresponds to a sample point of the 
m th spurious mode. Under special conditions this frequency may also 
correspond to a different sample point of the » th spurious mode (in gen- 
eral, this will not be so). From (288b), .4 lm] will depend on only a single 
Fourier coefficient (say the k th ) of the geometric imperfection. Similarly, 
A [„] will depend on only a single Fourier coefficient, but on a different 
one (say the / th ), since different spurious modes have different beat 
wavelengths. Since the different Fourier coefficients of the geometric 
imperfection are almost independent, 8A [m] and 8A [n] will thus also be 
almost independent at this frequency. Since in general the sample points 
corresponding to different spurious modes do not precisely coincide, the 
correlation coefficient between 5.4 [m] and 5,4 [n] at a single frequency will 
not be identically zero, but should be small. 

The correlation coefficient between the ac components of the TE i loss 
due to two different spurious modes generated by the same type of 
geometric imperfection is derived in Appendix G for the special case in 
which the geometric imperfection and hence the coupling coefficients 
have white power spectra. Numerical results are given for the important 
practical case of TE 12 + and TE„ + generated by random straightness 
deviations, in which the second derivatives of the rectangular co-ordi- 
nates of the guide axis are independent Gaussian random processes with 
white spectra. The normalized correlation coefficient for reasonable guide 
lengths is very small indeed. 

The cross terms in (321) will consequently make only a negligible 
contribution to the variance of the total TEoi loss. Therefore, <(8A ) 8 > is 
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given simply by the sum of the contributions of the individual spurious 
modes; 

<(SA) 2 > = Z<(«iwA (322) 

Fi 

where the <(8A [m] ) 2 > are given by (316a) (or by (304a) for a spurious 
mode with a single polarization, e.g., TE 0m )- 

4.3.4 Discussion 

The TEoi Joss in a given frequency band, resulting from a given spur- 
ious mode, depends only on the Fourier components of the corresponding 
geometric imperfection for a narrow range of mechanical wavelengths. 
This band of mechanical wavelengths corresponds to the range of beat 
wavelengths between TE i and the spurious mode over the frequency 
band of interest. The statistics of the TEoi loss are strongly dependent 
on the power spectrum of the geometric imperfection. 

The present results are strictly valid only for zero differential loss, 
Aa = 0, although they will remain approximately true so long as the 
differential loss over the length of guide remains small, | Aa | L « 1 . 
However, further study shows that moderate values of differential loss 
Aa will change the average TEoi loss very little, but will smooth out the 
fluctuations of the TE i loss, for the present case in which the coupling 
coefficient power spectrum is essentially flat in the range of interest. 

The following sections will present specific; numerical examples for 
the various types of geometric imperfections. 

4.4 TEoi Loss Statistics for Random Straightness Deviations 

4.4.1 Introduction 

In the present section we apply the results of Section 4.3 to the case 
of random deviations from straightness of the guide axis. In Section 4.5 
other types of continuous geometric imperfections are similarly treated. 
There are two reasons for treating straightness deviations separately : 

1. Straightness deviations introduced by the manufacturing process 
are almost entirely responsible for the additional TEoi loss of present 
2-inch I.D. copper guide. 

2. Experimental TE 0i transmission measurements have provided a 
fair idea of the shape of the power spectrum of straightness deviations 
for different types of guide, at least over a limited range. This power 
spectrum differs in some respects from the power spectra that might be 
assumed for other types of geometric imperfections, in that under certain 
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conditions it contains an infinite low-frequency (i.e., long mechanical 
wavelength) content; some additional discussion of this particular case 
seems appropriate. 

We must first specify the statistical properties of the coupling coeffi- 
cient; since the coupling coefficient c(z) is assumed to be a stationary 
Gaussian random process, its statistics are completely determined by 
its power spectrum S(f). The continuous case is inherently more com- 
plicated to discuss than the discrete case of Section III. For the discrete 
ease only the rms offset, tilt, or diameter change must be specified. In 
the continuous case, however, the TE i loss statistics are no longer de- 
termined only by the mean square value of straightness deviation or 
other geometric imperfection; the shape of the power spectrum of the 
mechanical imperfection strongly influences the resulting TE i loss. We 
must therefore know the power spectrum of the imperfection before we 
can predict the TE i loss statistics of a guide. Conversely, knowledge of 
the TEoi loss statistics enables us to estimate the power spectrum of the 
imperfections. Up to now there have been no existing mechanical meth- 
ods for measuring the straightness deviation to the required accuracy so 
that its covariance and power spectrum can be determined;! TEoi trans- 
mission measurements have provided the only means of determining 
the significant Fourier components of the straightness deviation. 

Present experimental measurements of the TEoi loss over a band ex- 
tending from 33 kmc to 90 kmc, made by A. P. King and G. D. Mande- 
ville, indicate that for one type of 2-inch I.D. copper guide, the radius 
of curvature of the straightness deviation has an approximately flat 
power spectrum over the range of interest. Thus, if x(z) and y{z) are 
the rectangular components of straightness deviation with power spectra 
X{$) and F(f), we ha vej 

X(r) = F(f) « I (323) 

for mechanical wavelengths lying in the beat wavelength range for TEn 
and TE !2 (the most important spurious modes), 1.4 to 4.4 feet for the 
35-90 kmc band. 

The power spectrum of (323) for x{z) and y(z), which corresponds to 
a white power spectrum for the radius of curvature of the guide axis or 

t Methods of making these mechanical measurements are currently under de- 
velopment by K. J. Dahms, W. G. Nutt, and R. B. Ramsey and their associates 
at Bell Telephone Laboratories. 

}f is the "mechanical frequency," having the dimension . -7-; the corre- 
sponding mechanical wavelength is -. [See for example (291).] 
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equivalently for the second derivatives x"(z) and y"(z), appears plaus- 
ible under certain conditions, considering one way in which the guide 
has been made. If we imagine a guide made by drawing a copper pipe 
with more or less random hardness or wall thickness variations through 
a die, it is not hard to see that the radius of curvature of the guide axis 
might be a random process with a very short correlation distance, or 
equivalently with a very wide power spectrum. The power spectrum for 
x" (z) and y" (z) must, of course, fall off for sufficiently high mechanical 
frequencies (or sufficiently short mechanical wavelengths). 

The coupling coefficients and hence x"{z) and y"(z), the second de- 
rivatives of the displacement of the guide axis, have been assumed to be 
stationary random processes. The displacements themselves, x (z) and 
y(z), will not in general be stationary random processes, unless the power 
spectrum of x"(z) and y"(z) (and the corresponding coupling coeffi- 
cients) has special properties. However, this situation seems to be in 
accord with the physical facts. As a simple example we may consider a 
guide made of pipes with random uniform bows, screwed together at 
random; we might further assume that the first pipe of the guide starts 
out with zero displacement and zero slope, x(Q) = y(0) = x'(0) = 
y'(0) = o. Then it is obvious that while the second derivatives x"{z) 
and y"{z) and hence the coupling coefficients are stationary random 
processes, the displacements x{z) and y(z) are not. The variances of the 
displacements, <x 2 (z)> and <y 2 (z)>, grow with distance z; the guide tends 
to wander more and more from the axis as z increases, unless additional 
mechanical constraints are imposed in laying the guide. 

Since our knowledge of the power spectrum X(f ) is limited, any ex- 
ample chosen to illustrate the order of magnitude of straightness toler- 
ance that will have a significant effect on the TEoi transmission must be 
arbitrary to a considerable extent. For the present numerical example 
we assume that X(f ) and F(f ) are as given by (323). For this power 
spectrum the displacement of the guide axis x(z) or y(z) is not a sta- 
tionary random process, and the integral of (323) or the "total power" 
is infinite. However, the principal spurious modes are TEi 2 + and TEn + , 
with beat wavelengths of about 2.2 and 2.7 feet at a frequency of 55 
kmc. In order to get a rough measure of the short-wavelength straight- 
ness deviations that are responsible for the additional TE i loss, we shall 
quite arbitrarily calculate a "total" mean square straightness deviation 
<x 2 (z) > + <y 2 (z)> by including only those components having mechanical 
wavelengths less than 5 feet, f > j. While this is a rather arbitrary choice, 
it makes some physical sense. The significant components for TEn + and 
TEi2 + lie between 1.4 and 4.4 feet for a band from 35-90 kmc. The 
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long-wavelength components (greater than 5 feet) do not affect the 
TEm transmission in this hand; in any case, these long-wavelength com- 
ponents will depend strongly on the random errors made in laying the 
guide, and perhaps very little on the straightness deviations introduced 
by the manufacturing process. The mean square value of the components 
of wavelength less than 5 feet will give us a rough idea of the order of 
magnitude of the tolerance on the short-wavelength "manufacturing" 
straightness deviations. 

For the numerical example presented below we assume a 20-mile total 
guide length of 2-inch I.D. copper guide with equally spaced mode filters, 
spaced either 200 or 1000 feet apart. These mode filters are assumed to 
have zero loss to the TE i signal mode, infinite loss to the spurious TE im + 
modes. TM n + is neglected, since it has been adequately treated in Sec- 
tion III. The differential loss is neglected even though it is not small in 
the distance between mode filters for all spurious modes, particularly for 
the 1000-foot mode filter spacing. The differential loss will not greatly 
affect the average TE i loss, but will reduce the TE i loss fluctuations 
below the values computed for zero differential loss. The treatment of a 
long line with ideal mode filters is given in Section 3.4. 

The .c and y components of the straightness deviation of the guide 
axis are assumed to have power spectra given in (323), discussed above. 
The magnitude of these power spectra is chosen to yield an additional 
average TE 01 loss (due to all the propagating TE,„, + spurious modes) of 
I dl) per mile, at a frequency of 55 kmc. The total rms straightness devia- 
tion for components having wavelengths less than 5 feet is stated, for 
reasons discussed above. In addition, the rms straightness deviation for 
components having wavelengths between 2 and 3 feet, corresponding to 
the TE 12 + and TEn + beat wavelengths for the 50-60 kmc band, is also 
given. The contributions of each of the spurious modes to the average 
TEm loss and to the TE 0] loss fluctuations are stated separately. 

4.4.2 Analysis 

Let A"(f) and F(f) be the power spectra of the rectangular com- 
ponents of the straightness deviation of the guide axis, x(z) and y(z). 
Then 

X(f ) - [ x <x(z)x(z 4- r)> e' i2rST ch 

"I (324) 

F(r) - / <y(z)y(*+ r)>e- i2 ' iT dr. 

J— 00 
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We will always assume identical spectra for x(z) and y(z), i.e., 

X(t) = ¥({). (325) 

Then the power spectra of x"{z) and y"(z), the second derivatives of 
the x and y components of the displacement of the guide axis, are given 
by 

(2rf ) 4 X(r) = £ <x"{z)x"(z + r)> 6-** dr 

(27rf) 4 F(r) - r </(«)/(« + r)>«-** p dr. 

J— oo 

Noting (308), the power spectra »S[m](f) of the coupling coefficients 
c lm] (z) [see (290) and (291)] are given by 

S w '(rt =C, [m] 2 -(27rf) 4 Z(f), 2 

S W *(f) = C t[w] 2 -(2 fl -f) 4 F(r), 

for the parallel and perpendicular polarizations of the TE lm spurious 
mode respectively. 

Far from cutoff C l[mi is approximately inversely proportional to the 
free-space wavelength X. It is thus sometimes convenient to write C t [ m \ 
as 

C IM =^- ] , (328) 

A 

where ©<[„,] is now approximately independent of X. 

From (314) the average TEoi loss due to the TEi m spurious mode 
(with two polarizations) is 

<.4,„,,(x)> = Le, w a z> w V x (jj£r M- ( 329 ) 

The constant D [m ] , defined in (173), is related to the beat wavelength 
B lm] by 

A'"l -. _ 1 _ ffo ~ Plm) (330) 

2tt B[ m ] 2w 

In (329), X(f) is the power spectrum of each of the rectangular com- 
ponents of straightness deviation [see (325)]; the only frequency- (or 
wavelength-) dependent terms are the factors X 2 and the mechanical 

power spectrum xl -—- X J . The rms fluctuation of the TE (U loss com- 
ponent due to TEi„, about its average is given by (316b) j 

V<(8A [ml y> = -^§-. (331) 
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As discussed above, for the present example we assume that over the 
range of interest (mechanical wavelengths less than 5 feet) X(f), the 
power spectrum for each component of the straightness deviation, has 
the shape given by (32.3). Therefore we take 

X(f) = Yin = 7^;, (332) 

(2xf) 4 

where X is a scaling parameter determining the magnitude of the 
straightness deviation. Then the average loss of (329) becomes 

<A M (\)> = Le <w~ Xo = LC tM 2 X . (333) 

A" 

The loss fluctuation of course remains as given by (331). 

The total TE i average loss <A> and mean square loss fluctuations 
<(<5.4) 2 > are given by (320a) and (322) respectively, summing over the 
contributions of all the propagating TEi„, + spurious modes. We have 

<.4> = LX Z<S'« [ml '\ (334) 

[m] 

<(8A) 2 > = \lIXi Z^[».i 4 - (335) 

Substituting numerical values from Appendix A and summing over the 
9 propagating TEi m + modes, for a frequency of 55 kmc and a 1-inch 
guide radius (334) and (335) become 

<.4> = 111.24 LX (33G) 

<(5.4) 2 > = 3721.64 L 2 X„ 2 . (337) 

From (333) and (334) the average TE W loss <A> is inversely propor- 
tional to X 2 , or directly proportional to/". This is approximately in agree- 
ment with the experimental results of A. P. King and G. D. Mandeville 
for one type of copper guide 39 and provides the reason for the particular 
choice of power spectrum for the straightness deviation that has been 
made here. 

Finally, the rms straightness deviation in a given range is obtained by 
integrating the straightness deviation power spectrum over the appropri- 
ate range. For the .r-component alone, we have for the mean square 
straightness deviation for mechanical frequencies lying between f a and 

ft, , or equivalently mechanical wavelengths - lying between — and — 



£♦£] ■■ 



<x 2 > = + / X(f)rff, (338 



1136 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1962 

where we recall that the straightness deviation power spectrum X(f ) has 
been defined for both positive and negative mechanical frequencies. 
Substituting the particular power spectrum given in (332), 

< ;C 2 > = 2X \ fl-lV (339) 

3(2x) 4 \ra 3 fbV 

From (325) and the fact that x and y are independent, we have for the 
total mean square straightness deviation in this range 

<* + <*>'&(&-$)• (340) 

4.4.3 Numerical Example 

We assume an added average TE i loss in 2-inch I.D. guide due to 
mode conversion of 1 db per mile at 55 kmc. From Section 3.4 and (336), 
we determine X . 

X = ^r = 1-060 X 10 -7 ft" 1 . (341) 

111.24X8.6859X5280 

We will specify the rms straightness deviation (including both the x and 
y components of the displacement of the guide axis) for components 
having wavelengths less than five feet, as discussed above. From (340), 

•v/<.x- 2 > + <if> = 1-737 mils ; - < 5 feet. (342) 

For purposes of illustration we shall emphasize the 50-60 kmc band. 
In this band the beat wavelength range for TE 12 + and TE U + , the most 
important spurious modes, is 2-3 feet. Since these components of 
straightness deviation are the principal contributors to the TE i loss in 
the 50-60 kmc band, it is of interest to give the rms straightness devia- 
tion lying in the 2-3 foot region. Again from (340), 

V<x 2 > + <y 2 > = 0.677 mils; 2 feet < - < 3 feet. (343) 

Table VI presents the transmission behavior of a 20-mile, 2-inch I.D. 
guide for two mode filter spacings, 200 feet and about 1000 feet, calcu- 
lated from the present results and those of Section 3.4. In addition to the 
total average loss and rms loss fluctuation, the contributions of the 
individual spurious modes are given. The average loss is of course the 
same for both cases; the rms total loss fluctuation is reduced for the 
shorter mode filter spacing. 
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Table VI — TE i Loss Statistics for Straightness 
Deviation with a Flat Curvature Power Spectrum 

/-■..-! = 20 miles, total guide length. 
8.6859 -',-" — = 1 db/mile, additional average loss at 55 kmc. 



I-> total 



/ = 55 kmc, midband frequency. 
o = l inch, guide radius. 



Case 
1 


L 
M 

8.6859\/<(5.-ltoi a i) 2 > 


= 200 feet, mode filter spacing. 
= 528, number of mode filters. 
= 0.4773 db, rms total loss fluctuation 
line. 


or 20-mile 


Case 
2 


L = 996.23 feet, mode filter spacing. 
M = 106, number of mode filters. 
8.6859\/<(5-U olu i) 2 > = 1.0652 db, rms total loss fluctuation for 20-mile 
line. 






Case 1 


Case 2 


Spurious 
Mode 


Beat Wave- 
length Range 
for 50-60 Kmc 
Hand 


RMS 
Straightness 
Deviation in 
Beat Wave- 
length Range 


Added 

Average 

Loss 


RMS Total 
Loss Fluctua- 
tion for 20 
Mile Line 


A/ldb; 3-db 
Bandwidth of 

Power 
Spectrum of 
Loss Fluctu- 
ation 


RMS Total 
Loss Fluctu- 
ation for 20 
Mile Line 


A/3,11,: 3-db 
Bandwidth of 

Power 
Spectrum of 
Loss Fluctu- 
ation 


TE,i+ 
TEi2 + 
TE, 8 + 
TEu + 
TEi» + 


feet 

2.453-2.949 
1.997-2.408 
0.464-0.562 
0.215-0.263 
0.123-0.152 


mils 

0.513 
0.381 
0.043 
0.014 
0.006 


db/ 
mile 

0.2649 
0.7290 
0.0055 
0.0005 
0.0001 


db 

0.1630 
0.4486 
0.0034 
0.0003 
0.0001 


mc 

1486 

1212 

282 

131 

76 


db 

0.3638 
1.0012 
0.0076 
0.0007 
0.0001 


mc 

298 

243 

57 

26 

15 



4.4.4 Discussion 

The above results show that for the assumed mechanical power spec- 
trum, the principal contributors to both the additional average loss and 
to the loss fluctuations arise from the TEi 2 + and the TEn spurious 
modes, as has been observed experimentally. As discussed above, the 
average loss measured over a very wide frequency band by A. P. King 
and (i. D. Mandcville 39 provides the primary experimental data on the 
power spectrum of the straightness deviations. 

For some guides both the experimental transmission data and con- 
sideration of the manufacturing process indicate that the power spectrum 
of the second derivative of the straightness deviation should be more or 
less flat up to some high mechanical frequency (short wavelength). As 
discussed above, this leads to a power spectrum for the straightness 
deviation itself with infinite power, so that the displacement of the guide 
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axis from a perfect straight line will not be a stationary random process 
However, in practice additional forces are imposed in laying the guide, 
so that the straightness deviation must become stationary. In order for 
this to be so, the power spectrum for x"(z) and y"{z) and for the cor- 
responding coupling coefficients must now fall to zero as f approaches 0, 
at least as fast as f 4 . However, it seems reasonable to assume that the 
modification in the coupling coefficient power spectrum will take place 
only for very small values of f (long wavelengths). In the important 
spectral region corresponding to the TEn + and TE J2 + beat wavelengths, 

for - less than a few feet, the power spectrum should be modified very 

little; therefore, very little change will take place in the TE i transmis- 
sion statistics. 

Since only components of the straightness deviation having wave- 
lengths between about 1.4 and 4.4 feet will significantly affect the TE IU 
loss (in a band from about 35-90 kmc), it is clear that random straight- 
ness deviations arising in the laying of the guide, or manufacturing im- 
perfections such as long bows, will have very little effect on the TEoi 
transmission, because such straightness deviations will have their prin- 
cipal components at much longer wavelengths (e.g., perhaps greater 
than ten feet). One model that is readily analyzed is the "random bow 
line," made of pipes of identical length with uniform bows, screwed to- 
gether at random. The x and y components of the coupling coefficient 
are simply random square waves, whose power spectrum is well known. 
The allowable tolerance is several orders of magnitude more lenient than 
the tolerance on short wavelength straightness deviations. 

Finally, it is possible that (mite different types of power spectra than 
those discussed here could arise for different manufacturing processes. 
For example, a process that resulted in a periodic straightness deviation 
in the beat wavelength range would result in a rather broadly peaked 
band-pass power spectrum for the coupling coefficient. Such things are 
of course to be avoided. In practice, different manufacturing processes 
have produced quite different straightness deviation power spectra. 

4.5 TEoi Loss Statistics for Random Diameter Variations, Ellipticity, and 
High er-Order Deformations 

4.5.1 Introduction 

In the present section we apply the results of Section 4.3 to random 
diameter variations, ellipticity, and higher-order deformations of the 
cross-section of the guide, using Morgan's coupling coefficients E [fiml . 
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Random straightness deviations may of course also be treated in this 
way, but arc omitted since they have been discussed in Section 4.4. 

We take the same model for the guide as in Section 4.4, i.e., a 20-mile 
total guide length of 2-inch I.D. copper with equally spaced mode filters, 
spaced either 200 feet or about 1000 feet apart. The mode filters are as- 
sumed to have zero loss for the TEoi signal mode, zero loss for the TE ,„ 
spurious modes, and infinite loss for all other TE„„, spurious modes. 
The differential loss is assumed small in the distance between mode filters 
for all spurious modes; in addition, for the TE 0m spurious modes the 
present analysis forces us to assume that the differential loss is negligible 
for the total guide length, 20 miles in the present example. 

The present analysis should provide a reasonable approximation for 
ellipticity and for higher-order deformations. The differential loss in these 
cases will not affect significantly the average TE i loss, but will reduce 
the TEoi loss fluctuations somewhat below the values computed for zero 
differential loss in those cases where the differential loss is not completely 
negligible in the distance between mode filters. 

However, for diameter variations the TE 0m differential loss is certainly 
not negligible in 20 miles, as required by the analysis. While this ap- 
proximation will not lead to an appreciable error for the average TEoi 
loss, it will certainly lead to serious error for the mean square TE i loss 
fluctuation, which is the really significant quantity, and for A/ 3l ib- The 
actual TEoi loss fluctuations will be much smaller than those computed 
here. 

The power spectrum for straightness deviations is known, at least ap- 
proximately, over a moderately wide range from TEoi transmission 
measurements, as discussed in Section 4.4. A very elementary considera- 
tion of the manufacturing process suggests the same shape for this 
power spectrum as is actually observed in certain cases. Unfortunately, 
for other types of deformation neither of these approaches has suggested 
the proper form for the power spectrum. Spurious modes, other than 
TE !2 + and TE n + , have not been observed at high enough levels to permit 
an estimate of the power spectrum of the corresponding mechanical im- 
perfection to be made. As yet, no simple picture of the manufacturing 
processes has yielded a guess as to the shape of the power spectra. About 
all that can be said is that the power spectra must eventually fall off 
in some manner for high enough mechanical frequencies f (short enough 
mechanical wavelengths), since the mean square values of the various 
cross-sectional tolerances are certainly bounded. 

Since the power spectra of the various deformations are not known, 
the numerical examples presented below are less specific than the example 
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for straightness deviation in Section 4.4. For example, ellipticity gener- 
ates the TE 2m + spurious modes. For each of these spurious modes the 
rms ellipticity in the mechanical wavelength range corresponding to the 
beat wavelength for the 50-60 kmc band is chosen to yield an additional 
average TEoi loss of 1 db per mile, assuming a flat power spectrum in 
this range. The same is done for trifoil and higher-order deformations. 

Diameter variations are treated separately, partly because these re- 
sults cannot be taken too seriously, as discussed above, and partly be- 
cause the equations are somewhat different for this case. Here for each 
of the spurious TE 0m + modes we choose the rms diameter variation in 
the mechanical wavelength range corresponding to the beat wavelength 
for the 50-60 kmc band to yield a 1-db rms loss fluctuation for the TE i 
loss component due to the spurious mode, for the entire 20-mile line. 

4.5.2 Random Diameter Variations 

From (179) to (181) the radius of the guide is given by 

r = a + ao(z). (344) 

From (207) the coupling coefficient is given by 

C[0m](3) - — E[Om]Oo(z). (345) 

Let 2JLo(f ) be the power spectrum of a (z), i.e., 

jr Q ( f ) = f <a (z)a Q (z + r)> e~ i2 ' iT dr. (346) 

J— CO 

Then the power spectrum S [0m ]U) of the coupling coefficient c [0m] (z) is 
WD - H [0m] 2 -^o(r). (347) 

Remembering that the spurious mode now has only a single polariza- 
tion, we have from (302) and (304b) 

<A [0m] > = V<(8A [0m] y> = \ W'^ x ) , < 348 ) 

where the constant D [0m ] is related to the beat wavelength B l0m ] by 

Aoml ^ _ 1 _ 001 — 00m ^ (349) 

2tt B[o m ] 2ir 

The coupling coefficient E[om] is approximately proportional to the free- 
space wavelength X. 

Finally, the mean square radius variation in the range of mechanical 

frequencies f a to f & ( or the range of mechanical wavelengths - from — 
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to — ) is given by 



<a 2 (z)> = 



:/;>£•] 



tfo(r)#, 



(350) 



again integrating over positive and negative frequencies. Assuming that 
the power spectrum is flat over the range of integration, 

A'o(r) =K , f„ < If I < ft, (351) 

(350) becomes 

<a 2 (z)> = 2Ko-(fi- fa). (352) 

These formulas are used to calculate the results in Table VII. 
The present analysis could alternately have been carried out in terms 
of the coupling coefficients dim] , as discussed in Section 2.3.9. 

Table VII — Radius Variation Yielding 1 db RMS TE i Loss 
Fluctuation and 1 db Average TE i Loss for the Loss Com- 
ponent A[ 0m ] Due to Each of the TE ,„ + Spurious Modes, at 

55 KMC 

Z/totni = 20 miles, total guide length. 
8.685 ( J\/<(5'l|omi) 2 > = 1 db, rnis loss fluctuation for each component, for 20-mile 
line. 
8.0859C1 [om i> = 1 db, additional average loss for each component, for 
20-mile line. 
/ = 55 kmc, midband frequency. 
a = 1 inch, guide radius. 
Differential loss assumed small over total guide length of 20 miles. 









A'o : Spectral Density of 
Radius Variation 


A/3db: 3-db Bandwidth of 


Mode 


Range for 50-60 
Kmc Band 


tion in Beat Wave- 
length Range 


Power Spectrum of Loss 
Fluctuation 




feet 


mils 


mils 3 / foot-' 


mc 


TEo2 + 


0.7886-0.9532 


0.087 


1.7291 X 10- 2 


0.907 


TE 03 + 


0.3003-0.3654 


0.097 


0.7942 X lO- 2 


0.347 


TEo4 + 


0.1588-0.1953 


0.102 


0.4397 X 10~ 2 


0.185 


TE 06 + 


0.0964-0.1206 


0.105 


0.2671 X 10- 2 


0.113 


TE 06 + 


0.0627-0.0805 


0.109 


0.1691 X 10-2 


0.075 


TE 7 + 


0.0418-0.0563 


0.115 


0.1064 X 10-2 


0.051 


TE 08 + 


0.0259-0.0402 


0.129 


0.0609 X 10-2 


0.035 



4.5.3 Random n-foils 

We now consider random cross-sectional deformations of higher order. 
For an "n-foil" Ihe radius of the guide is given by (179) to (181) as 



r — a + ci„(z) cos n<p + b„(z) sin rup. 



(353) 
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The n = case corresponds to diameter variations, studied in Section 
4.5.2. The n = 1 case corresponds to straightness deviations, studied in 
Section 4.4; the quantities ai(z) and bi(z) are equal to x(z) and y(z), 
the rectangular components of straightness deviation. The n = 2 case 
corresponds to ellipticity, the n = 3 case has been designated as "tri- 
foil," etc. All formulas in the present section hold only for n ^ 1; the 
n = case has been treated in the preceding section. 

The magnitude of the n-foil distortion at a given position z along the 
axis is specified by the maximum departure from a perfect guide, r = a. 
From (353) we have for an n-foil 



\r - a | max = Va n 2 (z) + b,;{z). (354) 

For diameter variations this definition yields the change in guide radius; 
for straightness deviations it yields the total displacement (in the x-y 
plane) of the guide axis. For ellipticity and higher-order deformations, 
(354) yields the maximum deviation from a perfect circle r = a. Note 
that for n = 2 this is only one-quarter as large as a commonly accepted 
definition of ellipticity, the maximum diameter minus the minimum 
diameter. 

The coupling coefficients for the two polarizations of each of the 
spurious modes are given from (207) by 

C[nm]"(2) = -E[n«]a„(2), (355) 

c lnm] x (z) = -S[n*]b.(«). (356) 

Let K H (t) be the power spectrum of each of the two components a„{z) 
and b n (z), i.e., 



v (>\ - fT <a »toa.(«+r)> 

KAn ~ LJ<b n (z)b n (z+T)>j 



e~ iUiT dr. (357) 



Then S [nm] ({), the power spectrum of each of the two components of 
the coupling coefficient, is given by 

S [nm] (t) = S [flm] 2 -if B (f)- (358) 

Using the results of (314) and (316b) for two polarizations, we have 

<A [nm] > = V2<(5A lnm] y-> = L'S lnm] (^f \j , (359) 

where 

glz^X = -J— = fa ~ /3 "'". (360) 

2lT B[ nm \ 2lT 
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The coupling coefficients S[„,„] are approximately proportional to the 
free-space wavelength A. 

Finally, the mean square n-foil magnitude (denned as the maximum 
departure from perfect circularity at a given cross section) in the range 

of mechanical frequencies ft, to ft, ( or the range of mechanical wave- 
lengths - from — to - J is given from (354) by 

<\r - ajLx> = <a,;(z)> + <h» 2 (z)> 

(361) 



2 |~( ""+ [\KMdt. 



Assuming that the power spectrum is flat in the range of integration, 

< I r - a |Lx> = 4/C„- (ft, - ft.). (362) 

These formulas are used to calculate the results in Tables VIII through 
XII for values of n ranging from 2 to 6. The case of straightness devia- 
tions, n = 1, is omitted; the same results as those of Section 4.4 would 
of course be obtained. 

4.5.4 Discussion 

Neither TE U i transmission measurements nor mechanical measure- 
ments have thus far yielded information on the shape of the power 
spectra for diameter variations, ellipticity, and higher-order cross-sec- 
tional deformations. This is true principally because these effects are 
very small in present guides. We are thus unable to predict what the 
relative contributions of the various spurious modes to the average 
TEoi loss and to the TE i loss fluctuation might be. 

In the numerical examples presented in this section, the spectral 
density of the geometric imperfection has been chosen in such a way 
that for each type of imperfection the contributions of each of the spuri- 
ous modes are equal. Plots of the logarithm of the spectral density K n 
vs the logarithm of the mechanical frequency f show that in each case 

K„(f) falls off approximately as -. Thus, if the spectral density of the 

imperfection is flat, the higher modes will become progressively more 
important as the second mode index increases. This is quite different 
from the observed behavior for straightness deviations, discussed in 
Section 4.4, where only TEi 2 + and TEn + have an appreciable effect on 
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Table VIII (n = 2) — Ellipticity Yielding 1 db/mile Addi- 
tional Average TEoi Loss for the Loss Component A [2m ] Due 
to Each of the TE 2 ,„ + Spurious Modes, at 55 kmc 



Ltotni = 20 miles, total guide length. 
g 685Q < — i ,2m) total = j db/mile, additional average loss for each 



Ltotnl 



component. 
/ = 55 kmc, midband frequency. 
a = 1 inch, guide radius. 



Case 
1 



Case 
2 



L = 200 feet, mode filter spacing. 
M = 528, number of mode filters. 

}859\/<(oVl [2mitotai) 2 > = 0.6155 db, rms total loss fluctuation for each 
component, for 20-mile line. 



L = 996.23 feet, mode filter spacing. 
^ M = 106, number of mode filters. 

8.6859\/<(^4 [2,„]totui) 2 > = 1.3736 db, rms total loss fluctuation for each 

component, for 20-mile line. 



Spurious 
Mode 



Beat Wavelength 

Range for 50-60 

Kmc Band 



TE 21 + 
TE22 4 
TE 23 + 
TE 24 + 
TE 2& 4 
TE 26 + 
TE 27 + 
TE 28 + 



feet 

5.1623-6 
0.9006 4 
0.3154-0 
0.1631-0 
0.0982-0 
0.0636-0 
0.0423-0 
0.0263-0 



.2110 
.0880 
.3836 
.2005 
.1226 
.0815 
.0569 
.0405 



RMS Ellipticity 
in Beat Wave- 
length Range 



0.3737 
0.5142 
0.5992 
0.6349 
0.6603 
0.6861 
0.7216 
0.8103 



Spectral 
Density of 
Ellipticity 



milsVfoot-' 

1.0676 
0.3457 
0.1592 
0.0882 
0.0536 
0.0340 
0.0214 
0.0123 



Case 1 



A/3cll>: 3-db Band- 
width of Power 
Spectrum of Loss 
Fluctuation 



A/3db: 3-db Bandwidth 

of Power Spectrum of 

Loss Fluctuation 



3128 
547 

102 
100 

(il 

40 
27 
19 



Case 2 



628 

110 
39 
20 
12 



5.5 

3.8 



the TE01 loss; there, however, the straightness deviation power spectrum 
1 

In any practical case the power spectrum K n {$) of the geometric 
deformation must eventually fall off as f increases (for n s^ 1). If the 
derivative of the imperfection exists (which seems a reasonable require- 

1^ 

r 2 

of course the real question is how large f must be for this behavior to 
dominate. The higher spurious modes have very short beat wavelengths 
— in the range of an inch or less; if the power spectrum K n (\; ) has begun 
to fall off appreciably at such wavelengths, the contribution of the higher- 
order spurious modes will be small. 



falls off very rapidly, as — . 



ment), K n (t ) must eventually fall off faster than — as f becomes large; 
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Table IX (n = 3) — Tri-Foil Yielding 1 db/mile Additional 

Average TEoi Loss for the Loss Component A [3m ] Due to 

Each of the TE,s to + Spurious Modes, at 55 kmc 



8.6859 



Z< total = 20 miles, total guide length. 

•=1 db/mile, additional average loss for each 



i»lotul 



component. 
/ = 55 kmc, midband frequency. 
o = l inch, guide radius. 



Case 
1 


8.6859-v / <(5A [ 3„ 


L = 200 feet, mode filter spacing. 
M = 528, number of mode filters. 
ltoini) 2 > = 0.6155 db, rms total loss fluctuation for each 
component, for 20-mile line. 


Case 


L = 996.23 feet, mode filter spacing. 
M = 106, number of mode filters. 


2 


8.6859V<(5.'1[ 3 „ 


ltoui) 2 > = 


.3736 db, rms total loss fluctuation for each 
component, for 20-mile line. 




Beat Wavelength 

Range for 50-60 

Kmc Band 


RMS Tri-Foil 
in Beat Wave- 
length Range 


Spectral 

Density of 

Tri-Foil 


Casel 


Case 2 


Spurious 
Mode 


A/31I ij : 3-db Band- 
width of Power 
Spectrum of Loss 
Fluctuation 


A/3 tib: 3-db Band- 
width of Power 
Spectrum of Loss 
Fluctuation 


TE 3 .+ 
TE 32 + 
TE 8S + 
TEi4 + 

TE SB + 
TE 3(i + 
TE 37 + 


feet 
9.2836-11.1800 
0.5464-0.6617 
0.2315-0.2827 
0.1285-0.1590 
0.0801-0.1011 
0.0526-0.0687 
0.0346-0.0485 


mils 
0.1876 
0.5364 
0.6045 
0.6389 
0.6661 
0.6964 
0.7449 


milsVfoot -1 

0.4817 
0.2257 
0.1168 
0.0684 
0.0428 
0.0272 
0.0167 


mc 

5628 

332 

141 

79 

50 

34 

23 


mc 
1130 

67 

28 

16 

10 
6.7 
4.6 



For cllipticity and higher-order deformations, the additional average 
loss in db per mile and the total loss fluctuation for a 20-mile guide are 
roughly comparable, for reasonable mode filter spacings. Thus if the 
additional average loss is small, as it is in present guides, the loss fluctua- 
tion will also be small. 

For diameter variations, the total additional average loss for a 20-mile 
guide is roughly comparable with the total loss fluctuation. Here the 
total loss fluctuation can remain serious even though the additional av- 
erage loss in db per mile remains small, as it does in present guides. 

The results for diameter variations given here are pessimistic, since 
the differential loss to the TE ,„ + spurious modes is neglected over the 
total guide length of 20 miles in the present analysis. The loss fluctua- 
tions in practice will be much smaller than those given here. The results 
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Table X in = 4) — 4-Foil Yielding 1 db/mile Additional 

Average TE i Loss for the Loss Component A [im] Due to 

Each of the TE 4m + Spurious Modes, at 55 kmc 



8.6859 



<A 



Lt„t. 



1< total = 20 miles, total guide length. 

1 = 1 db/mile, additional average loss for each 
component. 









/ = 
a = 


55 kmc, midband frequency 
1 inch, guide radius. 




Case 
1 


8.6859\/<(oA 


L = 

M = 

[4m]tutul) 2 > = 


200 feet, mode filter spacing. 
528, number of mode filters. 
0.6155 db, rms total loss fluctuation for each 
component, for 20-mile line. 


Case 


8. 




L = 
M = 


996.23 feet, mode filter spacing. 
106, number of mode filters. 


2 


5859 V< (S. 1 


4m]total) 2 > = 


1.3736 db, rms total loss fluctuation for each 
component, for 20-mile line. 






RMS 4-Foil in 
Beat Wave- 
length Range 


Spectral 

Density of 

4-Foil 


Casel 


Case 2 


Spurious 
Mode 


Bea 
Ra 


t Wavelength 
nge for 50-60 
kmc Band 


Ahdb'- 3-db Band- 
width of Power Spec- 
trum of Loss Fluctu- 
ation 


A/3,|l ( : 3-db Band- 
width of Power Spec- 
trum of loss Fluctua- 
tion 


TE 41 + 
TE* a + 
TE„ + 
TE« + 

TE 4 = + 
TE 46 + 

TE„ + 


feet 
2.0190-2.4344 
0.3757-0.4562 
0.1782-0.2186 
0.1039-0.1295 
0.0663-0.0847 
0.0438-0.0586 
0.0275-0.0416 


mils 

0.2991 
0.5433 
0.6059 
0.6414 
0.6719 
0.7091 
0.7890 


mils'/foot -1 

0.2647 
0.1571 
0.0885 
0.0540 
0.0344 
0.0218 
0.0127 


mc 

1225 

229 

109 

64 

42 

28 

19 


mc 

246 

46 

22 

13 
8.4 
5.7 
3.9 



for other cross-sectional deformations are of course valid, since the differ- 
ential loss will be small in the short distance between mode niters. 



4.6 Conclusions 

The TEoi loss of a long waveguide has been treated as a random proc- 
ess, and its statistics determined in terms of the statistics of the geo- 
metric imperfections. A statistical analysis is necessary because it would 
be impractical to make mechanical measurements of the complete geom- 
etry of any great length of guide, even if this were possible. 

The numerical results show that rms tolerances of the order of 1 mil 
are required, for any of the various types of imperfections, to yield an 
additional average TE i loss of the order of 1 db per mile. The rms di- 
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Table XI (n = 5) — 5-Foil Yieldixg 1 db/mile Additioxal 

Average TE i Loss for the Loss Component A[ 6m ] Due to 

Each of the TE 6 „, + Spurious Modes, at 55 kmc 



Ltotni = 20 miles, total guide length. 



8.6859 



i-tul.ll 



= 1 db/mile, additional average loss for each 
component. 
/ = 55 kmc, midband frequency. 
a = l inch, guide radius. 



Case 




L = 

M = 


200 feet, mode filter spacing. 
528, number of mode filters. 




8.6859\/<M 


Sm)totnl) 2 > = 


0.6155 db, rms total loss fluctuation for each 
component, for 20-mile line. 


Case 




L = 

M = 


996.23 feet, mode filter spacing. 
106, number of mode filters. 


2 


8. 6859 V< (6 A 


5„.,.ot nl ) 2 > = 


1.3736 db, rms total loss fluctuation for each 
component, for 20-mile line. 




Beat Wavelength 

Range for 50-60 

kmc Band 


RMS 5-Foil in 
Beat Wave- 
length Range 


Ks 

Spectral 

Density of 

S-Foil 


Case 1 


Case 2 


Spurious 
Mode 


■i/3«ll>: 3-db Band- 
width of Power Spec- 
trum of Loss 
Fluctuation 


&fidb'- 3-db Band- 
width of Power Spec- 
trum of Loss 
Fluctuation 


TE 5l + 
TE 52 + 

TE H + 
TE M + 
TE 55 + 
TE 58 + 


feet 

1.0317-1.2458 
0.2771-0.3375 
0.1417-0.1748 
0.0855-0.1076 
0.0554-0.0719 
0.0363-0.0503 


mils 

0.3297 
0.5439 
0.6057 
0.6436 
0.6787 
0.7269 


milsVfoot -1 

0.1632 
0.1144 
0.0687 
0.0432 
0.0277 
0.0172 


mc 

626 
169 
87 
53 
35 
24 


mc 

126 
34 

17 

11 
7.1 
4.8 



ameter variation of present copper waveguide is of the order of 0.1 mil. 
Consequently the diameter, ellipticity, tri-foil, and higher-order defor- 
mations must have comparable or smaller tolerances, so that they should 
have a negligible effect on present TEoi transmission measurements. 
Experimental observations support this conclusion. The only spurious 
modes ever observed in measurements on relatively short (i.e., a few 
hundred feet in length) waveguides are TE 12 + and TE U + , arising from 
straightness deviations. TE 0m + , TE 2m + , TE 3m + , and higher TE nm + modes 
have never been observed with significant magnitudes. 

Straightness deviation is the one tolerance about which in the past 
we have had no information at all, via mechanical measurements; it is 
the only significant tolerance in present measurements. Random straight- 
ness deviations are believed to account for substantially all of the addi- 
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Table XII (n = 6) — 6-Foil Yielding 1 db/mile Additional 

Average TE 0i Loss for the Loss Component A [6m] Due to 

Each of the TE 6m + Spurious Modes, at 55 kmc 



8.6859 



<A 



i'tolnl 



Ltotai = 20 miles, total guide length. 

— = 1 db/mile, additional average loss for each 

component. 
/ = 55 kmc, midband frequency. 
a = 1 inch, guide radius. 



Case 




L = 
M = 


« 200 feet, mode filter spacing. 
= 528, number of mode filters. 


1 


8. 6859 V< (« A 


6m]total) 2 > = 


= 0.6155 db, rms total loss fluctuation for each 
component, for 20-mile line. 


Case 




L = 

M = 


- 996.23 feet, mode filter spacing. 
= 106, number of mode filters. 


2 


8.6859V<(5.-1 


6m]totnl) 2 > = 


= 1.3736 db, rms total loss fluctuation for each 
component, for 20-mile line. 




Beat Wavelength 

Range for 50-60 

kmc Band 


RMS 6-Foil 

in Beat Wave- 
length Range 


Spectral 

Density of 

6-Foil 


Case 1 


Case 2 


Spurious 
Mode 


A/adb : 3-db Band- 
width of Power Spec- 
trum of Loss 
Fluctuation 


A/adl>: 3-db Band- 
width of Power Spec- 
trum of Loss 
Fluctuation 


TE„+ 
TE, a + 
TE M + 
TE 64 + 
TE 6r> + 
TEoc + 


feet 

0.6532-0.7902 
0.2138-0.2614 
0.1154-0.1433 
0.0714-0.0907 
0.0465-0.0617 
0.0295-0.0434 


mils 

0.3394 
0.5418 
0.6047 
0.6461 
0.6873 
0.7582 


mils'Vfoot -1 

0.1085 
0.0861 
0.0543 
0.0349 
0.0223 
0.0133 


mc 

397 
131 
71 
45 
30 
20 


mc 
80 
26 
14 

9 

6 

4.1 



tional loss observed in present copper waveguides, and to account for a 
substantial part of the additional loss in helix as well. 

Diameter variations in both copper and helix guide give rise to the 
TE 0m spurious modes, which cannot be satisfactorily attenuated by ex- 
isting structures. It may not be sufficient to have a good enough diam- 
eter tolerance to yield a small additional average TE i loss, since the 
TEoi loss fluctuations may still remain objectionable. The diameter 
tolerance must thus be substantially better than the tolerance for ellip- 
ticity, tri-foil, etc., in copper guide. However, the present results for 
diameter variations are too pessimistic, since the differential loss over 
the entire length of guide was neglected. 

Tolerances on diameter, ellipticity, tri-foil, etc. for drawn copper 
guide are tolerances at a single cross section, and are controlled pri- 
marily by the accuracy of the die through which the guide is drawn. A 



REFERENCES 1149 

good straightness tolerance requires accurate alignment between differ- 
ent cross sections separated by substantial distances (equal to the beat 
wavelength of TK 12 + , TEn + ). This depends on many factors other than 
the dimensional accuracy of the dies; for example, random variations in 
hardness or wall thickness may cause the axis of the guide to curve as 
it is being drawn. Consequently, a good straightness tolerance is more 
difficult to attain than any of the other cross-sectional tolerances, in 
drawing copper guide. Numerous other manufacturing processes for 
copper guide are currently under study. 42 

The variation of the TE i loss statistics with mode filter spacing is 
illustrated in the examples given above. The average loss is unaffected 
by the mode filter spacing, but the rms loss fluctuation is inversely pro- 
portional to the square root of the number of mode filters. As discussed 
in Section 3.4, the total TE i loss will be approximately a Gaussian 
random process. 

While the present analysis applies only to copper waveguide with the 
differential loss neglected, further study 41 shows that adding loss to the 
spurious mode has an effect similar to that of increasing the number of 
mode filters; for moderate values, as Aa increases the average loss changes 
very little, while the loss fluctuation will be progressively reduced. Thus, 
accurate tolerances will be important in helix (or in copper guide with 
a lossy dielectric lining) as well as copper guide, although in helix the 
principal effect of poor tolerances will be an increased average loss, the 
loss fluctuations remaining small (diameter variations excluded). 

Finally, the shape of the power spectra of the different mechanical 
imperfections is all-important in determining the resulting TE i loss due 
to mode conversion. Only components of the mechanical imperfection 
in the beat wavelength range of the important spurious modes have any 
effect on the TE i transmission. The short-wavelength straightness devi- 
ations "built in" to the guide in the manufacturing process will be prin- 
cipally responsible for additional loss due to mode conversion to modes 
such as TEi2 + and TEn + in copper guide; long bows or random straight- 
ness deviations due to imperfect laying of the guide will have only a 
very small effect. 
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APPENDIX A 

Coupling Coefficients for Tilts, Offsets, and Diameter Changes 1 *' 16 

General formulas for the coupling coefficients from TEm to the first- 
order spurious modes at offsets, tilts, and diameter changes, as deter- 
mined by S. P. Morgan, are given in Table XIII. Numerical values are 
tabulated at 55 kmc for 2-inch diameter guide in Tables XIV, XV and 
XVI; the computations here and in Appendix D were performed by Mrs. 
C. L. Beat tie. 

Notation 

C„[ m ] ± — coupling coefficient between TE i and forward ( + ) or back- 
ward ( — ) TEi« for an offset in copper guide. 

Ctim] — coupling coefficient between TE (ll and forward ( + ) or back- 
ward ( — ) TEi m for a tilt in copper guide. 

(',u]) + — coupling coefficient between TE 0J and forward TMn for a tilt 
in copper guide. 

Cdim* — coupling coefficient between TE i and forward ( + ) or back- 
ward ( — ) TEo„, for a diameter change in copper or helix guide. 

a — guide radius. 

A — free-space wavelength. 

/>•„„, — Bessel root given by J n '(k„ m ) = 0. 

"„,„ = -zr^- , cutoff factor for the TE„„, mode. 
'lira 

0„m = y" \/l — v„,„- , propagation constant for the TE„ m mode. 

A 
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Table XIII 



Coupling Coefficient 



Equation 



1 



kmki n 



V2 a (*!„' - fcoi 2 ) Vfcim 2 -1 V/Soi/5im 
a k ik lm 2 (/3oi ± <?i,„) 2 



" ml " V2 (&oi 8 - AwO'VW-l Vftift 
Ci(ii) 



+ _ V^TT d 



A'oi^ 

± _ 1 A-Ql^Dm 00m ± fill - j 

d|ml a Ar 0m 2 - hi 2 -\/p oi f5 0m ' 




Table XIV — Coupling Coefficients for Offset 
(/ = 55 kmc ; a = 1 inch) : 



Spurious Mode 


c olm] + 


COM' 


TE11 
TE U 

TE.3 
TEu 
TE 15 
TE I6 
TE 17 
TE, 8 
TE 19 


inch -1 

-1.052295 
2.140300 
0.800615 
0.520766 
0.392376 
0.317601 
0.269338 
0.238323 
0.230068 


inch -1 

0.0035023 
0.0087996 
0.0143221 
0.0204174 
0.0274958 
0.0362493 
0.0480943 
0.0668142 
0.1105813 



Table XV — Coupling Coefficients for Tilt 
(/ = 55 kmc; a = 1 inch) 



Spurious Mode 


C([ml + 


C,[m| 


TE„ 
TE12 
TE 13 
TEu 
TE 16 

ti-;,o 

TEit 
TEu 
TE 19 


radian -1 
5.428116 
9.004674 
0.784756 
0.237781 
0.103215 
0.0534059 
0.0306226 
0.0187490 
0.0121822 


radian -1 

0.000060129 

0.000152211 

0.000251132 

0.000365504 

0.000506837 

0.000695704 

0.000976413 

0.00147362 

0.00282165 




Cj(n) + 


Ci(u>~ 


TMn 


radian -1 
5.403 


radian -1 

0. 
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Table XVI — Coupling Coefficexts for Diameter Change 
(/ = 55 kmc ; a = 1 inch) 



Spurious Mode 


G*[.]+ 


C d \m\- 




inch -1 


inch -1 


TE 02 


1.556796 


-0.0162900 


TE 3 


0.878135 


-0.0244495 


TEo4 


0.627935 


-0.0336831 


TEo6 


0.493929 


-0.0447267 


TEo« 


0.410428 


-0.0589493 


TEo, 


0.355528 


-0.0795021 


TEob 


0.324255 


-0.1168700 


TEog 


0.425098 


-0.329569 
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Geometry of Discrete Tilts 

Let unit vectors directed along the guide axes in the two guide sec- 
tions adjacent to a discrete tilt be ti and t 2 . Then 



ti = i#i + jiji +kzi; | ti | - 1. 
ta = i-i'2 + m r- to* ; 1 1 2 1 = 1 . 



(B-l) 
(B-2) 



i, j, and k are unit vectors along the x, y, and z axes respectively. Since 
ti and ta are unit vectors, 



•'V 4- }i\ + z\ = 1, 

•<V' + !h + z-i = 1. 



(B-3) 

(B-4) 

Let the tilt have angle a, orientation 6 as defined in Section 2.1.2; fur- 
ther let the corresponding angles of the projections of the guide axes 
on the x-z and y-z planes he a x and a y , as in (38). Then by taking ap- 
propriate dot products we have: 



cos a = X X X% + 7/,.(/a + Z&, 
X1X2 + Zxz-i 



cos a, = 



cos a u = 



\/xi 2 + Zi 2 Vx 2 2 + Z? ' 



Ihlte 



Z\Zl 



V2/i 2 + z{- \/yr + z£ 



(B-5) 
(B-G) 

(B-7) 



where (B-3) and (B-4) of course hold true. If the angular deviation 
of the guide axis from the z-axis is small, we have 



xi «1, /yi«l, x, « 1, 



2/2 « 1. 



(B-8) 
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Under these conditions (B-5) to (B-7) yield to first order 

a « V(xo - a:,) 2 + (1J2 - Ui)\ ( B " 9 ) 

a x » x. 2 - Xl , (B-10) 

a u tty 2 -yi, (B-ll) 

where we have made use of (B-3) and (B-4). 
Now define the unit vector p as follows: 



U - tx tt- ti 

P ~ I t 2 - t! I 2 sin (a/2) 



(B-12) 



p lies in the plane of the tilt, i.e., in the tp = 6 plane, and bisects the 
angle made by the guide axes on the two sides of the tilt. The vectors 
t = (t x -f t 2 )/2 and p are the analogs for the discrete case of the tangent 
and principal normal vectors of differential geometry, introduced in 
Section 2.4, in the treatment of the continuous case. Denote the trans- 
verse component of p by p x , v ; then by (B-12), (B-l) and (B-2) 

*- = J^fei) lite " *> + ,( » - »» (B " ,3) 

Now if the angular deviation of the guide axis from the z-axis is small 
(on both sides of the tilt), as assumed in (B-8), and if a unit vector 
perpendicular to the guide axis and lying in the tp = plane is almost 
parallel to the z-axis, then the angle of p*.„ with respect to the .r-axis 
will be approximately equal to the orientation of the tilt; under these 
conditions we have from (B-13) 

cosfl^ /, X * ' l \ v 8 , ( B " 14 > 

■ Q ^ ——= V* ~ yi (B-15) 

B1D d ~ V(X2 - xi)« + (2/2 - i/i) 2 

From (B-9) to (B-ll), (B-14) and (B-15) we then have 

acosd ^ a x , (B-16) 

a sin 6 tt a u , (B-17) 

as stated in (38). 



(CM) 
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APPENDIX C 

Energy Relations for Guides with Real Coupling Coefficients 
Consider the coupled line equations given in (85) 
h'(z) = -TMz) +jc{z)h(z) 

h'(z) = jc(z)I (z) - TMz) 

To = «o + jfo ; r! = ai + jfa . (C-2) 

We assume in this appendix that c(z) is pure real; 

Im c(z) = 0. (C-3) 

Consider first the ease of ideal metallic guide, for which a„ = ai = 0- 
The total power P(z) flowing in the guide at the point z is simply 

P(z) = | h(z) | 2 + | l t {z) |- = h{z)h*{z) + h(z)h*(z). (C-4) 

Now in general, for guide whose walls are not perfect conductors, the 
a's will not be identically zero, and neither (C-3) nor (C-4) hold true. 
Helix waveguide furnishes an interesting example. Here the coupling 
coefficients are complex, so that (C-3) is not valid; further the powers 
in the various modes are not orthogonal, so that (C-4) is untrue. In this 
appendix we consider cases where a ^ but where (C-3) holds true, 
so that c(z) is pure real. We define a quantity P(z) by (C-4); however 
only for ideal metallic guide, where the a's are equal to zero, are we 
assured that P(z) really represents the total power. If the a's are not 
zero we have no reason to think that P(z) should be the total power; 
however the results given below render this plausible when c(z) is real. 
From (C-4) we write 



p'( z) = d £ i ^l = Uz)h*'(z) + h'(z)h*(z) 

(lz (C-5) 

+ h{z)h*'(z) + U{z)h*{z). 

Substituting (C-l) into (C-5), making use of (C-2) and (C-3), we find 

P'( Z ) = -2a„ | h(z) | 2 - 2a x | h(z) | 2 . (C-6) 

If P{z) is (he total power flowing along the guide, (C-6) has a simple 
physical interpretation. It says that each mode contributes to the de- 
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crease in power along the guide in proportion to the product of its at- 
tenuation constant and the power it carries. Equation (C-6) may be 
extended to any number of modes via straightforward matrix techniques. 
Finally, consider the case where 

ai) = ai = a } Aa = ao — ai = 0. (C-7) 

From the transformation of (88)-(89), (C-4) becomes 

P( f ) = e - 2az [ | G (z) | 2 + | flk(f) | 2 ]. (C-8) 

Similarly from (C-6), (C-7), and (C-4), wc find 

P'(z) = -2aP(z), (C-9) 

which has the solution 

P{z) = «- , "P(0) = e~ 2az [ | 7„(0) T + | /i(0) | 2 ]. (C-10) 

Assuming as usual that the guide is excited by a unit TEoi wave so that 
the initial conditions of (87) apply, i.e., 

7o(0) = 1, /i(0) - 0, (C-ll) 

(C-10) becomes 

P(z) = e" la \ (C-12) 

From (C-8) and (C-12) we have finally 

|Go(z)| 2 + |<?i(*)| 2 = 1, (C- 13 ) 

subject of course to the following conditions: 

(?o(0) = 1, ^(0) = 0. (CM) 

ao = ai , Aa = a u — ai — 0. (C-15) 

A similar treatment may of course be given for the coupled line equa- 
tions of (93). 

appendix d (See Section 2.3.8) 

Coupling Coefficients H [nm ] for General Continuous Deformations, and Beat 
Wavelengths B [nm] + , for Metallic Guide (Guide Diameter = 2 Inches) 
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n 


m 


Frequency = 50 Kmc 


Frequency = 55 Kmc 


Frequency = 60 Kmc 


Bi«i+ 

feet 


inch""' 


Bin.]* 
feet 


inch * 


fl(, lm ] + 
feet 


inch -1 




2 O.78861 


1.033665 


O.87IO!* 


0.935818 


0.95318 


O.855165 


3 O.30028 


1.5311*89 


0.33298 


1 . 380849 


0.36543 


1.25&058 


4 O.15879 


2.072235 


0.17717 


1.855771 


0.19531 


1.682663 


5 0.0961(1 


2.6891&0 


0.10863 


2 . 3808H 


0.12055 


2.142344 


6 0.06272 


3. 1*52631* 


0.07182 


2.992239 


0.08053 


2.658456 


7 0. 04176 


4. 574033 


0.01*936 


3-771644 


0.05630 


3.271486 


8 0.02590 


7.815578 


0.031*05 


4.986672 


0.04021 


4.077331 


9 CUT OFF 


0.02066 


10.765911 


0.02842 


5 . 422487 


1 1 -2.45264 


0.221*61+8 


-2.70093 


0.203996 


-2.94898 


0.186837 


2 1.99720 


0.56ml* 


2.20291 


O.508716 


2.40812 


0.465363 


3 0. 46389 


0.903729 


0.51323 


O.816787 


0.56232 


0.745451 


1* O.21503 


1.268585 


0.23908 


1.140518 


0.26289 


1.036959 


5 0.12293 


1.673662 


0.13773 


1.491619 


0.15228 


1.348056 


6 0.07762 


2.151892 


0.08805 


I.888713 


O.O9815 


1.690415 


7 0.051311- 


2.781*072 


0.05953 


2.368870 


0.06723 


2.084309 


8 0.03377 


:-. oiq- .-.>:■ 


0.04119 


3.029237 


0.04763 


2.572121 


9 CUT OFF 


0.02776 


4.333118 


0.03402 


3.272356 


10 CUT OFF 


CUT 


OFF 


0.02321 


4.869801 


2 1 -5.16225 


0.1*11*888 - 


5.6869'* 


0.376608 


-6.21100 


0.344833 


2 0.90059 


0.730808 


0.991*1*1* 


0.661825 


1.08797 


. 604918 


3 O.31537 


1.0811*01 


0.34959 


0.975358 


0.38358 


0.888839 


4 0.16311* 


L. 1+626 3b 


0.18196 


1.310395 


0.20053 


1.188501 


5 O.09816 


L.8972l»2 


0.11054 


1.680662 


0.12263 


1.512888 


6 0.06357 


2.1*33890 


0.C7271* 


2.III386 


O.O8153 


1 . 376880 


7 0.04225 


3.217515 


0.01*987 


2.659138 


O.05685 


2 . 308706 


8 0.02631 


5 . 387660 


0.03437 


3.507377 


0.04054 


2.874896 


9 CUT OFF 


0.02108 


7.134864 


0.02865 


3.812804 


1 9.28358 0.617958 10.23270 


O.560645 


11.18001 


0; 513138 


2 0.51*61*2 O.905720 


0.b0l+17 


0.819106 


O.66165 


0.747918 


3 0.2311*7 J 


L. 265707 


0.25719 


1.138773 


0.28266 


1.035917 


4 0.1281*9 1.667518 


O.14385 


1.487579 


O.I5895 


1.345276 


5 0.08006 I 


.11*1126 


0.09072 


I.881991 


0.10105 


I.685882 


6 O.O5262 I 


•763190 


o .06091 


2.357577 


O.O6872 


2.077242 


7 0.031*59 . 


.837927 


O.0 1*200 


3.006912 


0.04848 


2.560340 


8 CUT OFF 


3.02835 


4.252809 


0.03457 


3.248210 


9 CUT OFF 


CUT 


OFF 


0.02366 


4.759583 


J. 

1 2.01901* C 


.83421*6 


2.22697 


0.756350 


2.43441 


0.691899 


2 0.37568 1 


.087248 


D. 1*160 3 


0.981657 


O.45615 


0.895260 


3 0.17821 1 


.4581*01* 


3.19854 


1.308247 


O.21863 


1.187591 


4 0.10386 1 


.886313 


3. 11679 


I.67387O 


O.12945 


1.508455 


5 0.06627 2 


.1*1291*3 


3.07568 


2.099151 


0.08471 


1.868999 


6 0.01*377 3.169891 


3.05147 


2.036799 


O.05856 


2.295676 


7 0.02751* 5.053602 ( 


). 03538 


3.454117 


0.04157 


2.851181 


8 CUT OFF 0.02221 


6.258648 


0.02936 


3.751777 


9 




CUT OFF 


CUT OFF 


0.01672 


24.592342 
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H 


m 


Frequency = 50 Kmc 


Frequency = 55 Kmc 


Frequency = 60 Kmc 


B [„■.]* 

feet 


Hfnml 
inch"* 


B( 1 + 

feet 


E[„ml 
inch -1 


Bin*]* 
feet 


inch -2 



10 



1 


1.03171 


1.0631*96 


1.13893 


0.963362 


1.21*582 


0.880700 


2 


0.27706 


1.2765'*'* 


0.3071*0 


1 .150302 


0.33750 


1.047565 


3 


o.lirfU 


1.661453 


0.1581*2 


1.1*85036 


0.171*81* 


1. 31*1*734- 


It- 


0.08552 


2.123151 


0.0966^ 


1.871474 


0.10755 


1.679382 


5 


0.05538 


2.725027 


0.06388 


2.337156 


0.07193 


2 .061*780 


6 


0.03629 


3.721*150 


0.04371 


2.965282 


O.05029 


2.538314 


7 


CUT OFF 


0.02951* 


i+. 111601 


0.03571 


•3 .202782 


8 


CUT 


OFF 


CUT 


3FF 


.021*56 


4.570756 


l 


0.65318 


1.305633 


0.72180 


l.l8ll*5l* 


0.79015 


1 .079233 


2 


0.21370 


1.1*71*795 


0.23771 


1.325829 


0.26139 


1 .205391 


3 


0. 115*40 


1.87721*3 


0.1291*7 


1.070557 


0.14-326 


1.508273 


1* 


0.07137 


2.38393!* 


0.08123 


2.083255 


0.09071* 


1.859682 


5 


0.01*651 


3.0991*96 


0.051*37 


2.003525 


0.06168 


2.276663 


6 


0.0295** 


I+.66325I+ 


0.03715 


3.37'*330 


0.01*31*0 


2.614772 


7 


CUT 


OFF 


0.02386 


5.1*9238!* 


0.03060 


3.660164 


8 


CUT 


OFF 


CUT 


OFF 


0.01927 


7.278403 


1 


0. 1*5817 


1.56081*5 


0.50693 


1.1*10610 


0.5551*3 


I.287362 


2 


0.17021* 


1.683378 


0.18977 


1.509104 


0.20906 


! .369320 


3 


0.09560 


2.108882 


0.10773 


1.866529 


0.11957 


1.679271 


1* 


0.06011 


2.677695 


O.O6899 


2.313116 


0.0771*6 


2.051409 


5 


0.03903 


3.582532 


0.01*652 


2.9101*15 


0.05327 


2.509531 


o 


0.021 48 
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APPENDIX E 

Geometry of Continuous Bends 

Let x(z) and y(z) be the transverse displacements of the guide axis 
from the z-axis, in a rectangular co-ordinate system. Then we write 

r(«) = u-(z) + yj(z) + kz, (E-l) 

where i, j, and k are unit vectors along the x, y, and z axes respectively. 
Assuming that (215) holds true, from (216a) we have 

s^z. (E-2) 
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Thus, arc length along the bent guide axis is approximately equal to 
distance measured along the z-axis. From (211a) and (E-2), 

t=|«|=b'(*)+J2/'(*)+k, (E-3) 

where primes indicate differentiation with respect to z. From (211b) 

}p-^!~i*»M+te»M. (E-4) 

Therefore the curvature is approximately 

- » \A" 2 (*) + y"*(z) (E-5) 

p 



so that 



Then from (211c) 



Lr»(z) + ]ij"(z) 
P ~ V^(z) + y--(z) • (E " 6) 



b = tX P ^kXp^ W : W . (E-7) 

Vp 

Since p and b are approximately transverse, i.e., their z-components arc 
small, from (211) and (212) 

*«Zb-|«Zp, (E-8) 

since by (211c) b _|_ p. From (E-4) 

tan*«^. (E-9) 

X \Z) 

Thus, 

x"(z) 

"" 8 ^«) + a.) ^"' ( ' ) ' (E " 10) 

""" ^S^ *"^ (E - n) 

From (E-10) and (E-ll) we finally obtain the approximation of (210); 

— »»*(*), (E-12) 
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™i «¥*<•). (E-13) 

P 

The above analysis has been a little crude; it, is helpful to get a pre- 
cise estimate on the error in the approximate result of (E-4) for the 

vector -p, which in the above approximation is assumed to be purely 

P 
transverse. We write 

r(z) = \x{z) + ]y{z) + kz (E-14) 

as before. Next, 

ds = \/dx - + dif + dz\ (E-15) 

j = Vl + .v'Hz) + y' 2 (z) , (E-16) 

dz 

where as before we reserve the prime to denote differentiation with 
respect to z. Then 

_ dr _ dx dz 

-ds-Jzds (E17) 

VI + x n (z) + V 2 (z) 
Further 

1 rft dXdz 1 



?' = s = ss - r+?w+?w'f w + j ' y " w (E 18) 

Now we may write (E-18) as follows: 

i-p = ix"{z) + iu"{z) - A - B. (E-19) 

P 

A = i f ( £lf ( ^ S [k " (z) + ' VJ " {z)] - (E " 20) 

[1 + ;r' 2 (z) + (/'-(z)J- 

The first two terms of (E-19) are identical to the approximation of 
(E-4); the vectors A and B represent correction terms that we shall 
show to be small compared to the first two terms. 
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From (E-20) 

v^irhm - ifXtffw s X " M + «" M - (E - 22) 

From (E-21) 

, . | x'(z)x"(z) + !i'(z)y"(z) | r , 2 , . , 2 . 

|B!= [1 +*«(«)+ y«(*)P [1+3; U)+// W1 

^ !*'(*)«•(*) + y'(*V(*)| (E_23) 

^ W 2 (2) + // 2 (*) vV 2 U) + </" 2 (*), 
whore the lust stop follows from the Schwarz inequality. Then 

If (215) is satisfied, (E-22) and (E-24) show that the correction terms 
A and B of (E-19) are small compared to the first two terms, so that 
the approximation of (E-4) will he valid. 

APPENDIX F 

Rigorous Treatment of TE M Loss Statistics for the Discrete Case 

In treating the TEoi loss as a Fourier series with random coefficients, 
the frequency dependence of the Aa's and the C"s in (218) was neglected, 
since the principal frequency dependence occurs through the A/3's. While 
this provides a simple and accurate analysis, a rigorous treatment of 
the TEoi loss statistics as a function of frequency is of interest. 

We consider only the case of independent offsets or tilts, treated in 
Sections 3.3. 1 and 3.3.2 respectively. From (218), (23.5), and (236) 

A = A + 5,1 (F-l) 

A = i-l (F-2) 

it— i 
8 A = £ -U cos kAfik (F-3) 



,1, - e * A « J « £ ( x ? Xi+k n + Xi x x i+k x ) (F-4) 

</!„> = x % N, <.<!„-> = x*N(N + 1) (F-5a) 



..4 

X 



<A k > = 0, ol,-> = - (N - k) e k2Aa '°; l^k^N-l. (F-5b) 
— 
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<A k A,> = 0, k * I. (F-5c) 

The new quantity A is denned by (F-2) ; 54 remains the same as before. 
We have for the expected values of the various losses: 

<!> = §o4 >, <8A> = 0. (F-6) 

Therefore 

<A> = <A> = %<A > = | N. (F-7) 

Equation (F-7) gives the expected value of the TE i loss as a func- 
tion of frequency, since the rms conversion coefficient :t will in general 
vary with frequency (this variation will be small for offsets and diam- 
eter changes, approximately inversely with the free-space wavelength 
X for tilts). We may now average over wavelength (indicated by a bar) 
and obtain instead of the result given in (237): 

<!> = <I>=f AT. (F-8) 

Thus, the average of x over the band should be used instead of the value 
of # 2 at the middle of the band. Over reasonable bandwidths the error 
will be small. 

Next from (F-5c), A and 8 A are easily shown to be uncorrelated. 

<A(8A)> = 0. (F-9) 

We next find the mean square value of 8A. From (F-3), (F-5b), and 
(F-5c): 

AT— 1 N-l 

<(8A ) 2 > = £ Z) <A k A l > cos kAfik cos ZA/3/., 

A=l /=1 
N-l 

<(8A) 2 > = t T, (N - fc) c k ' 1Aa '° cos 2 fcA^o . (F-10) 

2 k=i 

The summation of (F-10) is easily written in closed form, but its gen- 
eral behavior is much more easily seen by examining the usual two 
special cases, small and large differential loss. 

1 . Small differential loss over total length L N = Nl Q . 

<(M)!> = |[ JV(JV _ 2) + (« ? y] ; _ N2Aah<<L mi) 
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2. Large differential loss over total length L N , small differential loss 
over pipe length l„ . 

j 1 / sin Agjo V 

<(*A)*> = f " /"fT/V ; - im * » X ' (F 19) 

4 — Aa/„ /sin Aj8/o\ (r-12j 



1 + 



— Aa/o 



-2Aa/ « 1. 



Sketches of the general behavior of <(5^4) 2 > vs A/3/ (proportional to 
the free-space wavelength X) are given in Fig. 11 for a single period, 
these functions being periodic of period -k. In both cases <(8A) 2 > is al- 
most constant except in narrow bands centered at A/3/ = mir, where 
it becomes twice as large; these peaks occur because we have assumed 
equally spaced mode converters. The half-width of these peaks initially 
decreases as 1.39/JV as long as the differential loss remains small, ap- 
proaching a limiting value of ( — Aa/ ) when the differential loss becomes 
large. Since these peaks are narrow, they may be neglected in averaging 
over X (or A/3/„), yielding for small and large differential loss respectively: 



<(5A) 2 > = <(5.4) 2 > = ^N(N - 2); -N2Aah« 1 (F-13) 



<(*)■> = <(IA)-> = (^ | ; -N2AaU » 1, (pi4) 

-2Aa/ « 1. 
Comparing with (241) and (242), we see that the average of x or of 

1* 

— — over the band should be used rather than their values at the middle 
— Aa 

of the band. Again, over reasonable bandwidths the error will be small. 

The minor difference between (F-13) and (241) for x 4 independent of 
frequency — i.e., the factor (N — 2) instead of (N — 1 ) — arises be- 
cause of our approximate integration of the function of Fig. 11(a), in 
which the narrow peaks and small ripples were ignored and the function 
set equal to N(N — 2). It is clear that this yields a result that is too 
small. An exact integration of (F-ll) yields a result identical to (241). 

The above calculations have considered only a single mode, the mode 
subscript being omitted as usual. For the total 8A, including all first 
order spurious modes, we have 

BA = T,8A lm]l (F-15) 

(,„] 

<(8A )'■'> = E Z <&i w M w >. (F-16) 
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<(<fA) 2 > 




N^-SN 



A/Jl ( 



<(<TA) 2 > 



(x 4 /4)Cn/-A«1 ) 



0.75 



0.50 



I 

-Aal 




A/3l c 



Fig. 11 — General behavior of <(SA) 2 > vs APU . 

In violation of our usual convention that [ ] indicate only TE spurious 
modes, in the present appendix this notation includes in addition the 
TMn + spurious mode in the case of tilts. 

The <(8A [m] ) 2 > have been discussed above; we need in addition the 
cross terms <8A [m] SA [lll >. These cross terms did not appear in the ap- 
proximate analysis, so we expect to find them negligible here. The differ- 
ent 8A lm] are of course not independent. From (218), (221) and (173) 
we see that 8A lm] (Xi) « 8A M (X 2 ) where \ t « x 2 ; i.e., the TE i loss 
component due to the m th spurious mode at one frequency is propor- 
tional to the TEo, loss component due to the n th spurious mode at a 
widely separated frequency. However, at the same frequency 8A [m] and 
SA[n] are almost uncorrelated, so that the cross terms in (F-16) may 
be neglected. 
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We have 

,V-1 N-\ 

<8A [m] SA[n}> = EE <A k[m] A llnJ > cos /••A/3 [ml / cos lA0 ln] / o (F-17) 
*-i 1=1 

whore; 

<A*^«.,> = ****** (N - k) ^*«tai+*«w»i; 

<,l, [ml .4 / |„)> = 0, A- ^ /. 



(F-18) 



Soiling 

cos /rA/3 ( „ ( |/ oos A-A/? ( , (1 /» = \ cos A-(A/3 (ml - A/? f „,)/ 

+ | cos/,-(A/3|, Ml + A/3,,,,)/,, 
(F-17) now becomes 

<5,l 1 „„5.4 [ „ l > = -Iwl^g ( N - fc) ^CA-w-M-uiM. ( 

•[cos A-(A/3,„,, - Aft Hl )Zo + cos />(A/3 [ml + Aft „,)/«.] 

This summation is of tho same general type as (F-10). In the special 
cases of small and large differential loss: 



A W U W > 8~L" 2iV + V^inl 

/ sin^V(Aff |ml +AgMJjp Y] . 
\ sini(^N) + A0 [Bl )/ o /J' 



I jVjAgw - A/3 [n] )/ „ 
i(A0, w| - Aft,,,)/,, 



(F-20) 



<54 tM ,SA [n ,> = 

1 - 



KAj8 w + A/3[ 



-iV(Aa lm] -f-Aa [n ])/„«l. 



.Y 



4 -(Aa [m] + Aa [((] )/ 

2 



— (At*!,,,] + Aa?[„] M» 



sin h( &&[„,] - A/3i„])/o]" 



?— -) [sin i(A/3 [m] - A/3 [n] )/„F 

r sin i(A0 tm] + A/3 [n] )/ ] 2 



1 + 



1 - 



-(A«, 



:f-2d 



i 



-(Aa [7 „| + Aa[„i)/ 



1 + 



[sin |(A(8 Im] + A/3 [B] )?o] 



— (Aa[ m] + Aa (M ])/ 
- AT(Aa( m] + Aa M )/„ » 1, -(Aa [ml + Aa [n) )/o « 1. 



<SA lm Mw> = - *™ A N > ~(*«M + Aa [nl )/„ « 1, 
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Equations (F-20) and (F-21) again exhibit narrow peaks of the same 
general type as illustrated in Fig. 1 1 ; away from these peaks the func- 
tions will be quite small. Thus if (A/3 [ffl) =F A0 [n] )/„ * mir, from (F-20) 
and (F-21) for either small or large differential loss: 

2 A 2 

N: -(Aa<„, + Aa<„,)L « 1. 

(F-22) 

(A/3 ( „ (] =F Afi[„])lo 9* mir. 

From (F-22) the correlation coefficient of 8A [m] and 8A [nl will be small 
for moderately large values of N. Therefore, for those modes that make 
a significant contribution to the total <(BA)b, the cross terms in (F-16) 
will be negligible. 

APPENDIX G 

Correlation Coefficient of TE i Loss Components due to Different Spurious 
Modes for the Continuous Case 

Consider the ac components 5.4 Im] and 5 A [n] of the total TEoi loss, 
due to two different spurious modes, each with two polarizations, gen- 
erated by the same type of geometric imperfection. The geometric im- 
perfection and thus the coupling coefficients are assumed to have white 

power spectra. 

As a specific example, consider random deviations of the guide axis 
from perfect straightness, which generate principally the forward TE I2 
and TEn spurious modes. The coupling coefficients to the two polariza- 
tions of the TE lm mode arc given in (308) in terms of the second deriva- 
tives of the rectangular co-ordinates of the guide axis, x"(z) and y"{z). 
We assume that x"(z) and y"{z) have white power spectra. 

We expand the geometric imperfection to which the coupling coeffi- 
cients are proportional in a Fourier series. For random straightness 
deviations : 

x"{z) = 2 yj e j2 ' n!lL ; 

,|=_0O 

>/(*)= E T. 1 .*"'! (CM> 

H=— 00 

The complex Fourier coefficients c n in Section 4.1 are simply proportional 
to the corresponding y„ . For straightness deviations, for the TEi ffl mode 



APPENDICES 



1169 



we have: 



C n [m] = Ctlm)y„ l , 
Cn[m) = Ct[ m ]'y n X , 



TE lm \ 



(G-2) 



Consequently, the 7„'s have the same statistical properties as the c„'s, 
given in Section 4.2. Since x"(z) and y"(z) have white power spectra, 



y n | > = 7 ■ 



<i7»' \> = < 

Further, the different 7's are strictly independent. 

We next define for convenience the following quantities: 



(G-3) 



VgJ(t) 1 




[aj "I 


g»Kt) 
gStt) 


= £ 

»=— CO 


A. 1 

J. 


_<//«)_ 




Jn\ 



(-1) 



„ sin ir(t — n) 
ir(t — n) 



t = 



A/3L 



L 

B 



(G-4) 



These quantities are proportional to the real or imaginary parts of / in 
(285b). These four quantities are approximately independent station- 
ary hand-limited Gaussian random processes in the practical case where 
the length L is large compared to the heat wavelength B, so that | t | » 1 . 
The autocorrelation function of each of these quantities, found directly 
from (G-4), is therefore 



r(r) ^ W A (*W ,,L (« + r)> = tf ^H. 



(CM) 



Next we define for convenience the quantity g 2 (t) as 

g\t) = gJ\t) + g^-(t) + g a ±2 (t) + g^Ht). (G-6) 

The autocorrelation function of g 2 {t) is easily found in terms of the 
autocorrelation function of the individual g's, given in (G-o). If x is a 
stationary Gaussian random process with autocorrelation R x (t) and 
y = x , then the autocorrelation of y, R u (t), is given by 37, 40 



Ry(r) = /?/(0) + 2R 2 (t), 



(G-7) 



where the first term corresponds to the dc component, the second to 
the ac component of y. Since the individual quantities whose squares 
appear on the right-hand side of (G-6) are independent random vari- 
ables, we have from (G-7) and (G-5) for the autocorrelation function 
R(t) oig 2 (t) 
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B(r) - <g\t)g\t + r)> = tf + 2? (^y, (G-8) 

where again the first term corresponds to the dc component, the second 
to the ac component of g (t). 

From the results of Section 4.1 we may now write the TE i loss due 
to the two polarizations of the m th spurious mode (TEi„, for straightness 
deviations) in terms of the function g 2 (t) denned in (G-6): 

. CtM'L* 2/. x , __ I A/3[„,] 1 L (Ci-f)) 

A [m ] = — -^ — g(h>»\), km] ^ — • v ; 

E(iuation (G-9) is appropriate for our present purposes because it places 
in evidence the relation between the different TE i loss components 
A [m] . From (G-9), (G-8), and the results of Section 4.1, the normalized 
correlation coefficient p„„, of 8A [m] and 8A [n] , the ac components of 
TEoi loss due to the m th and n th spurious modes, may be written in the 
following form: 

<8A [m] 8AM> 
p '"" = V<(8A [m] r> V<(5.4 |Hl ) 2 > 



/ sin ttt^ Y 
V m nm ) ' 



(|A/? [ffl] | ~ |A|9 lt ,i|)L 

Tmn = o 

Thus, 



(G-10) 



<(JLY. (cm: 

\irT mn / 



Since the different spurious modes have substantially different beat 
wavelengths (Appendix I)), p m „ « 1 for moderate values of length L. 
As a numerical example consider TE, 2 + and TE n + , for a total length 
between mode niters L = 200 feet. From (C-ll), p < 0.00036. In prac- 
tical cases the TE i loss contributions of the different spurious modes 
will be almost uncorrclatcd, so that <{8A) 2 > will be given by (322). 



